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Abstract 

In the first part of this Dissertation, I study the differences between LOCC (local operations and 
classical communication) and the more general class of separable operations. I show that the two 
classes coincide for the case of pure bipartite state input, and derive a set of important consequences. 
Using similar techniques I also generalize the no-cloning theorem when restricted to separable opera- 
tions and show that cloning becomes much more restrictive, by providing necessary (and sometimes 
sufficient) conditions. 

In the second part I investigate graph states and graph codes with carrier qudits of arbitrary 
dimensionality, and extend the notion of stabilizer to any dimension, not necessarily prime. I further 
study how and where information is located in the various subsets of the qudit carriers of arbitrary 
additive graph codes, and provide efficient techniques that can be used in deciding what types of 
information a subset contains. 
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Introduction and preliminary 
concepts 

1.1 Historical remarks 

Classical computers are indispensable in today's world and appear in almost all imaginable scenar- 
ios, ranging from simple MPS players to sophisticated supercomputers used in weather prediction. 
Although their architecture is strongly device-dependent, they all have one thing in common: every 
classical computer is a physical realization of a Turing machine, a mathematical model introduced 
by Alan Turing in 1937 |Tur37j which formalized the concept of computation. Intuitively, any com- 
puter with a certain minimum capability is, in principle, capable of performing the same tasks that 
any other computer can perform, if sufScient time and memory are provided. 

The fundamental processing unit of a classical computer is the bit: a two-state system commonly 
denoted by and 1. A computation is performed whenever some input bits arc processed by the 
computer, resulting in another sequence of bits that represent the output of the computation. Any 
two-state classical physical system can in principle represent a bit, and a classical computer able to 
operate on these bits and to produce a valid output can (in principle) be built using only classical 
devices, as billiard balls, pulleys and so on, although in real world applications everything tend to 
be miniaturized for efficiency purposes. 

A fruitful idea is that computation is physical: any physical system performs some kind of 
computation during its evolution. Consider for example a rock that is falling down from some 
height. The total falling time is directly proportional to the square root of the height, so in a sense 
by simply measuring the time one effectively computes the square root of the height. With a bit of 
imagination more sophisticated examples can be constructed. What if we make use of the intimate 
structure of quantum mechanics and quantum systems to perform computations? Are there any 
fundamentally new possibilities relative to the classical case, or, is a quantum computer potentially 
more powerful than a classical one? This idea was first introduced in 1982 by Richard Feynman 
[Fey82| , when the notion of a quantum computer was born. One may argue that a rock is also a 
(quite large) quantum mechanical system, so why should quantum mechanics be more "powerful" 
in performing computations than its classical counterpart? The main and fundamental difference 
between classical and quantum mechanics is that the latter has a much richer structure that allows 
for novel effects not present in classical physics due to its use of a Hilbert space. On larger scales 
the quantum effects tend to be smoothed out and the system becomes classical, or decoheres, but 
on smaller scales the effects can be quite significant. 

Although it was widely believed that a quantum computer can indeed be more powerful than a 
classical one, in cases such as simulation of complex quantum systems, the first quantum algorithm 



1. INTRODUCTION AND PRELIMINARY CONCEPTS 



2 



able to outperform any known classical version was discovered only in 1994 jSlio97j by the computer 
scientist Peter Shor. He invented a quantum algorithm able to factor large numbers in polynomial 
time and that did not (and still does not) have any efficient classical counterpart. 

Even though a genuinely good quantum algorithm existed, a great deal of skepticism was dis- 
played with respect to an actual physical realization. It was already known that quantum systems 
are very sensitive to external noise that induce errors in the computation, and the main problem 
seemed to be the inability to correct these errors. Classical error correction was well understood 
[MS77j . the basic idea behind the whole field being the usage of redundancy, or duplication of in- 
formation for better protection against errors. On the other hand, the "no-cloning" theorem of 
Wootters and Zurek [WZ82| . discovered in 1982, forbids the duplication of quantum information, 
e.g. non-orthogonal states cannot be duplicated. Therefore the perspectives for good quantum error 
correction schemes looked extremely unpromising. The solution was provided by the same Peter 
Shor in 1995 |Sho95| . who showed that good quantum error correction schemes, which are not simply 
based on duplication of quantum information, exist. Quantum error correcting codes arc of extreme 
importance in quantum information processing since they allow for high-fidelity transmission of 
quantum information and reduction of decoherence. 

A novel field of science. Quantum Computation and Quantum Information, suddenly became a 
hot topic at the intersection of physics, mathematics and computer science. It has developed rapidly 
since 1994 and remarkable theoretical as well as experimental progress has been achieved. Today the 
subject is still far from being well understood, and I hope this Dissertation contains some nontrivial 
contributions to it. 

1.2 Preliminary concepts 
1.2.1 Qubits and qudits 

The fundamental processing unit of a quantum computer is the qubit: a quantum system with two 
energy levels. The state of the qubit is often denoted by 



where a and /3 are complex coefficients satisfying |ap -I- |/3p = 1 and |0) and |1) are orthonormal 
basis vectors of a 2-dimensional complex Hilbert space Ji. The qudit is the natural generalization of 
the qubit to D level quantum systems and is represented by a complex Hilbert space of dimension 



1.2.2 Evolution and quantum channels 

The evolution of a quantum system interacting with an external environment is in general non- 
unitary. A standard way of describing such evolution uses the Kraus formalism: if the initial state 
of the system is given in terms of a density operator p, then after a time t the state evolves as 



a\0)+m, 



(1.1) 



D. 




(1.2) 



fe 



where the Fk are called Kraus operators that satisfy the closure condition 




(1.3) 



k 



where / denotes the identity operator. The Kraus representation (|1.2p can be formally derived by 
considering a unitary evolution of the combined system-environment, and then tracing away (or 
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measuring) the environment degrees of freedom. A non-trivial result is that any open quantum 
evolution can be represented in the form (jl.2[) . If the result of the measurement on the environment 
is known, e.g. equals k, then conditioned on this k the initial density operator of the system is 
transformed to 



P^ ^rt.. - (1-4) 



However, measurement is not a deterministic operation and the result k is obtained with some 
probability pk = TT:[FkpF^]. Therefore whenever the measurement results on the environment are 
not discarded, the initial density operator of the system is transformed to an ensemble 

If the system starts out in a pure state then (jl.5|) reduces to 

\'4>) ^ {Pk, l^fc)}, with FfelV;) - ^\iJk) and pk - {'^\FlF^\'ij:) . (1.6) 



Any evolution of the form ()1.2[) is also called a quantum channel, and ()1.2|) represents the 
Kraus representation of a quantum channel. Technically the map (|1.2|) is a completely positive 
trace preserving (CPTP) map, which intuitively means that it maps positive operators to positive 
operators, remains positive whenever the system is trivially enlarged to a larger one, and preserves 
the trace of an operator. The latter condition is imposed by the closure condition (|1.3p and ensures 
that probability is conserved. The study of quantum channels is the subject of intense theoretical 
investigations as they are much less understood than their classical counterpart. The interested 
reader can consult Chapter 8 of jNCOOj for a good introduction to the subject. 



1.2.3 Bipartite (multipartite) quantum systems 

In quantum theory a multipartite quantum system is described by a Hilbert space constructed 
as a tensor product of the individual Hilbert spaces, and any vector in this tensor product space 
represents a valid quantum state. This is one instance in which quantum mechanics has a much 
richer structure than classical mechanics, because the tensor product description allows the existence 
of entangled states, i.e. quantum states that cannot be written as a tensor product of individual 
states. For example, in a bipartite qubit quantum system represented by a Hilbert space "Ha ®'Hbi 
a state of the form 

\i;) = a\m)AB + fi\ll)AB (1.7) 

with |ap + = 1 is entangled as long as < |ap < 1. Otherwise it is called a product state. 
Whenever \a\ = \ = l/\/2 the state is called maximally entangled. Throughout this Dissertation 
I will use the shorthand notation |00)ab to denote the tensor product |0)^ ® \^)b- 

It turns out that any normalized bipartite entangled state € T-La ® T~Lb can be written in a 
canonical form known as the Schmidt form: there always exist ortho normal bases {|aj)} and {|6a;)} 
of T-La and T-Lbi respectively, in which the |^/;) has the form 

D-l 

IV) = VKW)A\br)B- (1.8) 

r=0 

Here D is the minimum of the dimensions of T-La and T-Lb- The Aj.'s can always be chosen to be 
positive real numbers that satisfy "^^.^r = 1 and are called Schmidt coefficients. An alternative 
definition is that a bipartite pure state is maximally entangled if and only if all Schmidt coefficients 
are equal. When all Schmidt coefficients except one are zero then the state is a product state. 
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otherwise it is partially entangled. The number of positive Schmidt coefRcients is called the Schmidt 
rank of \4>). If all Schmidt coefficients are strictly positive than we say that lip) has full Schmidt 
rank. There is no analog of the Schmidt form for multipartite pure states, and this constitutes a 
major obstacle in understanding them. 

Entanglement is often considered a precious resource, since it can be "consumed" by various 
non-classical protocols such as quantum teleportation [BBC"'"93] . quantum dense coding |BW92| 
etc. It also constitutes a key ingredient (although by no means the only one) in the construction 
of good quantum error codes and exponentially faster quantum algorithms, and therefore its study 
constitutes an important part of quantum information theory. 



1.2.4 LOCC and Separable Operations 

An important paradigm in quantum information is that of local operations and classical commu- 
nication (LOCC). Consider for example two spatially separated parties, traditionally named Alice 
and Bob, each having access to local quantum systems described by Hilbert spaces Ha and Hb, 
respectively. Both Alice and Bob are allowed to perform arbitrary quantum operations on their 
individual quantum systems, and can also make use of a classical channel to communicate. They 
are not allowed, however, to exchange quantum systems between themselves nor to use a quantum 
channel. What kind of tasks can be performed in this paradigm? What are the restrictions com- 
pared to the global case (i.e. when global quantum operations are allowed on the combined system 
T-La <8i 'Hb)'^ Understanding LOCC constitutes an extremely important research program, since in 
an actual realization of a quantum computer many qubits may be well separated in space, and per- 
forming a global operation on them, in contrast to LOCC, will often be challenging, at least from 
an experimental point of view. 

Every LOCC operation can be regarded as a composition of local operations conditioned on 
particular measurement results that may be communicated through the classical channel. It is not 
hard to see that the initial state p of a quantum system transforms under LOCC as 

p^Y.^Ak®Bk)p{Ak®Bk)\ (1.9) 

k 

with 



E 



{AIA^®bIb,,)=Ia®Ib. (1.10) 



k 



The concepts above generalize to more than two parties in a straightforward manner. 

One may be tempted to say that any quantum operation of the form (|1.9|) represents a valid 
LOCC, but this is not true! There are operations of this form that are not LOCC |BDF+99| . that 



is, cannot be implemented by an LOCC paradigm, and which are called separable operations. There 
exists a simple example |BDF+99] of a set of basis states in a bipartite qutrit [D = 3) system that 
cannot be distinguished by LOCC but can be distinguished by separable operations. Hence LOCC 
is a proper subset of this more general class of separable operations. Although separable operations 
are not always implementable by LOCC, studying them is worthwhile since in general they have 
a cleaner mathematical formulation than LOCC and any result proven to be true for all separable 
operations will automatically be valid for LOCC, since the latter is a subset of the former. 

In the bipartite setting, a maximally entangled state plays the role of a universal resource, i.e. 
from an operational point of view can be transformed deterministically by LOCC to any bipartite 
partially entangled state |Nie99| . This is not true anymore in the multipartite regime; there is no 
universal multipartite quantum state that can be transformed to any arbitrary multipartite state by 
LOCC |HHHH09| . 

Entanglement can be quantified by various measures, and the measure is called an entanglement 
monotone if it is non-increasing under LOCC. In general entanglement measures have a simple form 
only for pure bipartite states, and in this case depend only on the Schmidt coefficients of the state. 
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For multipartite pure states or even for bipartite mixed states such a Schmidt decomposition does 
not exist, and entanglement in these cases is far from being understood. For a good introduction to 
the theory of entanglement see ^HHHH09) . 

1.2.5 Graph states and graph codes 

Quantum states and quantum entanglement are much better understood in the bipartite setting than 
in the multipartite setting. Two main difficulties that arise in the latter case are: i) the exponential 
growth, with the number of constituent parts, of the number of complex amplitudes used to describe 
a multipartite quantum state, and ii) the non-existence of a Schmidt representation. 

However, graph states form a class of multipartite states with a fairly simple structure. Given a 
graph G = (V, E) with n vertices V , each corresponding to a qubit, and a collection E of undirected 
edges connecting pairs of distinct vertices (no self loops are allowed) , a graph state is obtained by 
preparing a set of initial qubits in the 1+) = (|0) + |l))/-\/2 state, then applying controlled-phase 
gates between any two neighbors that are connected by an edge in the corresponding graph. Graph 
states can be generalized to higher dimensional qudits in a direct manner, using graphs with multiple 
edges, as described in detail in Chapters [S] and [HI 

Graph states were first introduced by Raussendorf and Briegel |RB01| in their measurement-based 
computation morfe/ (often called the one-way model). This new model of quantum computation 
is fundamentally different from the well-known circuit model, and is also universal. Any desired 
quantum computation can be achieved by only performing local measurements (in a basis that is 
conditioned on previous measurement results) on the qubits of a sufficiently large cluster state, 
a graph state in which the graph is a finite part of a lattice such as the square lattice. Graph 
states are therefore universal resources for quantum computation and intense theoretical as well 
as experimental work has been dedicated to their study. For a comprehensive introduction to the 
subject see |HDE+] . 

Graph states are also an instance of the so-called stabilizer states |Gotb| . A stabilizer state is 
a multipartite quantum state that can be described by an Abelian group (the stabilizer group) of 
Pauli-like operators on the Hilbert space of the carriers, each of which leaves the state invariant. 
Instead of describing a stabilizer state of n qubits by 2" complex amplitudes, it is enough to specify 
the generators of its stabilizer group of which there are no more than n. Hence stabilizer states allow 
for a very compact description. Stabilizer states play an extremely important role in the theory of 
quantum error-correction codes and they extend the notion of linear classical error correcting codes 
[MS77j to the quantum domain. 

1.2.6 Location of information 

First let us define what we mean by (classical) information. Information is embedded in correlations 
between two systems, e.g. we say that information about a system A is "located" in another system 
B if the statistical correlations between A and B are such that information about some properties 
of A can be recovered from B. For example, the photons that bounce off the Sunday newspaper hit 
the reader's retina and correlations between the letters on the newspaper page and reader's brain 
are established: the information about the latest news is now located in the reader's brain. If the 
systems are perfectly correlated then we say that all information about A is perfectly present in _B, 
and if they are totally uncorrelated than we say that no information about A is present in B (or, 
equivalently, that all information about A is absent from B). Traditionally A is considered to be 
the input of a communication channel and B its output at a later time, but information theory is 
not restricted to channels. 

Whereas the concept seems to be quite natural and simple, the rigorous mathematical theory 
of information was founded by Claude Shannon only in 1948 |Sha48| and has been continuously 
developed since then (see |CT05] for a comprehensive introduction). Information theory is a very 
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important subject in modern communication, cryptography, classical error-correcting codes etc, with 
applications ranging from audio CD error-correction to military satellite communication. 

Unlike classical information, quantum information can be present in more than one "type" , a 
terminology introduced by Griffiths |Gri07| . and various types can be incompatible, e.g. associated 
to operators that do not commute. Formally, a type of information is associated with a projective 
decomposition of the identity 



We also associate a type of information with a normal operator through its spectral decomposition, 
and, for example, call the information corresponding to the x component of an angular momentum 
of a spin one-half particle, represented by the Pauli operator a^, usually denoted by X, the X-type 
of information, and the information corresponding to the z component as the Z-type of information. 
Of course the X and Z types are incompatible, since their corresponding operators do not commute. 

1.3 Overview of the Dissertation 

1.3.1 Separable operations 

An LOCC acting on a pure bipartite state \iP)ab will in general produce an ensemble of states 
{Pki \4'k)AB}- where is the probability of obtaining the result k through a measurement of the 
environment. Finding necessary and sufficient conditions for when such a transformation is possible 
represents an important problem, and a complete solution for bipartite systems was first provided 
by Nielsen for an ensemble with only one output state |Nie99| and generalized by Jonathan and 
Plenio to ensembles with a finite number of states jJP99| . Both of these necessary and sufficient 
conditions are given in terms of majorization relations [Bha97| between Schmidt coefficients, and 
they completely characterize LOCC operations acting on pure bipartite states. 

In Chapter[5]and Chapter[3]we study separable operations acting on a pure bipartite state, trying 
to generalize previously known results. The most important result is that any output ensemble 
produced by a separable operation acting on a pure bipartite state can in fact be produced by some 
LOCC acting on the same state. Our result effectively says that LOCC and separable operations 
are the same class of quantum operations when acting on pure bipartite states. In particular, we 
prove that the majorization conditions of Jonathan and Plenio jJP99| are necessary and sufficient 
in the more general case of separable operations. 

Our result also has important consequences in the theory of entanglement, implying that a large 
number of mixed-state entanglement monotones remain monotone under separable operations. Since 
such monotonicity under LOCC has long been considered a necessary, or at least a very desirable 
condition for any "reasonable" entanglement monotone, one wonders whether monotonicity under 
separable operations, in principle a stronger condition, might be an equally good or even superior 
desideratum. 

An interesting question that follows from our result is: are separable operations and LOCC the 
same class of quantum operations when acting on multipartite pure states? It might be, but proving 
it would require very different methods than the ones we used, since there are no simple analogs 
of the Schmidt decomposition and majorization conditions. Necessary and sufficient conditions are 
not known even for LOCC. 

1.3.2 Local cloning of bipartite entangled states 

In Chapter |4] we consider the slightly different problem of cloning orthogonal entangled states by 
separable operations, a problem that belongs to the more general framework of deterministic mapping 
of an ensemble of pure states (and not just one state, as before) into another ensemble of pure states 




(1.11) 



3 
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by a separable operation. As summarized by the "no-cloning" theorem of [WZ82j , any set of quantum 
states can be deterministically cloned if and only if the states in the set are mutually orthogonal. 
When the states are not orthogonal, there is no deterministic apparatus capable of performing such 
a cloning. However, probabilistic cloning may still be possible and a significant amount of work has 
been dedicated to studying this case |SIGA05] . 

Formally, a set of quantum states is cloned whenever there exist a quantum operation 

that performs 

1^,) ® 10) 1^,) ® IV,), yj. (1.12) 

If the transformation is deterministic for all j, then the cloning is deterministic. Otherwise is 
probabilistic. The state |</)) plays the role of a resource or a "blank state" in which the copy of \tjjj) 
is to be imprinted. 

When the set consists of bipartite entangled states, and the cloning is restricted to LOCC (or to 
separable operations), the problem becomes much more difficult, and further restrictions have to be 
imposed. The mere orthogonality of the states no longer implies that they can be (locally) cloned. 
An LOCC analog of the "no-cloning" theorem was not yet found, and finding it may prove useful. 

It turns out that any two (and no more than two) orthogonal maximally entangled two-qubit 
states can be locally cloned by LOCC, using a maximally entangled "blank state" |ACP04| (on 
which the copy is to be imprinted) . A generalization to D maximally entangled states of two qudits 
of prime dimension D was given in |OH06j , which showed that a set of D such states can be locally 
cloned using a maximally entangled resource if and only if the states in the set are locally (cyclically) 
shifted 

r=0 

where the © symbol denotes addition modulo D. Kay and Ericsson |KE06j extended the above 
results to the LOCC cloning of full Schmidt rank partially entangled states using a maximally 
entangled blank state. They presented an explicit protocol for the local cloning of a set oi D x D 
cyclically shifted partially entangled states of the form 

D-l 

|V;.) = ^ VA;|r)-^|r®i)^ (1.14) 

r=0 

using a maximally entangled blank state, but failed to prove that any clonable set of states must be 
of this form. 

We investigate the conditions under which a set of pure bipartite quantum states on sl D x D 
system can be locally cloned deterministically by separable operations when at least one of the states 
is full Schmidt rank. We do not assume that D is necessarily a prime number and we also allow for 
the possibility of cloning using a resource state that is less than maximally entangled. We derive a 
set of necessary conditions, that are also sufficient in the case of qubits. In this latter case we proved 
a long-standing conjecture that a maximally entangled state is a necessary resource for such local 
cloning, even if the states to be cloned are partially entangled. We also generalize the protocol of 
Kay and Ericsson and show that any set of partially entangled "group-shifted" states 

1^/) = E V^l5)^l/.9)'', (1-15) 

where G is a group of order D and the elements of the group label an orthonormal basis of the 
(local) Hilbert space, can be locally cloned using a maximally entangled blank state, by providing 
an explicit LOCC circuit. Our protocol reduces to the one of Kay and Ericsson in the case of cyclic 
groups, since the latter is isomorphic to the additive group of integers mod D and the set of states 
defined by (|1.14p and (jl.lSp are the same. 

Our results significantly extend previous work in the literature (limited only to LOCC) jGKR04[ 
IACP041 [OH061 IKE061 ICEKROT] . 
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1.3.3 Graph codes 

Using graph states provides a fruitful approach for constructing good quantum error correcting codes, 
called graph codes. A graph code is a subspace of the Hilbert space of n carrier qubits spanned by 
a collection of graph basis states: a quantum state obtained from a fixed graph state by applying 
local Z operators on some of the qubits. See Chapter [5] for a detailed introduction to graph codes. 

There have been extensions to higher dimensional qudits, but all of them considered only qudits 
of prime dimensionality D. The main difficulty in the non-prime case is that Z/j is a ring, not a 
finite field, and the lack of the multiplicative inverse operation poses some technical problems, which 
we have successfully solved. 

Most known quantum error correcting codes are graph codes or arc equivalent to graph codes 
under local unitary transformations, hence understanding them is extremely important. In Chapter[5] 
we present an elegant method for constructing such graph codes, allowing for carrier qudits of 
arbitrary dimensionality, not necessarily prime. Our method allows for additive (or stabilizer |Gotb| ) 
as well as non-additive graph codes and was simultaneously developed by Cross et al |CSSZ09| for 
qubits and Chen et al [CZCOSj for higher dimensional qudits. We use simple graphical methods and 
computer searches to construct both additive and non-additive quantum error correcting codes, but 
computer searches are much faster for additive codes. In a number of cases we have been able to 
construct what we call quantum Singleton codes that saturate the quantum Singleton bound |KL97| . 
Our numerical techniques are based on finding a maximum clique |GJ79j in some related graph. The 
maximum clique problem on a general graph is known to be NP-complete, but our method may still 
be useful for constructing quantum codes with relatively small number of carriers. We also generalize 
the concept of a stabilizer group to the non-prime case and derive an elegant duality (known before 
only in the prime case) between the coding space and its corresponding stabilizer group. 

1.3.4 Location of quantum information 

In Chapter [H] we develop a mathematical formalism that can be successfully applied in studying how 
quantum information is encoded in additive graph codes and where is it located, a subject closely 
related to the one of Chapter[5l Studying additive graph codes is worthwhile since the vast majority 
of quantum error correcting codes are stabilizer codes, and stabilizer codes are locally equivalent to 
additive graph codes [Schb| . 

We show how to encode some input quantum information in the carrier qudits of an additive 
graph code and demonstrate how to use the concept of types of information to study the location of 
quantum information in arbitrary subsets of the carrier qudits. What types and how much informa- 
tion about the input can be then recovered? To various types of information we associate a collection 
of operators on the coding space which form what we call the information group. It represents the 
input information through an encoding operation constructed as an explicit quantum circuit (hence 
generalizing the encoding methods developed before only for prime dimensional qudits). Our for- 
malism is very general and works for arbitrary additive graph codes of arbitrary dimension (not 
necessarily prime). We also present an efficient numerical algorithm that can be successfully used 
in deciding where information is located and which types of information are present. As a side 
remark, note that we have not studied the "recovery problem" , i.e. finding the decoding operation 
that effectively "extracts" the quantum information from some carrier qudits, but in principle such 
a decoding always exist, provided all quantum information is located in these qudits. This recovery 
operation is interesting, but is not included in this Dissertation. 

The methods presented here allow for a better understanding of the intimate nature of quantum 
codes and may be of use in constructing better quantum error-correcting codes or quantum secret 
sharing schemes [MS08| . 
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1.3.5 Bipartite equientangled bases 

Chapter [7] is not closely related to the other chapters but presents a solution to a problem posed 
in jKM06j of constructing a family of "equientangled bases" for a bipartite system of two qudits 
of arbitrary (but equal) dimension: (i) The basis continuously changes from a product basis to a 
maximally entangled basis, by varying a parameter t, and (ii) for a fixed t, all basis states are equally 
entangled. 

We actually construct two solutions to the problem, one based on quadratic Gauss sums and the 
other using qudit graph states. These bases may find applications in various quantum information 
protocols including quantum cryptography, optimal Bell tests, investigation of the enhancement of 
channel capacity due to entanglement and the study of multipartite entanglement. 

1.4 The structure of the Dissertation 

All chapters of this Dissertation arc self contained and consist of published (or accepted for publi- 
cation) articles in refereed journals. The contents of each chapter is almost the same as that of the 
published paper, with minor modifications made for the sake of consistency of notation throughout 
the Dissertation. Most of the chapters represent collaborative work with different persons in our 
research group, as described below. 

• Chapter [2J Entanglement transformations using separable operations: published in 
Physical Review A (GG07| . Collaboration with Robert B. GrifRths. Both authors made major 
contributions. My most important contributions were: i) the application of map-state duality 
formalism ii) the idea of using an inequality by Minkowski in deriving necessary conditions 
for pure state transformations under separable operations, and iii) the investigation of random 
separable unitary channels. 

• Chapter [U Separable operations on pure states: published in Physical Review A as a 
Rapid Communication iGGOSi . Collaboration with Robert B. Griffiths. Both authors made 
major contributions. My most important contributions consisted in: i) extensive numerical 
studies that led us to the conjecture that separable operations are implemcntablc by LOCC in 
the case of pure bipartite states; ii) the application of map-state duality formalism; iii) parts 
of the proof of the main majorization theorem that was conjectured by Griffiths; iv) derivation 
of the consequences of our result, the most important being that all convex- roof mixed state 
entanglement measures remain monotone under the more general class of separable operations. 

• Chapter m Local cloning of entangled states by separable operations: accepted for 
publication in Physical Review A and available on arXiv |GYC| . Collaboration with Scott M. 
Cohen and Li Yu. I made major contributions to this work, including: i) the map-state duality 
formalism used in our investigation of the problem; ii) various necessary conditions, such as the 
necessary form of qubit entangled states, equality of G-concurrence, information-theoretical 
observation etc.; and iii) various proofs of theorems. Scott Cohen introduced the idea of using 
finite groups to study sets of clonable states, and Li Yu proved that a maximally entangled 
state is necessary for the local cloning of group shifted states in Z) = 2 and D ^ 3. 

• Chapter O Quantum error correcting codes using qudit graph states: published in 
Physical Review A [LYGG08] . Collaboration with Shiang Yong Looi, Li Yu and Robert B. 
Griffiths. This work was not one of my main projects. My main contributions consisted in 
developing the stabilizer formalism for non-prime qudits in Sec. 15.51 and proof of the X — Z 
rule for qudit graph states in Sec. 15. Al I was not involved in the numerical work for searching 
good quantum error-correcting codes. 
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• Chapter [SI Location of quantum information in additive graph codes: published in 
Physical Review A jGLGlOj . Collaboration with Shiang Yong Looi and Robert B. Griffiths. I 
had a major part in this paper, together with my co-authors. My most important contributions 
were: i) introducing a set of useful Clifford gates for arbitrary dimensions; ii) the use of 
Smith diagonal forms over rings of integers, which allows one to deal with problems that 
appear in the non-prime case; iii) the discovery of a general encoding operation in terms of an 
explicit quantum circuit that extended previous work restricted to qudits with prime D; iv) 
implementation of an efficient linear-algebra algorithm used to decide where and which types 
of informations are present in a given subset of the carriers; v) proofs of various theorems. 

• Chapter [71 Bipartite equientangled bases: accepted for publication in Physical Review A 
and available on arXiv [GLj . Collaboration with Shiang Yong Looi. This work was split into 
two parts, between myself and my co-author. I found the solution based on Gauss sums (first 
part), whereas my co-author found the one based on graph states (second part). I also studied 
the entanglement properties of the second solution in terms of G-concurrence, and proved its 
monotonicity as a function of t. 
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Entanglement transformations 
using separable operations 



2.1 Introduction 



A separable operation A on a bipartite quantum system is a transformation of the form 



p' = Hp) = I] (An <^ B^)piAl <E> Bl), 



(2.1) 



where p is an initial density operator on the Hilbert space Ha ^'Hb- The Kraus operators A^ ® B, 
are arbitrary product operators satisfying the closure condition 



The extension of (|2.ip and (|2.2p to multipartite systems is obvious, but here we will only consider 
the bipartite case. To avoid technical issues the sums in (|2.ip and (|2.2p and the dimensions of Ha 
and Hb are assumed to be finite. 

Various kinds of separable operations play important roles in quantum information theory. When 
m takes on only one value the operators Ai and Bi are (or can be chosen to be) unitary operators, 
and the operation is a local unitary transformation. When every Am and every _B,„ is proportional 
to a unitary operator, we call the operation a separable random unitary channel. Both of these are 
members of the well-studied class of local operations with classical communication (LOCC), which 
can be thought of as an operation carried out by Alice on Ha with the outcome communicated 
to Bob. He then uses this information to choose an operation that is carried out on Hb, with 
outcome communicated to Alice, who uses it to determine the next operation on Hai and so forth. 
For a precise definition and a discussion, see [ |HHHH09) . Sec. XI]. While any LOCC is a separable 
operation, i.e., can be written in the form (|2.ip . the reverse is not true: there are separable operations 
which faU outside the LOCC class |BDF+99| . 

Studying properties of general separable operations seems worthwhile because any results ob- 
tained this way then apply to the LOCC subcategory, which is harder to characterize from a math- 
ematical point of view. However, relatively little is known about separable operations, whereas 
LOCC has been the subject of intensive studies, with many important results. For example, an 
LOCC apphed to a pure entangled state (i.e., p = \ip) {ip\ in (|2.ip ) results in an ensemble of pure 
states (labeled by to) whose average entanglement cannot exceed that of [ [HHHH09] . Sec. XV 
D]. One suspects that the same is true of a general separable operation A, but this has not been 
proved. All that seems to be known is that A cannot "generate" entanglement when applied to a 




(2.2) 
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product pure state or a separable mixed state: the outcome (as is easily cheeked) will be a separable 
state. 

If an LOCC is applied to a pure (entangled) state \^), Lo and Popescu |LP01| have shown that 
the same result, typically an ensemble, can be achieved using a different LOCC (depending both on 
the original operation and on \tJj)) in which Alice carries out an appropriate operation on Ha and 
Bob a unitary, depending on that outcome, on Hb- This in turn is the basis of a condition due 
to Nielsen |Nie99| which states that there is an LOCC operation deterministically (probability 1) 
mapping a given bipartite state to another pure state \(j>) if and only if \(j>) majorizes |^/;) 0. 

In this chapter we derive a necessary condition for a separable operation to deterministically 
map to 10) in terms of their Schmidt coefficients, the inequality (|2.5p . While it is weaker than 
Nielsen's condition (unless either T-Ia or "Hb is two dimensional, in which case it is equivalent), it is 
not trivial. In the particular case that the Schmidt coefficients are the same, i.e., \ip) and |0) are 
equivalent under local unitaries, we show that all the Am and Bm operators in (|2.ip are proportional 
to unitaries, so that in this case the separable operation is also a random unitary channel. For this 
situation we also study the conditions under which a whole collection of pure states are 

deterministically mapped to pure states, a problem which seems not to have been previously studied 
either for LOCC or for more general separable operations. 

The remainder of this chapter is organized as follows. Section 12.21 has the proof, based on a 
inequality by Minkowski, p. 482 of |H J99| . of the relationship between the Schmidt coefficients of 
\ip) and when a separable operation deterministically maps ji/') to |0), and some consequences 
of this result. In Section [2.31 we derive and discuss the conditions under which a separable random 
unitary channel will map a collection of pure states to pure states. A summary and some discussion 
of open questions will be foimd in Section \TM 

2.2 Local transformations of bipartite entangled states 

We use the term Schmidt coefficients for the nonnegative coefficients {Aj} in the Schmidt expansion 

D-l 

1^) = E (2.3) 

of a state |'0) € Ha^'Hb, using appropriately chosen orthonormal bases {1%)} and {|&j)}, with the 
order chosen so that 

Ao > Ai > • • • > Ai3_i > 0. (2.4) 

The number r of positive (nonzero) Schmidt coefficients is called the Schmidt rank. We call the 
subspace of T-La spanned by |ao), |ai) . . . \ar-i), i.e., the basis kets for which the Schmidt coefficients 
are positive, the Ha support of {ip), and that spanned by |6o), l^i) • ■ • \br-i) its Hb support. 
Our main result is the following: 

Theorem 2.1. Let \ip) and |(/)) be two bipartite entangled states on Ha ®Hb with positive Schmidt 
coefficients {Aj} and {/ij}, respectively, in decreasing order, and let r be the Schmidt rank of \^) . If 

can be transformed to by a deterministic separable operation, then 

i) The Schmidt rank of \(j}) is less than or equal to r. 

ii) 

1 — 1 r— 1 

J]A,>nM.- {2-5) 

J=0 3=0 

^ By "!</)) majorizes \^)" we mean that the vector of eigenvalues of the reduced density operator p((p) of \(p) on Ha 
majorizes that of the reduced density operator pitp) of in the sense discussed in INie99l . or in f INCOOI . Sec. 12.5.1]: 
the sum of the k largest eigenvalues of p(0) is never smaller than the corresponding sum for p{ip). A helpful discussion 
of majorization is also found in [HJ99| (see the index), with, however, the opposite convention from Nielsen for "A 
majorizes B" 
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in) If (j2.5p is an equality with both sides positive, the Schmidt coefficients of \ip) and \(f)) are 
identical, Xj = fij, and the operators Am and Bm restricted to the "Ha o,nd Hb supports of \Tp), 
respectively, are proportional to unitary operators. 

iv) The reverse deterministic transformation of\(j)) to ji/;) by a separable operation is only possible 
when the Schmidt coefficients are identical, Xj ~ fij . 

Proof. For the proof it is convenient to use map-state duality (see |ZB04[ [GWYC06] and [ [BZ06j . 
Chap. 11]) defined in the following way. Let b e an orthonormal basis of Hb that will remain 

fixed throughout the following discussion. Any ket \x) G "Ha ®'Hb can be expanded in this basis in 
the form 

|X>-EI"^)®I^^)' (2.6) 
j 

where the {Ictj}} are the (unnormalized) expansion coefficients. We define the corresponding dual 
map X ■'Hb ^ Ha to be 

X = J2\a,){b,\. (2.7) 

j 

Obviously any map from "Hb to Ha can be written in the form (j2.7p . and can thus be transformed 
into a ket on 'Ha®'H b by the inverse process: replacing {bj \ with \bj). The transformation depends on 
the choice of basis but this will not matter, because our results will in the end be independent 

of this choice. Note in particular that the rank of the operator x is exactly the same as the Schmidt 
rank of |x)- 

For a separable operation that deterministically maps \^) to \(f>) (or, to be more specific, 
to 10) (01) it must be the case that 

(An®S„)|^) = ^10), (2.8) 

for every m, as otherwise the result of the separable operation acting on would be a mixed 
state. (One could also include a complex phase factor depending on to, but this can be removed 
by incorporating it in Am — an operation is not changed if the Kraus operators are multiplied by 
phases.) By using map-state duality we may rewrite (|2.8p in the form 

AmlpBm = y/p^(f>, (2.9) 

where by B^ we mean the transpose of this operator in the basis — or, to be more precise, 

the operator whose matrix in this basis is the transpose of the matrix of Bm- From l|2.9p one sees 
at once that since the rank of a product of operators cannot be larger than the rank of any of the 
factors, the rank of cj) cannot be greater than that of 0. When translated back into Schmidt ranks 
this proves (i). 

For the next part of the proof let us first assume that Ha and He have the same dimension D, 
and that the Schmidt ranks of both \ip) and \(f)) are equal to D; we leave until later the modifications 
necessary when these conditions are not satisfied. In light of the previous discussion of (|2.9p . we sec 
that Bm has rank D, so is invertible. Therefore one can solve (|2.9p for Am, and if the solution is 
inserted in (|2.2|) the result is 

/® / = Y,Pmy^-'^Bm'\<l>U)Bm'^~'] ® [^Jn^m] (2-10) 
m 

The Minkowski inequality f jHJ99| . p. 482) for a sum of positive semidefinite operators on a 
5-dimensional space is 

det(EQ„) >E(detQ,„) , (2.11) 
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with equality if and only if all Qm^s arc proportional, i.e. Qi ~ fijQj, where the fij are positive 
constants. Since Al^A^ (g) Bl-^Bm is a positive operator on a S* = dimensional space, (|2.10p and 
yield 



1 > 



> 



E 



det (p„[V'"'^S,-it(0t^)s^-i^-i] ^ [BlB^] 

det(0V)i/^ _ det((/.t(/.)V^ 
"det(V^V^)i/-° ^ det(^/'t^)i/i3' 



1/D^ 



which is equivalent to 



det('0^-0) > det((/)1' 



(2.12) 



(2.13) 



The relation det(j4 ® _B) = (det 74)''(det , where a, b are the dimensions of A and B, was used in 
deriving (|2.12p . Since p.l3p is the square of p.Sp . this proves part (ii). 

If (|2.5p is an equality with both sides positive, det{(f>^ (f>) / dot {il'^ip) = 1 and the inequality (j2.12p 
becomes an equality, which implies that all positive operators in (|2.1ip are proportional, i.e. 



AlAn ® BlB.^ = IranAiAn ® B^B^, (2.14) 

where the /m„ arc positive constants. Setting n = 1 in (I2.14p and inserting it in (|2.2p one gets 



{Y,fr.^l)A\Al®BlBl^I®I. 



(2.15) 



This implies that both and b\b^ are proportional to the identity, so Ai and Bi are pro- 

portional to unitary operators, and of course the same argument works for every m. Since local 
unitaries cannot change the Schmidt coefficients, it is obvious that \ip) and \(f)) must share the same 
set of Schmidt coefficients, that is \j = /ij, for every j, and this proves (iii). 

To prove (iv), note that if there is a separable operation carrying j?/;) to 10) and another carrying 
10) to I?/;), the Schmidt ranks of \'\\)) and |0) must be equal by (i), and (|2.5p is an equality, so (iii) 
implies equal Schmidt coefEcients. 

Next let us consider the modifications needed when the Schmidt ranks of 1-0) and |0) might be 
unequal, and are possibly less than the dimensions of T-La or Hb, which need not be the same. As 
noted previously, (j2.9p shows that the Schmidt rank of |0) cannot be greater than that of j^/'). If it is 
less, then the right side of (|2.5p is zero, because at least one of the /ij in the product will be zero, so 
part (ii) of the theorem is automatically satisfied, part (iii) does not apply, and (iv) is trivial. Thus 
we only need to discuss the case in which the Schmidt ranks of {ip) and |0) have the same value r. 
Let Pa and Pb be the projectors on the Ha and Hb supports Sa and Sb of jV') (as defined at the 
beginning of this section), and let Ta and 7b be the corresponding supports of |0). Note that each 
of these subspaces is of dimension r. Since {Pa <8) Pb)|V') = IV^), p.Sp can be rewritten as 



{A'^^B'^)\i;) 



where 



AL 



-^iT7 Pa 



BL 



B„,Pr, 



(2.16) 



(2.17) 



are the operators A^ and i?,„ restricted to the supports of {ip). In fact, AJ„ maps Sa onto Ta, and 
-B,'„ maps Sb onto Tb, as this is the only way in which (|2.16p can be satisfied when |0) and lip) have 
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the same Schmidt rank. Finally, by multiplying (j2.2p by Pa ® Pb on both left and right one arrives 
at the closm'e condition 

E '^'n'^'m ® B'Jb:^ ^Pa® Pb. (2.18) 

m 

Thus if we use the restricted operators A'^ and -B^ we are back to the situation considered previously, 
with Sa and Ta (which are isomorphic) playing the role of T-La^ and Sb and Tb the role of "Hs, and 
hence the previous proof applies. □ 

Some connections between LOCC and the more general category of separable operations are 
indicated in the following corollaries: 

Corollary 2.2. When \^) is majorized by \4>), so there is a deterministic LOCC mapping Itjj) to 
10), there does not exist a separable operation that deterministically maps \(f)) to unless these 
have equal Schmidt coefficients (are equivalent under local unitaries). 

This is nothing but (iv) of Theorem 1 applied when the \'ip) to map is LOCC, and thus 
separable. It is nonetheless worth pointing out because majorization provides a very precise char- 
acterization of what deterministic LOCC operations can accomplish, and the corollary provides a 
connection with more general separable operations. 

Corollary 2.3. If either "Ha orJis is 2-dimensional, then |^/;) can be deterministically transformed 
to 10) if and only if this is possible using LOCC, i.e., is majorized by \<j)) . 

The proof comes from noting that when there are only two nonzero Schmidt coefficients, the 
majorization condition is Xq, and this is equivalent to ()2.5|) . 



2.3 Separable random unitary channel 
2.3.1 Condition for deterministic mapping 

Any quantum operation (trace-preserving completely positive map) can be thought of as a quantum 
channel, and if the Kraus operators are proportional to unitaries, the channel is bistochastic (maps 
/ to /) and is called a random unitary channel or a random external field in Sec. 10.6 of |BZ06j . 
Thus a separable operation in which the Am and Bm are proportional to unitaries and Vm, so 
(|2.ip takes the form 

p' = A(p) = Ep™(^™ ® Vra)p{Um ® Vra)\ (2.19) 
m 

with the pra > summing to 1, can be called a separable random unitary channel. We shall be 
interested in the case in which T-La and "Hb have the same dimension D, and in which the separable 
unitary channel deterministically maps not just one but a collection {|?/'j)}, 1 < j < of pure states 
of full Schmidt rank D to pure states. This means that (|2.8p written in the form 

{U^®V„,)\^j) (2.20) 

must hold for all j as well as for all m. The dot equality ^ means the two sides can differ by at most 
a complex phase. Here such phases cannot simply be incorporated in Um or Vm, because (|2.20p must 
hold for all values of j, even though they are not relevant for the map carrying \'ipj) \ to | ) {4>j\. 

Theorem 2.4. Let {\'ipj)}, 1 < j < N be a collection of states of full Schmidt rank on a tensor 
product T-La ® "Hb of two spaces of equal dimension, and let A be the separable random unitary 
channel defined by (j2.19p . Letijjj and (f)j be the operators dual to \Tpj) and \(f)j) — see (|2.6p and (|2.7p . 
i) If every \tjjj) from the collection is deterministically mapped to a pure state, then 

ULU^i'ji'i = ^M^LU^ (2-21) 



2. ENTANGLEMENT TRANSFORMATIONS USING SEPARABLE OPERATIONS 16 



for every m,n,j, and k. 

a) If (j2.2ip holds for a fixed m and every n^j, and k, it holds for every m,n,j, and k. If in 
addition at least one of the states from the collection {\ipj)} is deterministically mapped to a pure 
state by A, then every state in the collection is mapped to a pure state. 

Hi) Statements (i) and (ii) also hold when (|2.2ip is replaced with 

ViV^i^]i>k = ^li^kVlV^. (2.22) 

Proof. Part (i). By map-state duality p.20|) can be rewritten as 

U^^jV^, = (Pj, (2.23) 

where Vm is the transpose of Vm — see the remarks following (|2.9p . By combining (|2.23p with its 
adjoint with j replaced by k, and using the fact that Vm is unitary, we arrive at 

Since the right side is independent of m, so is the left, which means that 

Un^ji'iul = Umi'.^Ul (2.25) 

Multiply on the left by U^^ and on the right by Un to obtain p.2ip . 

Part (ii). If p.25|) . which is equivalent to (|2.2ip . holds for m = 1 it obviously holds for all values 
of m. Now assume that jV'i) is mapped by A to a pure state |0i), so (|2.23|) holds for all m when 
j — I. Take the adjoint of this equation and multiply by Vm to obtain 

i^\ul = VmCt>\. (2.26) 

Set /c = 1 in (|2.25p . and use p.26p to rewrite it as 

U^i>,V,A = Um^,VmA- (2-27) 

Since by hypothesis \ipi) has Schmidt rank D, the same is true of ipi, and since C/,„ and Vm in (j2.23p 
are unitaries, 4>i and thus also (/)| has rank D and is invertible. Consequently, (|2.27p implies that 

C/„^,K = Um^jVm, (2.28) 

and we can define 0^ to be one of these common values, for example UiipjV. Map-state duality 
transforms this (t>j into \(t>j) which, because of (|2.28p . satisfies (12.201) . 

Part (iii). The roles of Um and Vm arc obviously symmetrical, but our convention for map-state 
duality makes ipj a map from Hb to Ha, which is the reason why its adjoint appears in ()2.22p . □ 

2.3.2 Example 

Let us apply Theorem 12.41 to see what pure states of full Schmidt rank are deterministically mapped 
onto pure states by the following separable random unitary channel on two qubits: 

A(p)=pp+{l-p){X(g)Z)p{X(g}Z). (2.29) 

The Kraus operators are I (E) I and X (E) Z , so Ui ^ I and U2 = X. Thus the condition (|2.2ip for a 
collection of states {IV'j)} to be deterministically mapped to pure states is 

Xi:,4 = i',4x. (2.30) 
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It is easily checked that 

IV^i) = (I+)|0) + |-)|1))/a/2 
is mapped to itself by (|2.29p . If the corresponding 

^-KJ -0 

is inserted in p.30|) with k — 1, one can show that (|2.30|) is satisfied for any 2x2 matrix 

having Cj ~ ±aj and dj ~ T^j, and that in turn these satisfy ()2.30|) for every j and k. 
states of the form 

\tp±) = a\00) + b\Ol) ± a\lO) T 
with a and b complex numbers, are mapped by this channel into pure states. 

2.4 Conclusions 

Our main results are in Theorem 1 2. II if a pure state on a bipartite system Jij^ is dcterministi- 

cally mapped to a pure state by a separable operation {A„ (X)i?„j}, then the product of the Schmidt 
coefficients can only decrease, and if it remains the same, the two sets of Schmidt coefHcients are 
identical to each other, and the Am and Bm operators are proportional to unitaries. (See the de- 
tailed statement of the theorem for situations in which some of the Schmidt coefficients vanish.) 
This product condition is necessary but not sufhcient: i.e., even if it is satisfied there is no guarantee 
that a separable operation exists which can carry out the specified map. Indeed, we think it is likely 
that when both T-La and T-Lb have dimension 3 or more there are situations in which the product 
condition is satisfied but a deterministic map is not possible. The reason is that (|2.5p is consistent 
with \4>) having a larger entanglement than 1-0), and we doubt whether a separable operation can 
increase entanglement. While it is known that LOCC cannot increase the average entanglement 
[ |HHHH09] . Sec. XV D], there seems to be no similar result for general separable operations. This 
is an important open question. 

It is helpful to compare the product condition (j2.5p with Nielsen's majorization condition, which 
says that a deterministic separable operation of the LOCC type can map |^/;) to \<j)) if and only if \<j)) 
majorizes !■(/;), see[TJ Corollarv 12.31 of Theorem 12. II shows that the two arc identical if system A or 
system B is 2-dimcnsional. Under this condition a general separable operation can deterministically 
map I"!/;) to \(f>) only if it is possible with LOCC. This observation gives rise to the conjecture that 
when either ^ or i? is 2-dimcnsional any separable operation is actually of the LOCC form. This 
conjecture is consistent with the fact that the well-known example [BPF+QO] of a separable operation 
that is not LOCC uses the tensor product of two 3-dimensional spaces. But whether separable and 
LOCC coincide even in the simple case of a 2 x 2 system is at present an open question (see note 
added in proof). 

When the dimensions of A and B are both 3 or more the product condition of Theorem 12.11 is 
weaker than the majorization condition: if |0) majorizes then (j2.5p will hold[f|, but the converse 
is in general not true. Thus there might be situations in which a separable operation determinis- 
tically maps 1-0) to 10) even though |0) docs not majorize |V^). If such cases exist, Corollarv 12.21 
of Theorem 12.11 tells us that ji/;) and \4>) must be incomparable under majorization: neither one 
majorizes the other. Finding an instance, or demonstrating its impossibility, would help clarify how 
general separable operations differ from the LOCC subclass. 

^ The general argument that 112. 5t is implied by (though it docs not imply) majorization will be found in f lNie| . 
Sec. 4], or as an exercise on f lHJ99l . p. 199] 
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When a separable operation detcrministically maps \tjj) to \(f>) and the product of the two sets of 
Schmidt coefficients arc the same, part (iii) of Theorem 12.11 tells us that the collections of Schmidt 
coefficients arc in fact identical, and that the Am and i?,„ operators (restricted if necessary to 
the supports of \^)) are proportional to unitaries. Given this proportionality (and that the map 
is deterministic), the identity of the collection of Schmidt coefficients is immediately evident, but 
the converse is not at all obvious. The result just mentioned can be used to simplify part of the 
proof in some interesting work on local copying, specifically the unitarity of local Kraus operators in 
[ |ACP04j . Sec. 3.1]. It might have applications in other cases where one is interested in deterministic 
nonlocal operations. 

Finally, Theorem 12.41 gives conditions under which a separable random unitary operation can 
detcrministically map a whole collection of pure states to pure states. These conditions [see (j2.2ip 
or (|2.22p ] involve both the unitary operators and the states themselves, expressed as operators 
using map-state duality, in an interesting combination. While these results apply only to a very 
special category, they raise the question whether simultaneous deterministic maps of several pure 
states might be of interest for more general separable operations. The nonlocal copying problem, as 
discussed in [ACP041 IKE061 [GKR041 fOHOGj . is one situation where results of this type are relevant, 
and there may be others. 

Note added in proof. Our conjecture on the equivalence of separable operations and LOCC for 
low dimensions has been shown to be false |DFYj . 
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Separable operations on pure 
states 



3.1 Introduction 

A separable operation A on a bipartite quantum system is a transformation of the form 

JV-l 



p' = A(p) = J2 (^fc ® Bk)p{Ak ® BkY, (3.1) 

fc=0 

where p is an initial density operator on the Hilbert space Ha ® "Hb- The Kraus operators Ak ® 
are arbitrary product operators satisfying the closure condition 

N-l 

Y,Al^k®BlBk=lA®lB, (3.2) 

A;=0 

with Ia and Ib the identity operators. The extension to multipartite systems is obvious, but here 
we will only consider the bipartite case. To avoid technical issues the sums in p.ip and p.2p as well 
as the dimensions Da and Db of Ha and Hb are assumed to be finite. 

Local operations with classical communication (LOCC) form a subset of separable operations in 
which the Kraus operators A/j (g) Bk are restricted by the requirement that they be generated in the 
following fashion. Alice carries out an operation J2i ^i^''^ ^i^^ ~ ^A, in the usual way with 

the help of an ancilla, the measurement of which yields the value of i, which is then transmitted to 
Bob. He uses i to choose an operation {B^^'^^}, the result j of which is transmitted back to Alice, 
whose next operation can depend on j as well as i, and so forth. While it is (fairly) easy to see that 
the end result after an arbitrary number of rounds is of the form (|3.ip , it is difficult to characterize 
in simple mathematical or physical terms precisely what it is that distinguishes LOCC from more 
general separable operations. Examples show that separable operations can be more effective than 
LOCC in distinguishing certain sets of orthogonal states |BDF+99 , even in a system as simple as 



two qubits |DFYj . but apart from this little is known about the difference. 

What we demonstrate in Sec. 13.21 of this chapter is that the ensemble {pk, |0fc)} produced by a 
separable operation acting on a pure state \ip)^ see (j3.5l) . satisfies a majorization condition p.7p . 
which is already known to be a necessary and sufficient condition for producing the same ensemble 
from the same by LOCC. Among the consequences discussed in Sec. 13.31 are: a separable opera- 
tion acting on a pure state can be "simulated" by LOCC; a necessary condition for a deterministic 



3. SEPARABLE OPERATIONS ON PURE STATES 



20 



transformation \ip) — > \(f)) given in |GG07| can be replaced by a necessary and sufficient majoriza- 
tion condition; and certain entanglement measures arc nonincreasing under separable operations. 
Section [3.41 summarizes our main result and indicates some open questions. 

3.2 Ensembles produced by separable operations on pure bi- 
partite states 

3.2.1 Majorization conditions 

Let {Ak ^ -BfejfcLi be a separable operation on Ha <E)'Hb, specified by N Kraus operators satisfying 
the closure condition (|3.2p . Let \ifj) be a normalized entangled state on Ha^Hb with Schmidt form 

D-l 

W = E V^l«^)l^j)' (3-3) 

where D = Db, and we assume without loss of generality that Da > Db- Here {!%)} and 

are orthonormal bases chosen so that the Schmidt weights (coefficients) Xj are in increasing order, 

i.e. 

Ao Ai ••• Ai5-i. (3.4) 
The separable operation acting on will produce an ensemble {pk, \'t'k)}k=i' where 

[Ak ^ BkM = ^\(bk) (3.5) 

and 

Pk = {MAUk®BlBk\^P). (3.6) 

In |JP99] it was shown that such an ensemble {pfc, |0fc)}^^^ can be produced from |-0) by a 
suitable LOCC if and only if the majorization inequalities 

J2 PkE^ilcpk)) Enili;)) (3.7) 

hold for ^ n ^ Z? — 1, where 

n 

Enim = X„(Ti-a|V') (^I ) = E (3.8) 

3=0 

and similarly for the |(/)fc). Here Tr^(|'0) is the reduced density operator of on Bob's side, 

and Xn (-) is defined to be the sum of the first n smallest eigenvalues of its argument. Note that we 
are assuming that D = Db < Da, because Db were greater than Da the extra zero eigenvalues 
in Tr^il?/') (?/'| would cause confusion when using Xn- 
Our main result is the following. 

Theorem 3.1. The ensemble {pk, \4>k)}k=i) '''^'^ produced by a bipartite separable operation acting 
on the normalized state if and only if the majorization condition defined by the collection of 
inequalities in (j3.7p is satisfied. 

Proof. To simplify the proof we assume that Da = Db = D. If Da is larger, one always modify 
each Ak by following it with a suitable local unitary which has the result that as long as the Kraus 
operators are acting on a fixed \ip) the action on the A side takes place in a subspace of Ha of 
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dimension D. These local unitarics do not change the Schmidt weights of the \(j)k) or alter the 
closure condition p.2p . For more details about this "decoupling" see jGG07| . 

When the majorization condition p.7p holds the result in |JP99j guarantees the existence of 
an LOCC (hence separable operation) which will produce the ensemble out of lip). The reverse 
inference, that the ensemble {pk,\'t'k)}k=a defined in p.Sp and p.6p satisfies p.7p . follows from 
noting that 

PkEnil^k)) ^ Xn{TrA[Ak Bk\^) (VI 4® 4])' (3-9) 
and applying Theorem 13.21 below with R ^ Ia ® Ib, corresponding to (|3.2[) . so = 1. □ 

3.2.2 A majorization theorem 

Theorem 3.2. Let Ha andT-Ls have the same dimension D, let \^) he some pure state on 'Ha®'Rb, 
and let {A^ ® Skj^S^ be any collection of product operators on T-La ®T~Lb- Then for every ^ n ^ 
D- 1 

N-l 

Y^Xn{TrA[Ak®Bk\i:){^\Al®Bl]) ^ |li?|lx„(Tr^|V)(V|), (3.10) 

fe=0 

where \\R\\ = s^'P\\u)\\=i the largest eigenvalue of the positive operator 

N-l 

i?= (3.11) 

fe=0 

Proof. By map-state duality [GG071 [ZB041 [GWYC06j . using the Schmidt bases of ji/i), we transform 
the state Ak (E) Bkl^j) to a map AkipBk, where 

D-l 

^= E (3-12) 

denotes an operator mapping 'Hs to T-La, and Bk — B"[ is the transpose of Bk- The matrix of ip 
using the Schmidt bases of \ip) is diagonal, with the entries on the diagonal in increasing order. (See 
Sec. II of jGG07] for more details on map-state duality.) Upon writing the partial traces as 

Tr^lV) (VI = TYA[Ak ® BfciV) (VI 4 ® = AfeVB^B^V^^, (3.13) 

the inequalities p.lOp become: 

^ Xn{Ak^BkBl^^ Al) ^ ||i?||x„(VV'^). (3.14) 

For some n between and D — \ write the diagonal matrix V as 

V = V'n+Vri, (3.15) 

where Vn is the same matrix but with A„, A„+i, . . . set equal to zero, while Vn is obtained by setting 
Ao, Ai, . . . A„_i equal to zero. Lemma 13.31 below, tells us that for each k, 

XniAk^BkBlx};^ Al) ^ Tr(AfcVnBfeBiV,\4)- (3-16) 

By map-state duality, 

Tr(AfcVnBfeBl^J,4) = {'^n\AlAk ® BlBfc|V„) (3.17) 
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where \ipn), the counterpart of ipn, is given by (|3.3p with D replaced by n. Inserting p.l7p in p.l6p 
and summing over k, see yields 

Af-l 

^ XniAki^B.Blij^Ai) < (V„|i?|^„) < \\RmM = WRWxnii^H). (3.18) 

This establishes p.l4p . which is equivalent to (|3.10p . □ 

Lemma 3.3. Let A, B , and -0 be D x D matrices, where "0 is diagonal with nonnegative diagonal 
elements in increasing order, and for some ^ n ^ D — 1 let ipn be obtained from ip by setting all 
but the n smallest diagonal elements equal to 0, as in (|3.15p . Then 

XniA^pBB^'^A^) TiiAiPnBB^iPlA''). (3.19) 

Proof. The inequality 

XniAiJjBB'f'ip'iA'i) ^ TxiPnAi^BB^'^ A'^ Pn), (3.20) 

where P„ is a projector (orthogonal projection operator) of rank at least n, follows from the fact that 
for any Hermitian operator T the sum of its n smallest eigenvalues is the minimum of Tr(P„rP„) 
over such P„, see page 24 of [Bha97j . Choose P„ to be the projector onto the orthogonal complement 
of the range of Aijjn, where ipn = V' ~ V'ru in (|3.15p . The rank of Aipn is no larger than the rank 
of Tpm which is smaller than or equal to D — n. Thus the dimension of the range of Aipn cannot 
exceed D — n, so the rank of P„ is at least n. By construction, PnAtpn = 0, so 

PnA4> = PnA{lPn + i>n) = P„^^„. (3.21) 

Thus with this choice of P„ the right side of p.20p is 

Tr(P„A^.„Ppt^,t^tp^)^ (3.22) 

and this implies p.l9p . since P„ ^ / and AtpnBB^ipl^A^ is positive. □ 



3.3 Consequences 

The following are some consequences of Theorem 13.11 

i) An ensemble {p^, |0fe)} can be produced by a separable operation acting on a bipartite en- 
tangled state \ip) if and only if it can be produced by some LOCC acting on the same state 

IV'). 

ii) For a given bipartite \Tp) and separable operation {Ai; (g) Pfcj^J^^ there is another operation 
of the form {Ai (g) Ui}ff^Q^, where the Ui are unitary operators (and the closure condition is 
J2iLo^ ~ ^a), which produces the same ensemble when applied to \ip). Here M could be 
different from N, as two Kraus operators might yield the same \(f>k)- For more details about 
the relation between the {Ak, Bky^^SQ set and the {Ai (g) Ui}^J^^ set see [LPOlj . 

iii) A deterministic transformation ji/i) — |0) by a separable operation is possible if and only if 
En{\4')) ^ En{\ip)) for every n between and D — 1, with En{.) defined in (|3.8p This is often 
written as A^, -< A^, where A^ and A^ are vectors of the corresponding Schmidt weights. (This 
extends Theorem 1 in |GG07| .) 

iv) The maximum probability of success for the transformation |?/') — > \(f>) by a separable operation 
is given by 



^^^1^(1^)^!^))= nrin JiMl, (3.23) 



where En{-) was defined in 
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v) An entanglement measures E defined on pure bipartite states is nonincreasing on average 
under separable operations, whieh is to say 

N-l 

if and only if it is similarly nonincreasing under LOCC. 

vi) Let 

E{p)=m^Y,P^Em), (3.25) 

i 

with the infimum over all ensembles {pi, \'ipi)} yielding the density operator p ~ ^^vA'^i) {^iV 
be the convex roof extension of a pure state entanglement measure E that is monotone on 
pure states in the sense of p.24p . Then E is monotone on mixed states in the sense that 

E{p) ^ J2 PkE{<Jk) (3.26) 

k=0 

for any ensemble {pfc,(Tfc} produced from p by separable operations. 

The result (i) is an immediate consequence of Theorem 13.11 as the same majorization condition 
applies for both separable and LOCC. Then (ii), (iii), and (iv) are immediate consequences of 
known results, in }LP01| . |Nie99| . and }Vid99| . respectively, for LOCC. The result (v) is an obvious 
consequence of (i), whereas (vi) follows from general arguments about convex roof extensions; see 
Sec. XV.C.2 of |HHHH09j . 

3.4 Conclusion 

We have shown that possible ensembles of states produced by applying a separable operation to 
a bipartite entangled pure state can be exactly characterized through a majorization condition, 
the collection of inequalities (|3.7p for different n. These have long been known to be necessary 
and sufhcient conditions for producing such an ensemble using LOCC, so their extension to the 
full class of separable operations is not altogether surprising, even if our proof is not altogether 
straightforward. 

Connecting the full set of separable operations with the more specialized LOCC class immediately 
yields several significant consequences for the former, as indicated in the list in Sec. 13.31 because much 
is already known about the latter. Of particular significance is that various entanglement measures 
are monotone, meaning they cannot increase, under separable operations — something expected on 
intuitive grounds, but now rigorously proved. Since such monotonicity under LOCC has long been 
considered a necessary, or at least a very desirable condition for any "reasonable" entanglement 
measure on mixed states (see Sec. XV. B of |HHHH09] ) . one wonders whether monotonicity under 
separable operations, in principle a stronger condition, might be an equally good or even superior 
desideratum. 

Our results apply only to bipartite states, but separable operations and the LOCC subclass can 
both be defined for multipartite systems. Might it be that in the multipartite case the ensemble 
produced by applying a separable operation to a pure entangled state could also be produced by 
some LOCC applied to the same state? It might be, but proving it would require very different 
methods than used here. There are no simple multipartite analogs of the Schmidt representation 
p.3|) . the majorization condition p.7|) . or map-state duality. 

Even in the bipartite case we still know very little about separable operations which are not 
LOCC, aside from the fact that they exist and can be used to distinguish certain collections of 
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orthogonal states more effectively than LOCC. The results in this chapter contribute only indirectly 
to a better understanding of this matter: looking at what a separable operation does when applied 
to a single entangled state will not help; one must ask what it does to several different states. 
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4 

Local cloning of entangled states 
by separable operations 

4.1 Introduction 

As summarized by the "no-cloning" theorem of |WZ82j , any set of quantum states can be determin- 
isticahy cloned if and only if the states in the set are mutually orthogonal. When the set consists of 
bipartite entangled states, and the cloning is restricted to local operations and classical communica- 
tion (LOCC), the problem becomes much more difficult, and further restrictions have to be imposed. 
The mere orthogonality of the states no longer implies that they can be (locally) cloned. 

The local cloning protocol of a set of bipartite entangled states S = {IV'i)'^^} is schematically 
represented as 

1^.)^^ ® 1^.)^^ ® \Ar\ vz, (4.1) 

where the letters A, a label Alice's systems and B, b label Bob's systems. Both parties are assumed to 
have access to ancillary qudits and may share a classical communication channel, so that in principle 
any LOCC operation can be performed. The state |0) is shared in advance between the parties, and 
it plays the role of a "blank state" on which the copy of \tjji) is to be imprinted. 

The local cloning problem has recently received a great deal of attentio n |GKR04l IACP041 IOH061 
IKE06[ ICKRR07] . and was partially extended to tripartite systems in |CKK+07] . The question 
addressed in all previous work was which sets of states S can be locally cloned (by LOCC) using a 
given blank state 

Note that if one can use LOCC to transform \(f)) into three maximally entangled states of sufficient 
Schmidt rank, then the local cloning of any set of bipartite orthogonal entangled states becomes 
trivially possible, using teleportation: Alice uses one maximally entangled state to teleport her part 
of I'i/'i) to Bob, who then distinguishes it (i.e. learns i), and next communicates the result back to 
Alice. Now both Alice and Bob know which state was fed into the local cloning machine. Finally 
they transform deterministically the two remaining maximally entangled states into [ipi) \ipi) by 
LOCC, which is always possible, according to jNie99| . 

Another possible scenario that uses only two entangled blank states involves using LOCC to 
deterministically distinguish which state \^pi) was fed into the local cloning machine, which can 
always be done if there are only two states in the set S |WSHVOO] . Then, knowing the state, one 
can deterministically transform the two blank states into \ipi) (S) \ipi) (by LOCC). In this case, one 
needs at least two maximally entangled resource states, one for each of the two copies that must 
now be created, since in general the entanglement of the original state will have been destroyed in 
the process of distinguishing the states |Coh07] . 
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One might hope, however, that local cloning can be performed using even less entanglement. As 
first shown in |GKR04] . this hope is sometimes correct. Any two (and not more) two-qubit Bell 
states can be locally cloned using a two-qubit maximally entangled state. 

This result was further extended in jACP04j and jOH06| . which considered local cloning of 
maximally entangled states on higher-dimensional D x D systems using a maximally entangled 
resource of Schmidt rank D. First, necessary and sufficient conditions for the local cloning of two 
maximally entangled states were provided in |ACP04j . which also proved that for _D = 2 (qubits) 
or I? = 3 (qutrits) , any pair of maximally entangled states can be locally cloned with a maximally 
entangled blank state. Whenever D is not prime the authors showed that there always exist pairs of 
maximally entangled states that cannot be locally cloned with a maximally entangled blank state. 
A generalization to more than 2 states but prime D was given in |OH06| . which showed that a set 
of D maximally entangled states can be locally cloned using a maximally entangled resource if and 
only if the states in the set are locally (cyclically) shifted 



where the © symbol denotes addition modulo D. 

Kay and Ericsson [KE06| extended the above results to the LOCC cloning of full Schmidt rank 
partially entangled states using a maximally entangled blank state. They presented an explicit 
protocol for the local cloning of a set oi D x D cyclically shifted partially entangled states 



and asserted that (|4.3p is also a necessary condition for such cloning; that the states to be cloned 
must be of this form. Unfortunately, the proof is not correct 0, and therefore finding necessary 
conditions when the states are partially entangled remains an open problem. 

In this chapter, we consider a set S = {I'fpi}^'^} of full Schmidt rank qudit (of arbitrary dimension) 
partially entangled states. Actually, we will begin by considering sets S in which only one state is 
required to be full Schmidt rank, and then we will see that in fact, all states in S must be full rank. 
Previous work assumed the blank state |0) to be maximally entangled, but in the present chapter 
we do not impose any a priori assumptions on |0) and find that its Schmidt rank must be at least 
that of the states in S. Furthermore, we do not restrict to LOCC cloning, but allow for the more 
general class of separable operations — all the necessary conditions we find for separable operations 
will also be necessary for LOCC since the latter is a (proper) subset of the former [BDF+Qij] . 

The remainder of the chapter is organized as follows. In the next section we give a preliminary 
discussion and define some terms that will be used. Then, in Sec. 14.31 we turn to the characterization 
of clonable sets of states, where we show that \<j>) and all states in S must be full Schmidt rank, 

^The matter was discussed with Kay [Kay06| . The fact that the argument is not correct can be observed 
after a careful reading of the paragraph following Eq. (3) in IKE06I . The authors claim that the local cloning 
of partially entangled states is equivalent to the cloning of maximally entangled states, but this statement is 
incorrect, because the authors implicitly modified the Kraus operators that defined the local cloning, i.e. changed 
j4j. to = AfeMo, where Mq (defined in Eq. (3) of IKEOGI ) is the operator that transforms the maximally 

entangled state (1/v^) Z^^J^^ |r->'*|r->^ to the partially entangled state I'/'o) = Zl^Jo^ v^k)'^k>^- The new 
Kraus operators do not satisfy the closure condition anymore (necessary for a deterministic transformation), since 
J2k ^'k^^k ® Bk^'Bk = Y.k ^li^k^ ^k)Mo (gi Bk^Bk = M^Mo ® I I ® I, because Mo is not a unitary operator 
(unless \tl)o) is maximally entangled, case excluded). 

Another way of seeing that the argument is not correct is to observe that, if the Bk operator performs the 
cloning of a maximally entangled state using a maximally entangled blank, as it is claimed, then Bk must be 
proportional to a unitary operator, see Theorem l(iii) of |GG07) and Sec. 3.1 of |ACP04| . It then follows that the 
closure condition for the Kraus operators is not satisfied, with Ak as defined in Eq. (3) of IKE06| . 




(4.2) 



D-l 




(4.3) 



r=0 
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provide additional necessary conditions on S, and then prove the group structure of these sets. From 
this group structure, it is then shown that the number of states in S must divide D exactly, and this 
is followed by a proof of a necessary ( "group-shifted" ) condition on the local cloning of a set of 13 x D 
maximally entangled states. Then, in Sec. 14.41 we further consider group-shifted sets, now allowed to 
be not maximally entangled, showing that a maximally entangled blank state is sufficient by giving 
an LOCC protocol that clones these states. This demonstrates that the necessary condition found 
in the previous section for cloning maximally entangled states is also sufficient for LOCC cloning. In 
Sec. 14.51 we provide necessary conditions on the minimum entanglement in the blank. In addition, 
we obtain necessary and sufficient conditions for local cloning of any set when D = 2 (entangled 
qubits), and for any group-shifted set for Z3 = 3 (entangled qutrits); in both these cases we find 
that the blank state must be maximally entangled, even when the states to be cloned are not. For 
higher dimensions with these group-shifted sets, we also show that the blank must have strictly more 
entanglement than the states to be cloned. Finally, Sec. 14.61 provides concluding remarks as well as 
some open questions. Longer proofs are presented in the Appendices. 

4.2 Preliminary remarks and definitions 

A separable operation A on a bipartite quantum system "Ha ® is a transformation that can be 
written as 

M-l 

p' = A{p) = ^ (A™ ® B,n)p{A^ <E) B^)^ (4.4) 

where p is an initial density operator on the Hilbert space T-La <E) Hb- The Kraus operators are 
arbitrary product operators satisfying the closure condition 

Ai,A^<»Bl^Bra^lA(8)lB, (4.5) 

with Ia and Ib the identity operators. The extension to multipartite systems is obvious, but here 
we will only consider the bipartite case. To avoid technical issues the sums in (|4.4p and (|4.5p . as 
well as the dimensions of Ha and H b , are assumed to be finite. 

The local cloning protocol is described as follows. Suppose Alice and Bob are two spatially 
separated parties, each holding a pair of quantum systems of dimension D, with Alice's systems 
described by a Hilbert space Ha <8) Ha and Bob's by Hb <8) Hb- Let S = {|V'i)'^^}£o"^ be a set 
of orthogonal bipartite entangled states on Ha ®Hb- Let e Ha ® Hb be another bipartite 

entangled state that plays the role of a resource, which we call the blank state, and is shared in 
advance between Alice and Bob. Their goal is to implement deterministically (i.e. with probability 
one) the transformation 

m^"" ® W 1^,;)-^'' ® \i^^T\y^ = Q...N-l (4.6) 

by a bipartite separable operation. Alice and Bob know exactly the states that belong to the set 
S and also know the blank state ^ but they do not know which state will be fed to the local 
cloning machine described by (|4.6p — the machine has to work equally well for all states in 5! Note 
that local cloning is defined up to local unitarics, i.e., a set S ~ can be locally cloned 

if and only if the set S' ~ {U^ ® \'4'i)^^}f=o^ can be locally cloned, where and are local 
unitaries. This is true because local unitaries can always be implemented deterministically at the 
beginning or at the end of the cloning operation. 

The Schmidt coefficients of \il>i)^^ are labelled by Af*"* and by convention are sorted in decreasing 
order, with \ f ^ A^''^ > ■ ■ ■ ^ A^^i and Y.r=^ ^r*' = 1, for aU i = . . . TV - 1, and the Schmidt 
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coefficients of are labelled by 7^, with 70 ^ 71 ■ • • 5^ 7_d-i and ^f^^ 7^ = 1. To remind the 
reader that the components of a vector A are ordered in decreasing order we use the notation A^. 

The Schmidt rank of a bipartite state is the number of its non-zero Schmidt coefficients. We say 
that a state of a _D x _D dimensional system has full Schmidt rank if its Schmidt rank is equal to D. 

We use the concept of majorization, which is a partial ordering on £)-dimensional real vectors. 
More precisely, if af = (xq, . . . ,xd-i) and y = (yo, • ■ • tVd-i) are two real D-dimensional vectors, 
we say that x is majorized by y and write a; ^ y if and only if X)j=o ^ '^j=n Vj holds for all 
fc = 0, . . . , I? — 1, with equality when k = D — 1. 

For two D X D bipartite pure states |x) and I?;), we use the shorthand notation |x) ^ I??) to 
denote the fact that the vector of Schmidt coefficients of \x) is majorized by the vector of Schmidt 
coefficients of |?7). See [Nie99j or Chap. 12.5 of [NCOOj for more details about majorization. 

The entanglement of a -D x _D bipartite pure state \x) can be quantified by various entanglement 
measures 0, the ones used extensively in this chapter being the entropy of entanglement 

D-l 

E{\x}) = -J2^r\0gaXr (4.7) 

and the G- concurrence |Gou05j 

/D-l 

CG{\x))=D\r[Xrj , (4.8) 

where Xr denotes the r-th Schmidt coefficient of The base D in the logarithm in (|4.7p as well as 
the prefactor D in (|4.8p appear for normalization purposes, so that the entropy of entanglement as 
well as the G-concurrence of a maximally entangled state are both 1, regardless of the dimension. 

4.3 Characterizing sets of clonable states 
4.3.1 Preliminary analysis 

Mathematically, the local cloning problem can be formulated in terms of a separable transformation 
on a set of pure input states S = {|V'i)'^^}£o^' using a blank state !</>)"''. 

If a set of states S can be locally cloned using the blank state |(/>)°^, then there must exist a 
bipartite separable operation A for which 

A(|V.)('^.|^'' ® I0)(0r') = I^.X^,:!^"" ® \^^){^^\^'', yi = O...N-l (4.9) 

(note here that an overall phase factor in the definition of the individual states is of no significance). 
Since A is separable, it can be represented by a set of product Kraus operators, 

M-l 

J2 {Am ® B,nm^){iJ^\^'' ® |0)(0r^)(A™ ® B^)^ 
m=0 

= \^^){H'''' yi^O...N-l, (4.10) 

where operators Am act on Ha ® Ha, and Bm ohHb ® Ht- The above equation is equivalent to 

Am ® Bmrn^"" ® = vp^e'^- ® lij^r"), 

Vi = 0...7V-l, Vm = O...M-l (4.11) 

■^Oftcn called entanglement monotones, i.e., non-increasing under local operations and classical communication 
(LOCC). 
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where e"''™' is a complex phase that may depend on m and i, and pmi are probabihties for which 

= Vi = 0...iV-l. (4.12) 

By map-state duahty in the computational basis El fZBOil IGWYCnei iggoti Iggos] one can 
rewrite (|4.1ip as 

^m(V'*'»0)S^ = \^™e'^'-V«®V'z, Vi,m, (4.13) 

where ipi and (j) are now operators obtained from the corresponding kets by turning a ket into a bra, 
and is the transpose of Bm- 

The superscripts in (|4.13p that label the Hilbcrt spaces have been dropped for clarity, since 
now one can regard everything as abstract linear operators, or matrices in the computational basis. 
Although map-state duality is basis-dependent, our results will not depend on the choice of a specific 
basis. 

We now state our first result characterizing sets of states S that can be locally cloned. 

Theorem 4.1 (Rank of states in S). Let S = {|V'i)"^'^}i!=o^ bipartite orthogonal states 

on "Ha ®'Hb with one state, say IV'o); having full Schmidt rank. If the local cloning of S is possible 
by a separable operation using a blank state \4>) , then and all states in S must be full rank. 

Proof. This result follows directly from (|4.13p . If |V'o) has full Schmidt rank, then -0o is a full rank 
operator. Then, since the rank of a tensor product is the product of ranks, it must be that (j) is also 
full rank, as are A„i and i?,„ for each m (a product of operators cannot have rank exceeding that 
of any of the individual operators in the product). Since the rank of a product of two operators is 
equal to that of the first whenever the second is full rank, (|4.13p with i ^i) directly implies that -04 
has full rank for every i, and we are done. □ 

In this chapter, we are considering sets S in which at least one state is full rank. Therefore by 
this theorem, we may instead restrict to sets in which every state is full rank, and we will do so 
throughout the remainder of the chapter. 

As just argued in the proof of the previous theorem, all operators in (|4.13p are full rank, hence 
invertible. Now take the inverse of (|4.13p . replace i by j, and right multiply (|4.13|) by it to obtain 

A™(V.V''' ® I)A-;n' - ,/^e'(^-'-'^-)(V.V7' ® i^ri^j')- (4.14) 
■' V P™i 

Define 

^ V Pmj ^ 

so that (|4.14p can be written more compactly as 

A„,{T^f ® I)A-^ = T^f ® T\f. (4.16) 

Since for every i, is full rank, we see that det(-0i) ^ 0, so det(T^^™'') is also non- vanishing. Thus, 
taking the determinant on both sides of (|4.16p yields 

det(4"')'' = 1, (4.17) 

'^As an example of map-state duality, a bipartite state |x)^^ S ® Wb, Ix)'^^ = E lO'^li)^! transformed 
into a map x : V.B — ^ ^A, X = E 01^. Note that the rank of the operator x the Schmidt rank of |x)^^i 

and the squares of the singular values of x (or, equivalently, the eigenvalues of XX^) s-i'c the Schmidt coefficients of 
I^^AS Pqj, mQj-e details about map-state duality see Sec. II of |GG07) . 
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where we have used the fact that det(A ® i?) = det(A)^^ det(i?)^, for A and B being N x N and 
M X M matrices, respectively. RecaUing the definition of t/^™'' in (|4.15p . this condition becomes 



\Prnj) 



Pmj Pmi 



det(V'.) 



det(?/'j 



det(V',) 



(4.18) 



(4.19) 



det(?/'j) 

|det(^,)l = |det(V'j)| , (4.20) 



Summing (|4.19p over m yields 
implying 

Pmi = Pmj, (4-21) 

hence these determinants and probabilities are independent of the input state. As a consequence, 
we may write T^^^ in the simpler form, 

T^p) = e'('^--'^'"^V.V'7'- (4.22) 

Observation: The fact that Pmi = Pm, independent of i, implies that the cloning apparatus 
provides no information whatsoever about which state was input to that apparatus, nor can any such 
information "leak" to an external environment that might be used to implement the local cloning 
separable operation. This is not without interest, since it rules out the possibility of local cloning by 
locally distinguishing while preserving entanglement [Coh07| . This result turns out to be valid in the 
much more general setting of one-to-one transformation of full Schmidt rank pure state ensembles 
by separable operations, but a discussion of these broader implications will be presented in a future 
publication. 

We can now provide additional conditions that must hold in order for 5 to be a set of states 
that can be locally cloned by separable operations. These are stated in the following theorem, which 
holds under completely general conditions, applicable for any N and D. 

Theorem 4.2 (Necessary conditions). Let S ~ il'ipi)^^} be a set of full Schmidt rank bipartite 
orthogonal entangled states on 'Ha®'Hb ■ If the local cloning of S using a blank state S Ua^Ub 
is possible by a separable operation, then the following must hold: 

i) All states in S are equally entangled with respect to the G-concurrence measure, 

CGm^'')^CGm^^), Vz,j. (4.23) 

ii) Any two states in S must either share the same set of Schmidt coefficients or he incomparable 
under majorization. 



Hi) 



Spec(i;f ) ® /) = Spec(7;^) ® 4"^), j, (4.24) 



where Spcc{-) denotes the spectrum of its argument and T^^™'' is defined as in ()4.22|) . 

Proof. Proof of i) This follows at once from (|4.20p . the definition (|4.8p of G-concurrence, and the 
fact that for any state \x) the product of its Schmidt coefficients is equal to | det(x)p. 

Proof of ii) The proof follows from Theorem 1 (ii,iii) of |GG07| which states that any two bipartite 
states Ix) and I77) that are comparable under majorization (i.e. |x) ^ Iv) or Iv) ^ Ix)) and have 
equal G-concurrence must share the same set of Schmidt coefficients. 

Proof of iii) The proof follows at once from (|4.16p . □ 
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4.3.2 Characterization of clonable sets in terms of finite groups 

We next show that to any set S of states that can all be cloned by the same apparatus, there can 
be associated a finite group, and the set is essentially generated by this group. 

Theorem 4.3 (Group structure of S). Let S = a set of full Schmidt rank bipartite 

orthogonal entangled states on 'Ha®T~Lb- If the local cloning of S is possible by a separable operation, 
then the set S can be extended to a larger set such that {T'/^™''} of (|4.22p for fixed j,m constitutes 
an ordinary representation of a finite group, G. Since the states in S are related as e"'''"' iTAi) ~ 
Q^Vmj ^j^ip) (g) Ig'^^ipj'j^ then the larger set, with N = \G\ members, is generated by the action of the 
group G on any individual state in the set. 

Proof. The starting point of the proof is to multiply (|4.16p on the left of (|4.13|) (with index k) to 
obtain 

ArniT^fi'k ® <l^)Bl = VP^^c^'^-'T^f^^Pk ® rlp^^bk- (4.25) 
Using (|4.22p this becomes 

AU^.^j'^k <E) <l>)Bl - ^e'('^'"'-'^-+'^'"'')V^^7Vfe ^^^7 Vfe, (4.26) 

which by map-state duality implies that the state It/jitl^J^tpk) is cloned by the same apparatus as all 
the states in the original set S. Therefore lipiipj^tpk) — which, by considering the version of (|4.26p 
that corresponds to states (as in (|4.1ip ). taking the squared norm of both sides and summing over 
TO, is seen to be normalized — must cither (i) be orthogonal to the entire set S, or (ii) it is equal 
to one of those original states up to an overall phase factor. If this state is orthogonal to S, then S 
can be extended by including this state as one of its members. So assume S has been extended to 
its maximal size (since we arc working in finite dimensions, this size will be finite), and then we can 
conclude that for every i,j, k, 

^^i^~'^k = e'('^'"'-'^'-+'^--'^'"'=)V'/, (4.27) 

for some /, where the phase in the above expression has been determined by comparing (|4.26l) to 
(|4.13p . Next multiply this latter expression on the right by e~"^""'(/'n^ to obtain 

Tt^Ttf=Tt\ (4.28) 

Hence the collection of Pj^^^ is closed under matrix multiplication, which is associative. In addition, 
Tj^™^ = / for every i and T^p^TjP^ = I for every «, j, so we see that the identity element and inverses 

are present, which concludes the proof that the set {T^^^} with fixed to form a representation of 
a group, G. Now, the number of index pairs (i, j) is N'^, where N is the number of states in S. 
However, we will now show that in fact the order \G\ of this group is equal to N and not A^^. 
Setting n = j in (|4.28p , we have 

^Hj ^kj -^ij 

so the product is closed even when the second index is constrained to be the same. If we set I = j, we 
see that with tJJ"^ = /, then for each i there exists k such that T^™^ = (Ti™^)^^- Hence, for every 

fixed j the set Tj ~ {Tij™"*} also is a representation of G. Similarly, one can show the same holds if 
instead it is the first index that is held fixed. Note now that by multiplying (|4.28p on the right by 
(T^™-*)"^, and given that (|4.28p holds for any i,j,k,n, we sec that for every i,j, Tij is a member 
of the group formed by the T^n for fixed n. That is, the group of the Tkn for fixed n contains all 
elements Tij. 



T^fTif=Tlf, (4.29) 
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Could two or more of the T^^^ be equal, for fixed j? We will now show this is not the case by 
demonstrating the linearly independence of the set Tj . Indeed, 

N-l N-l 

A;=0 k=0 
N-l 

^ Cfce''^'"'=Vfc- (4.30) 

fc=0 

However, the ipk are mutually orthogonal, TT{ipj^ipj) = Sjk, so this can only be satisfied if all the c^. 
vanish, implying that Tj is linearly independent, and hence, that \G\ = N: the (maximal) number 
of states in S is equal to the order of G. □ 

For the remainder of the chapter, we will use labels f,g,h instead of i,j,k, where the former 
represent elements of the group G; the group multiplication is denoted as fg, with e the identity 
element. For example, instead of ipo wc will now write ipc and in place of T^^™"* we will simply write 

rpim.) 

■ 

We may now utilize the powerful tools of group theory to study sets S of clonable states, obtaining 
a very strong constraint on how many states any given apparatus can possibly clone. Any group 
G is characterized by its irreducible representations, which we denote as r^")(/), f G, and any 
representation of G may be decomposed into a direct sum of irreducible representations with a 
given irreducible representation r'"-'(/) appearing some number Ua times in that sum. In general, 
a given representation may have Ua = for some a, but since here our representation is linearly 
independent, we know that every irreducible representation must appear at least once |YGC10j . 

We can use character theory |IIam89j to calculate n^,. Defining characters as x(^/™^) = Tr(Tj'"'') 
and x'"H/) = Tr(r(")(/)), we have that 

"a = ^Ex<"H/)*^(^r^)- (4.31) 

However, by taking the trace of (|4.16p and recalling that the trace of a tensor product is equal to 
the product of the traces, we see that x(2^/'"^) equal to either or D. Since the identity element 
of the group, e, is always in a conjugacy class by itself (e = sfg^^ if and only if / = e), then we 
may conclude that x(2^/™'') vanishes except when / = e, in which case xiTe"^^) — D. Thus, (j4.3ip 
reduces to 

Ddr, 

(4.32) 



\G\ 



where da = X^^H^) is the dimension of the a'^^ irreducible representation. Since in every finite 
group there is always the trivial irreducible representation of all ones, F*-*' (/) = 1 V/ S G, where 
this irreducible representation has dimension dt = 1, we have immediately that rit = D/\G\ is an 
integer, implying that N =\G\ divides D. Thus, 

Theorem 4.4 (Number of clonable states). // an apparatus can locally clone more than one state 
on a D X D system, where at least one ( and therefore all, see Theorem |^. J[ ) of these states has full 
Schmidt rank, then that apparatus can in fact clone a number of states that divides D exactly. In 
particular if D is prime, then any such apparatus can clone exactly D states, no more and no less. 

Now we see from ()4.32p that Ua is an integer multiple of da. If |G| = D so that Ua = da, we 
have what is known as the regular representation of G. Otherwise, our representation is a direct 
sum of an integer number nt — D/\G\ of copies of the regular representation. As is well known. 
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there is always a choice of basis in which the matrices in a unitary regular representation appear 
as permutation matrices L{f), with each row (column) having only a single non-zero entry equal to 
one. In this basis, denoted as {|g)}gGG, we have that L{f)\g) — \fg). The representation £(/) is 
called the left regular representation. One can as well use the right regular representation R{f) with 
R{f)\g) ~ |.g/~^), but without loss of generality in the rest of the chapter we restrict only to L{f), 
since for finite groups the right and left regular representations are equivalent |Mic95| . 

In our case the representation will generally not be unitary, so when \G\ ~ D wc will have that 

r)™^ = SL{f)S-\ (4.33) 

for some invertible matrix S. 

In the remainder of the chapter we restrict consideration to \G\ ~ D (or, equivalently, to nt — 1), 
and note that all results obtained in the remainder of the chapter are valid (with small modifications) 
also when \G\ < D. However, the notation becomes a bit cumbersome, so we defer detailed discussion 
about the |G| < D case to Appendix 14. Bl 



4.3.3 Form of the clonable states when all are maximally entangled 

It was shown in |ACP04j that when at least one of the states in S is maximally entangled, then all 
states in S must also be maximally entangled. In this section, we consider such sets, in which case 
the Tj™^ must all be unitary. This follows directly from the fact that when -06 is proportional to 

the identity then ipf is proportional to Tj™\ and also that \^pf) is maximally entangled if and only 
if "0/ is proportional to a unitary. 

We have seen that when N = D, then t]™^ = SL{f)S'^ for some invertible S, and L{f) is the 
permutation form of the regular representation of group G. However, we have 

Lemma 4.5 (Unitary equivalence). For any two unitary representations Tf and L{f) of a finite 
group G, which are equivalent in the sense that Tf ~ SL[f)S~^ for some invertible matrix S, then 
these two representations are also equivalent by a unitary similarity transformation, Tf = WL{f)W^ , 
with W unitary. 

A proof of this lemma is given in Chap. 3.3 of |Ma07j . and we provide an alternative proof in 
Appendix I4.A.11 

What this lemma tells us is that ipf is proportional to WL{f)iljeW^ (since by local unitaries, ipe 
can be made proportional to the identity, we will assume here that this is the case, and then ijje 
commutes with W^), or 



\i.f)=CfiWLif)^W*)Y.\9)^\9)'' 

geG 

i=(W^®Ty*)^|/g)^|5)^, (4.34) 



96G 



where W* is the complex conjugate of W, the states {\g)}g^a are some orthonormal basis, {g\h) = 
Sg^h, and we have omitted an unimportant overall phase (from c/, of magnitude D^^^^) in the last 
line. Note that up to unimportant local unitaries and relabeling of group elements, the set of states 
(|4.34p can be written either as 

m = 4=Y.\f9)%)'' (4.35) 



^geG 



\^f) = 4^Y.\9)^\-f9)''- (4.36) 



DgeG 
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The states above are of a form that we will refer to as "group-shifted" . 

In Section we provide an explicit LOCC protocol that accomplishes cloning of such shifted 
sets of states. Thus, we have 

Theorem 4.6 (Maximally entangled states). A set of maximally entangled states on a D x D system 
can be cloned by LOCC if and only if there exists a choice of Schmidt bases shared by those states 
such that they have a group-shifted form, as in (|4.35p or (j4.36p . 

This extends the result of |OH06j . which applied only for prime D. 

Additionally, we remark that in our protocol presented in Sec. 14.41 there is no need for classical 
communication (the measurement Mr and the additional corrections Q,. appearing in that protocol 
can be omitted when the states to be cloned are maximally entangled) . This result was first proven 
in |ACP04| ■ where it was shown that the Kraus operators implementing the cloning of maximally 
entangled states have to be proportional to unitary operators. A completely different proof of this 
fact was later provided in |GG07j . in which it was shown that a separable operation that maps a 
pure state to another pure state, both sharing the same set of Schmidt coefficients, must have its 
Kraus operators proportional to unitaries; in our case jV'/) ^ \4>) and {ipf) (E) {ipf) do share the same 
set of Schmidt coefficients, since they are maximally entangled. We here have another simple proof 
of this result, since we have proved in Theorem 14.61 that a set of maximally entangled states must 
be group-shifted in order that they can be cloned, and since our protocol in Sec. 14.41 clones any set 
that is group-shifted without using communication. 



4.3.4 Form of the clonable states when D = 2 (qubits) 

Here, we restrict our attention to local cloning of qubit entangled states, D ~ 2. As Z) is prime, 
we know from Theorem 14. 31 that exactly two states can be cloned, S = {\^pe)^^ , li^g)^^}- Both are 
assumed to be entangled (non-product), but not maximally entangled. 

Since there is only one independent Schmidt coefficient for a two-qubit state, any two such states 
are comparable under majorization, and then from part ii) of Theorem 14.21 it follows at once that 
these states have to share the same set of Schmidt coefficients. This is already a surprising result, 
implicitly assumed (but not proved) in recent work on local cloning of qubit states [CKRR07] . We 
can actually prove a stronger condition: not only do the states have to share the same set of Schmidt 
coefficients, but they must also share the same Schmidt basis and be of a shifted form, as summarized 
by the following theorem. 

Theorem 4.7 (Entangled qubits). Let S — {jV'e)'^^, IV'g)'^^} ^ ■sei of 2 orthogonal two-qubit 
entangled states and let A be the largest Schmidt coefficient oflipe)"^^, assumed to satisfy 1/2 < A < 1. 
// the local cloning of S using a two-qubit entangled blank state | is possible by a separable 
operation, then, up to local unitaries (that is, the same local unitaries acting on both states), the 
states must either be of the form 

^ yxio)^|o)^ + \/r^|i)'^|i)^ 

\^g)AB ^ VA|o)^|1)S ^ \/T^|l)'4|0)^ (4.37) 



iV^e)^^ = VA|0)^|0)^ + a/T^|1)'4|1)^ 

\^g)^^ = VA|1)-^|0)^ + VT^\0)^\1)^. (4.38) 

Note that a relative phase e"' may be introduced into jiAg), without altering |V'e), by Alice and Bob 
doing local unitaries on systems A and B, U^-^ ~ |0)(0| -I- e*"'/^|l)(l| (one of them chooses the 
upper sign, the other does the lower, which accomplishes the task up to an unimportant overall 
phase). Therefore, the theorem allows cloning of states with these phases. 
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Proof. First note that without loss of gcncrahty one can always assume that the first state \ipe)^'^ 
is already in Schmidt form, 

= \/A|0)^|0)^ + x/T^|l)^|l)^, (4.39) 

since this can be done by a local unitary map U^'^ ® V^''. Therefore, the operators -0e and ipg 
obtained by map-state duality can be assumed to have the form 

*• ^ ( :: ) • 

where A is the largest Schmidt coefficient of \i^e)^^ and are complex numbers with ^ = 1, 
which is equivalent to the requirement that lijjg) be normalized. 
Orthogonality between these two states implies that 

= VAaoo + Aaii. (4.42) 



Since the only group of order 2 is cyclic with elements e,g and g = e, we have from Theorem 14.31 
that (rg'™^)2 = SL{gfS-^ = I. Thus, we require 

{^9^e'f - ( ) , (4.43) 

where the factor of e"' arises from the phases that appear in the definition of Tg™'\ see (j4.22p . Thus, 
implies 

and either (i) aooVl — A = —an^/X] or (ii) ooi = = aio- The condition that ipg be normalized in 
the latter case (ii), along with (|4.42p and ()4.44p . can only be satisfied if A = 1/2, a case we are not 
considering here. The former case (i) along with ()4.42p implies that ooo = = an (again, assuming 
A ^ 1/2). This concludes the proof, since it implies that has to have either the form ()4.37|) 

or the form (j4.38p . up to an unimportant global phase. □ 

Now one can immediately see that one of the families of states considered in |CKRR07] , of the 
form iV^e) = \/A|0)'4|0)^ +'\/r^|l)^|l)^ and = \/r^|0)^|0)^ - \/A|l)^|l)^ cannot be 
locally cloned with a blank state of Schmidt rank 2, unless they are maximally entangled, case 
already studied in |ACP04) . 



4.4 Local cloning of group-shifted states: explicit protocol 
using a maximally entangled blank state 

Consider now a set of group-shifted partially entangled states S = {|-0/)'^-^}/gG on 'Ha^'Hb, where 
the dimension of both Hilbert spaces Ha and "Hs is equal to D, 

IV'/)^^ = E Va;|5)^I/5)^, (4.45) 
see 

and we remind the reader that throughout this section we restrict to the \G\ = D case (see Ap- 
pendix |4B] for the \G\ < D case). 

In the following we present a protocol that locally clones S using a maximally entangled blank 
state of Schmidt rank D. Our protocol, which works for any group G, is a direct generalization of 
the one presented for the special case of a cyclic group in |KE06| . 
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A 



Pn 



Mr 



Qr 



IV-/) ® \4>) 



B 



Pn 



I :__ 



Qr 



B 



B 



Figure 4.1: Circuit diagram for the local cloning of group-shifted states with a maximally entangled 
blank state. There is no need to perform the measurement Mr and the corrections Qr whenever the 
states to be cloned are maximally entangled. 



Theorem 4.8 (Group shifted states). Let S ~ f)^^} feG be a set of group- shifted full Schmidt 
rank bipartite orthogonal entangled states on Ha 'SS'Hb'is defined by (j4.45p . The local cloning of S 
is always possible using a maximally entangled blank state 10)°^ of Schmidt rank D. 

Proof. Without loss of generality the maximally entangled blank state can be written as 

W^^Y.\^y\h)'- (4.46) 
^ ^ heG 

The local cloning protocol is summarized below and the quantum circuit is displayed in Fig. 14.11 
1. Starting with \ipf)^^ ® 10)°'', both Alice and Bob apply the "controlled-group" unitary 

J2\9){9\^Pa. withP, = J2\9h){h\, (4.47) 

geG heG 

where the permutation Pg acts on system a (6) and is controlled by system A {B), to obtain 



geG 



E ^j\9)^\f9)''4= E \9m.f9h)'' 

3 heG 
Ev^l5>-'l/3)''^El'^)"l/'^)' 



(4.48) 



geG ^ heG 

2. Next Alice performs a generalized measurement on system a with Kraus operators 

Mr = E VX^I^> {h\ , E ^'^r^'Mr = /, (4.49) 

heG reG 

and communicates the result r to Bob. Conditioned on the result r, the output state is 



geG heG 



(4.50) 
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3. Both Alice and Bob apply the unitary correction 
on systems a and 6, respectively, to obtain 

gGG /iSG 

= l^;)^^0l^/r^ (4.52) 

which is the desired output. 

□ 

Note that from symmetry considerations states of the form J2geG \/^\f9)^\9)^ (with the term 
fg appearing now on Alice's side instead of Bob's side) can also be locally-cloned, by interchanging 
the roles of Alice and Bob in the protocol, e.g. performing the measurement Air on system b instead 
of a, then sending the result back to a. Therefore in the following, when discussing group-shifted 
states, we will restrict to the states of the form ()4.45p . 



4.5 Local cloning of group-shifted states: minimum entan- 
glement of the blank 

Here again, we restrict for simplicity to the \G\ ~ D case, and discuss the extension of the results 
for |G| < £> in Appendix |H1 

4.5.1 Necessary conditions for arbitrary D 

We now turn our attention to the task of characterizing the blank state, which essentially amounts 
to determining the amount of entanglement it must have in order for the local cloning to be possible. 
Wc first give a very general lower bound as. 

Theorem 4.9 (Minimum entanglement of the blank). Let S = {|i/'/)"^^}/6G be a set of full Schmidt 
rank bipartite orthogonal entangled states on Ha ®'Hb- If the local cloning of S using a blank state 
10)°'' € T-La ®T~Lb is possible by a separable operation, then it must be that 

Ent{\(j))''^) ^ max£;ni(|V'/)^^), (4.53) 
where Ent(-) denotes any pure-state entanglement measure. 

Proof. Wc recently proved in jGGOSj that any pure state entanglement monotone is non-increasing 
on average under the general class of separable operations. The theorem follows directly, since 
otherwise the local cloning machine increases entanglement across the Aa/Bb cut. □ 

Providing a more detailed lower bound appears to be difficult in general, but turns out to be 
possible in the special case of group-shifted states. 

Consider again the set of D group-shifted entangled states (|4.45p . and allow for arbitrary phases, 

m^"" = Y.^9<^'''-\9)^\f9f- (4.54) 

geG 
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Without loss of generality, the blank state can be written as 



^""-"^ /Wmh)", (4.55) 



heG 



where are its Schmidt coefficients, J2heG^>^ ~ ^• 

All states in S have the same Schmidt coefficients, and hence the same entanglement. As shown 
above, the local cloning of the above set of states is possible using a maximally entangled blank state 
when all phases e"*^" are chosen to be f , but it is not yet known if one can accomplish this task 
using less entanglement. One might hope that the local cloning of S is possible using a blank state 
having the same entanglement as each of the states in 5, which could be regarded as an "optimal" 
local cloning. However we prove below that such an optimal local cloning is impossible with these 
states. Indeed we find a sizeable gap between the entanglement needed in the blank state and the 
entanglement of the states of S. For D = 2 and D = 3, we prove that a maximally entangled blank 
state is always necessary. 

In the rest of this section we will use the rearrangement inequality (see Chap. X of |HLP99j ). 
which states that 

XnUi H h xiyn ^ 2:^(1)2/1 H h a;„(„)J/„ ^ xij/i H \- a;„y„ (4.56) 

for every choice of real numbers xi ^ • • • ^ a:„ and yi ^ • • • ^ ?/„ and every permutation 

Xct(1)i ■ ■ ■ iXa(n) of , . . . , Xn ■ 

The following Lemma is the most important technical result of this section (note that in the 
statement of this result, we will use 'g for inverses of elements in the group G, which will make 
the notation somewhat more readable). 

Lemma 4.10 (Majorization conditions). Let S ~ {\i^s)^^} fee be a set of D group-shifted full 
Schmidt rank bipartite orthogonal entangled states on Ha ®'Hb o,s defined by (j4.54p and considered 
to he not maximally entangled. If the local cloning of S using a blank state 10)"'' is possible by a 
separable operation, then 

i) The majorization condition, 

a<l3, (4.57) 
must hold. Here, a and j3 are vectors with components indexed by elements g,h £ G, 

ag.h = 7/1 51 Pa^i^ H l^J^f9^fh^ (4-58) 

/eG /GG 

and {A*/}/eG an arbitrary set of non-negative real coefficients that satisfy X]//^/ ~ 1- 

it) The smallest Schmidt coefficient 7inin of the blank state has to satisfy 

Ym\ig^heGllfeGl^J^f9>'fh 

7mi„ ^ max : =4 . (4.59) 

{m/} mnigeG 2^ f^G l^J^fg 

Hi) In particular, a good choice of {fif} is given by 

fif = ^, with ry-i = l/Ag, (4.60) 
seG 

for which (j4.59p becomes 

7min ^ ^ ■^^f9>^fh- (4-61) 

^' f^G f 
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The niajorization relation (j4.57p restricts the possible allowed Schmidt coefhcients for the blank 
state and can easily be checked numerically, but an analytic expression is difficult to find, since there 
is no simple way of ordering (|4.58p . That is why parts ii) and iii) of the Lemma have their importance, 
since they focus only on the smallest Schmidt coefficient of the blank state. In particular, the bound 
iii) is crucial in deriving the necessity of a maximally entangled blank state for the local cloning of 
qubit and group-shifted qutrit states. 

The proof of the Lemma is rather technical and is presented in Appendix I4.A.2I However, the 
main idea of the proof consists of adding an ancillary system He of dimension D on Alice's side 
and then considering a superposition X]/6G •v/^l'/'/)"^^ ^ 1'/')"^ ® \f}^ ^ that will be mapped by 
the deterministic separable operation to an ensemble {pm: l^m.out)'*"^'''^}, with \'^m,out)^°''^''^ = 
E/eG '^"^'"^ -v/M/IV^/)"^^ ® ® 1/)"^, and wc have used the fact discovered above that pmf ~ p,n, 

independent of /. The average Schmidt vector of the output ensemble over the AaE/Bb cut has 
to majorize the input Schmidt vector, see [GGOSj . and this yields i). Parts ii) and iii) are direct 
implications of i). 

4.5.2 Qubits and Qutrits 

When D = 2 or 13 = 3, one can easily show that the minimum in (|4.6ip is exactly one, and therefore 
Theorem 4.11 (Necessity of maximally entangled blank). The following must hold. 

i) A maximally entangled state of Schmidt rank 2 is the minimum required resource for the local 
cloning of 2 entangled qubit states. 

ii) A maximally entangled state of Schmidt rank 3 is the minimum required resource for the local 
cloning of 3 group-shifted entangled qutrit states. 

The proof of both i) and ii) follows easily from Lemma |4.10[ iii), by applying the rearrangement 
inequality to (|4.6ip , and is presented in Appendix I4.A.3I 

When £> = 2, or when D = 3 and all phases e"'^" = 1, an explicit protocol for cloning these states 
exists |KE06j (alternatively, see the proof of our Theorem l4.8p . and therefore Theorem 14.111 becomes 
a necessary and sufficient condition for the local cloning of such states. In particular, together with 
Theorem 14. 7[ it provides a complete solution to the problem of local cloning when D = 2. 

4.5.3 D > 3, finite gap in the necessary entanglement 

For D > 3, preliminary numerical studies indicate that the minimum (|4.6ip in Lemma 14.101 iii) is 
often equal to one, with few exceptions. It might be the case that a better choice of {fif} in (j4.59p 
of Lemma |4.10[ ii) may provide the 1/D lower bound, but we were unable to prove this. 

However, for any set of group-shifted states, we can prove that there is a rather sizeable gap 
between the entanglement needed in the blank state and the entanglement of the states of S, as 
stated by the following theorem. 

Theorem 4.12 (Finite gap). Let S = {|V'/)^^}/eG be a set of D group-shifted full Schmidt rank 
bipartite orthogonal entangled states on "Ha ^'Hb cls defined by (|4.54p and considered to be not 
maximally entangled. If the local cloning of S using a blank state \<f>)'^^ is possible by a separable 
operation, then the entanglement of the blank state has to be strictly greater that the entanglement 
of the states in S, often by a wide margin. Specifically, 

Eil^r") ^ H{{qr}) > Ei\^jf)^^),yf e G, (4.62) 

where E{-) denotes the entropy of entanglement andH{{qr}) is the Shannon entropy of the probability 
distribution {qr}, qr -"^ J2 fee ^f-^fr ' EreG*' = ^■ 

The proof follows by setting i^ij = 1/D in Lemma 14.101 i) , but is rather long and is presented in 
Appendix I4.A.4I 
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4.6 Conclusion and open questions 

We have investigated the problem of local cloning of a set S of bipartite D x D entangled states 
by separable operations, at least one of which is full Schmidt rank. We proved that all states in S 
must be full rank and that the maximal set of clonable states must be generated by a finite group 
G of order N, the number of states in this maximal set, and then we showed that has to divide 
D exactly. We further proved that all states in S must be equally entangled with respect to the 
G-concurrence measure, and this implied that any two states in S must either share the same set of 
Schmidt coefficients or otherwise be incomparable under majorization. 

We have completely solved two important problems in local cloning. For D = 2 (entangled 
qubits), we proved that no more than two states can be locally cloned, and that these states must be 
locally-shifted. We showed that a two-qubit maximally entangled state is a necessary and sufficient 
resource for such a cloning. In addition, we provided necessary and sufficient conditions when the 
states are maximally entangled, valid for any dimension D, showing that the states must be group- 
shifted, and then we also provided an LOCC protocol that clones such a set of states. 

We have studied in detail the local cloning of partially entangled group-shifted states and provided 
an explicit protocol for local cloning of such states with a maximally entangled resource. For D = 3 
(entangled qutrits) we showed that a maximally entangled blank state is also necessary and sufficient, 
whereas for > 3 we proved that the blank state has to be strictly more entangled than any state 
in S, often by a sizeable amount. 

The necessary form of the clonable states for D > 2 remains an open problem. One might guess 
that the states have to be of a group-shifted form, but a proof of such a claim is not presently 
available. Although we proved the necessity of a maximally entangled resource for the D = 2 case 
and for group-shifted states in the D ~ 3 case, in higher dimensions it is still not clear if a maximally 
entangled state of Schmidt rank D is always necessary. Finally it would be of interest to investigate 
the local cloning of less than full Schmidt rank states, a problem that is likely to bring in additional 
complications, such as the possibility of first distinguishing amongst the states in S while preserving 
the states intact jCoh07| . and then once the state is known, the cloning becomes straightforward 
with a blank state having Schmidt coefficients that are majorized by those of each of the states in 
S [Nk99l[GG08] . 

4. A Mathematical proofs 
4.A.1 Proof of Lemma 14.51 

Consider the singular value decomposition of 5, S* = VDU with T) diagonal and positive definite, 
and V and U unitary operators. Using this expression for S in Tf = SL(f)S^^ shows that 

V^TfV ^V{UL{f)U^)V^^, (4.63) 

or with ff = V^TfV and L{J) = UL{f)W, 

ffV^VLif). (4.64) 

Since Tf and L{f) are both unitary, it is not difficult to see from this that each commutes with 
ptp = V'^. That is, 

Vf[Lif)],, = [Lif)l,,Vl (4.65) 

from which we conclude that when Vi ^ Vj, [Tf]ij = = [L{f)]ij. By a judicious choice of U and 
V, we may arrange for I? to be a direct sum of scalar matrices (some may be one-dimensional). 
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That is, 2? ~ (B^a,^!^, and then we see that Tf and L{f) share the same block-diagonal structure, 
with blocks corresponding to this direct sum decomposition of V. 
We also have directly from (|4.64p that 

[ff].,,V,^V.,[L{f (4.66) 

Therefore, when Vj = T>i, [Tf]ij = [L{f)]ij, and we see that the blocks of Tf are identical to those 
of L{f). In other words, we have shown that Tf = L{f) or equivalently, Tf = WL{f)W^ with 
W = VU, completing the proof. 



4. A. 2 Proof of Lemma 14.101 

Proof of i) Let us introduce an ancillary system He oi dimension D on Alice's side and construct 
the superposition 

\^.^^yBabE ^ ^li^f)^^ ® 10)-" ® 1/)^, (4.67) 

with {fif}f^G 3-11 arbitrary set of non- negative real coefficients that satisfy '^ffJ'f = 1- The proof 
is based on the fact that if {ipf)^^ (8 is deterministically mapped to e^'^"^f\ipf)^^ (g) IV'/)"^ 

(see (|4.1ip ). then |^'i„)^-^"''-^ will be deterministically mapped to an ensemble {p™, |^'m,out)^°^^^}, 
where 

|*™,out)^'^^''^ = E e^'^-'Wl^/)^'' ® I'/'/)"' ® I/)""- (4-68) 

Note that this conclusion rests crucially on the fact, discovered in the main text, that p„if = Pm, 
independent of /. 

Let us now wr ite l^'in)^^''''^ in Schmidt for m over the AaE/Bb cut. One has (again we use 

7 = ./-') 

fee \g,h<£G J 

= E e"'^''^y^^7A;7;^l5>^|/^>1/>''® 1/5)^1/^)' 



E E c''''"yM/A7,7/.|/5)-'l/)'' I \hr ® \9)''\h) 
cj.heG \feG 



= E E ^''''■'^^^^l>^fMf9)^\f)'' I W ® \9)''\h)\ (4.69) 

g,h(^G \feG 

where we used the group property of G and replaced g by fg and summation over / by summation 
over / where necessary. The states on the AaE system are orthogonal for different pairs of g, h, 
and therefore (|4.69p represents a Schmidt decomposition, with Schmidt coefficients ctg^h given by 
the squared norm of the states on the AaE system, 

ag.h = jh E M/A/ff- (4.70) 

/6G 

A similar calculation yields for the Schmidt coefficients /3gji of l^'m.out)"^^"'''^ the expression 

l^g.h = E M/^/ff^/'^' (4-71) 
fee 
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independent of m, which means that the average Schmidt vector of the output ensemble under the 
Aa/BbE cut is the same as the Schmidt vector of an individual state \'i'm,out)^^°'''^ ■ 

We have proven in |GG08| that the average Schmidt vector of the output ensemble produced 
by a separable operation acting on a pure state has to majorize the input Schmidt vector, and this 
concludes i). 

Proof of ii) The proof follows as a direct consequence of i). A particular majorization inequality 
imposed by Lemma 14.101 i) requires that the smallest Schmidt coefficients amin and /3,nin have to 
satisfy 

Ckmin -5^ /^min: (4.72) 

where a and /3 were defined in (|4.70p and (|4.7ip . respectively. This is equivalent to 

7min i> . ^ , ■ 

mingeGZ^/ecAiJ-^/s 

The above equation must hold regardless of which set of {/^/} was chosen, hence taking the maximum 
over all possible sets {^/} concludes the proof of ii). 

Proof of iii) Inserting the expression (|4.60p for {/i/} in (|4.73p yields 

ming,,,eGE/eG XT'^^/sV 
7min ^ : (4-74) 



(4-75) 



D g.neu • — ■ Af 
fee J 

where (j4.75p follows from applying the rearrangement inequality to the denominator in (|4.74p , which 
in this case reads as 



/SG J fee J 



4.A.3 Proof of Theorem \4A1] 

Proof of i) In this case the group G is the cyclic group of order 2, and we identify its group elements 
by {0, 1}. We proved in Theorem 14. 71 that the qubit states have to be locally shifted. The minimum 
in (|4.6ip of Lemma I4.10[ iii) becomes explicitly a minimum over 4 quantities that correspond to 
all possible pairings of g, h; a straightforward calculation shows that 3 out of these 4 quantities are 
equal to 1, except for g = ft, = 1, in which case the sum in (|4.75p equals A^/Aq + Ag/Ai. Order the 
A's such that Ao ^ Xi and note that 

^ s$ and (4.77) 
Ao Ai 

Ai Aq. (4.78) 
From the rearrangement inequality applied to (|4.77p and (|4.78p it follows that 

X2 \2 \2 \2 
M Ao Ai 

and hence the minimum in case i) equals 1. 

Proof of ii) Now the group G is isomorphic to the cyclic group of order 3 and again we identify 
its elements by {0, 1, 2}. We order the A's such that Ao ^ Ai ^ A2. The minimum in (|4.75p is now 
taken over 9 possible pairs g, h. Again straightforward algebra shows that most expressions sum up 
to 1, except for the following three cases for which we show that the sum exceeds 1. 
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1. g = h = 1, ioi which the sum in (j4.75p equals A^/Ao + A|/Ai + Aq/A2; 

2. g = h = 2, for which the sum in (|4.75p equals Aj/Ao + Ap/Ai + Aj/A2; 

3. g = l,h — 2 or g — 2,h — 1, for which the sum in (|4.75p equals A1A2/A0 + A2A0/A1 + A0A1/A2. 
Note first that 

Aq Ai A2 

A2^ < Ai^ < Ao^ and (4.81) 

A1A2 ==; A2A0 AoAi. (4.82) 

From the rearrangement inequality applied to (|4.80p and ()4.8ip it follows that 

—Xl + — A2 + — Aq ^ 

Ao Ai A2 

^ + 1a? + ^A^ = 1, (4.83) 



Ao Ai A 



which proves case 1, and 



2 



— A2 + — Aq + —Xj > 

AO Al A2 

^ + + ^A^ = 1, (4.84) 



Ao Al A2 

which proves case 2. 

Next apply the rearrangement inequality to (|4.80p and (|4.82p to get 

^(AiAa) + ^(AsAo) + ^(AqAi) 

Ao Al A2 



^ i^AoAi + -!-AiA2 + t!-AoA2 = 1 (4.85) 
Ao Al A2 



and this proves case 3. 



4.A.4 Proof of Theorem 

By setting Hf ^ 1/D in Lemma [4.101 i), for all / £ G, the majorization relation (|4.57p reads as 

-^7x1^/3, (4.86) 

where {1/D)j x 1 represents a component vector with components ^h/D, each component 
repeated D times; here 7 is the Schmidt vector of the blank state The components 

I3g,h of P are given by 

Pa,h ^fa^ff^ = E ^/^/s'^- (4.87) 

fee fee 

Note that it is also the case that /3 has D components each repeated D times, so the majorization 
relation (|4.86p implies a majorization relation between 2 D-component vectors 

7 -< (f, (4.88) 



4. LOCAL CLONING OF ENTANGLED STATES BY SEPARABLE OPERATIONS 



44 



where the r-th component of q is given by 

D ■ l3g^h\gh=r^J2^f^f'- ("^-^^^ 
fee 

Note that both 7 and q are normahzed probabihty vectors. Since the Shannon entropy is a Schur- 
concave function, (|4.88p impHes at once that 

EQcj^r") ^ HHqr}). (4.90) 

We now show that the second inequahty in (|4.62p is strict. First we will prove that the ordered 
vector of probabilities g-^ with components defined in (j4.89p and decreasing magnitudes of entries 
down its column, is majorized by A^, the ordered vector of the A/, 

-< A^. (4.91) 

Since the Shannon entropy is not just Schur-concave, but strictly Schur-concave, this will imply at 
once that 

H{{qr}) ^ H{{Xf}) = Ei\4,f)^B), V/ e G, (4.92) 

with equality if and only if q^ equals A"^ (or, equivalently, if and only if the unordered vector q is 
the same as A up to a permutation). One can see that q is not a permutation of A unless all A's are 
equal, case that we exclude. Hence, once we show the majorization condition (|4.9ip holds, the proof 
will be complete. 

We will actually show that A"^ majorizes every vector q of the g^'s no matter how q is ordered. 
Denote by Sn, with \Sn\ = n and n = I,-- - ,D — 1, the subset consisting of those elements f E G 
such that A/ is one of the largest n of the A's. Then, we need to show that for each n, 

ses„ ges„ ges„ fee 

where cr is an arbitrary permutation of the group elements. Since A/ = 1, this is equivalent to 

> 0. (4.94) 



E^i 

/6G 



E -^9 ^ E ^f'^ia) 
aes„ ges,. 



However, given the way we have defined Sn, it is always true that the quantity in square brackets 
is non-negative. The reason is that the first term in this quantity is the sum of the n largest of the 
A's. Therefore the second term, which is also a sum of n of the A's, cannot possibly be greater than 
the first. In fact, it is clear that for general sets of Schmidt coefficients {A/}, the quantity in square 
brackets will not be particularly small, implying that the gap between the required entanglement of 
the blank state and the entanglement of the states in S will be sizable. This ends the proof. 



4.B |G| < D case 

In the main body of the current chapter, we restricted our consideration to the \G\ ~ D case. All 
of our results remain valid also when \G\ < D, with minor modifications. Briefly, when \G\ < D, 
t'j"^^ is a direct sum of rit = D /\G\ copies of L{f), and the following Theorems/Lemmas have to be 
modified accordingly. 
Theorem 
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Since Lemma 14.51 holds for any two unitary representations, it will hold when the regular repre- 
sentation £(/) is replaced by a direct sum of a number of copies of L{f). In this case, the maximally 
entangled group-shifted states (|4.35p and (|4.36p of Theorem 14.61 have the form 

-■ "t 

IV-/)^^ = T^E E \f9.ri)%,n)^, (4.95) 

or 

IV'/)^^ = ^E E \9.n)^\f9,n)^. (4.96) 



respectively. Here the states {\g,n)}gi^Q_n=i,...,nt are an orthonormal basis, {g,n\h,m) = Sg,h5n,m- 
The symbols f,g& G label the group elements and m, n = 1, . . . , label the copies of the regular 
representation. 
Theorem HTSl 

When the family of partially entangled group-shifted states (|4.45p is replaced by 

nt 

\^f)^'' = E E V^l5,")^l/5,^)'' (4.97) 

n=l geG 

and the maximally entangled blank state ()4.46p is modified to 

1 "* 

l<^)"'' = 7^EEl^'"^)"l'^'™)'' (4-98) 



the local cloning protocol of Theorem 14.81 continues to work, provided that 
1. The controlled-group unitary (|4.47p is replaced by 



E E 1-9' ^-9' "I ® ^^^'^ 

n=l geG 

Pa=J2Y.\9h,m){h,m\. (4.99) 

m=l heG 

2. The measurement (|4.49p Alice performs is changed to 

= E 7?^ 1 u7? E m) {h, m\ . (4.100) 

m=l Kl^keG ^k,m) 

where the factor involving the sum over k is needed to insure that this set of measurement 
operators corresponds to a complete measurement. 

3. Finally the unitary correction (|4.5ip Alice and Bob perform is modified to 

nt 

Qr =^^\hr,m) {h,m\. (4.101) 



m=l heG 



Lemma 14.101 

First the blank state has to be modified to 



>'^'' = 7=E EVt^I'^'™)"!^"^)'- (4-102) 

^ ^ m=l heG 
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Next wc follow the line of thought in Appendix I4.A.21 Even though there are only \G\ < D states 
in the clonablc set S, we still use a D dimensional ancillary system He on Alice's side, with a basis 
now given by {\f,n)^}f^G,n=i,...,nf Restricting to an ancillary system of dimension \G\ leads to 
unnecessary complications, since the rearrangement inequality can no longer be applied in part ii) 
to obtain iii). 

We consider again an input superposition 

E E VJ^l^f)^"" ® l^r' ® I/' (4-103) 

n=l /eG 

and look at the Schmidt vector of the output ensemble produced by the separable operation acting 
on (|4.103p . where {/i/,n} is an arbitrary set of coefficients satisfying X]"Li S/eG ^/.s ~ ^- then 
have: 

i) The majorization condition a ~< /3 corresponding to (|4.57p holds, provided the vectors a and 
/3 in (|4.58p are redefined as 

= 7/i,m 2^ 2^ f^J,s^fg.n, 
s=l f^G 

nt 

/^;T=EE/^7,s^/^'."V.'- (4.104) 

s=i fee 



ii) The smallest Schmidt coefficient 7inin of the blank has to satisfy 

min™,„ming,/ieGEr4iE/GGM7,^^/ff,nA//>,m ,,,n^. 

7i„i„ ^ max : ; ' . (4.105) 

{m/,.} mm„mmggG L,s=iT,feGf^7.s^fg,n 



iii) A good choice of {/i/.s} is given by /i/.s = 1/ Aj^ (ig nore the normalization, since ^f^s appears 
both on the numerator and denominator of (|4.105|) V Then (|4.105p becomes 



1 1 
7min ^ TTTrmin min - — Xfg.nXfh.m- (4.106) 

I I " s=l fee 



Theorem \4A2[ 

Theorem 14.121 still provides a finite gap between the entanglement needed in the blank state 
and the entanglement of group shifted states ()4.97p . The proof follows the same ideas as before, by 
setting fif^s = for all / S G and s = 1, . . . , rii in the majorization relation of the "modified" 

Lemma [4.101 1) above. 
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5 

Quantum error correcting codes 
using qudit graph states 

5.1 Introduction 

Quantum error correction is an important part of various schemes for quantum computation and 
quantum communication, and hence quantum error correcting codes, first introduced about a decade 
ago [Sho95 . K L97. .Ste96| have received a great deal of attention. For a detailed discussion see Ch. 10 
of |NCOO| . Most of the early work dealt with codes for qubits, with a Hilbert space of dimension 
D = 2 , but qudit codes with D>2 have also been studied jRaiQQi HMTl ISWO 1 [ [8^^! [Schal IGBR04[ 
lAKPj . They are of intrinsic interest and could turn out to be of some practical value. 

Cluster or graph states, which were initially introduced in connection with measurement based 
or one-way quantum computing |RB01| . are also quite useful for constructing quantum codes, as 
shown in [SWOl, Schb, Scha] in a context in which both the encoding operation and the resulting 
encoded information are represented in terms of graph states. In the present chapter we follow [Schaj 
in focusing on qudits with general D, thought of as elements of the additive group Zjj oi integers 
mod D. However, our strategy is somewhat different, in that we use graph states and an associated 
basis (graph basis) of the n-qudit Hilbert space in order to construct the coding subspace, while not 
concerning ourselves with the encoding process. This leads to a considerable simplification of the 
problem along with the possibility of treating nonadditive graph codes on exactly the same basis as 
additive or stabilizer codes. It also clarifies the relationship (within the context of graph codes as we 
define them) of degenerate and nondegenerate codes, though in this chapter we focus mainly on the 
latter. The approach used here was developed independently in [CSSZ09| and jYCOj for D = 2, and 
in [HTZ+OS] for D > 2; thus several of our results are similar to those reported in these references. 

Following an introduction in Sec. 15.21 to Pauli operators, graph states, and the graph basis, as 
used in this chapter, the construction of graph codes is the topic of Sec. 15.31 In Sec. I5.3.1l we review 
the conditions for an ((n, K, 5)) d code, where n is the number of carriers, K the number of codewords 
or dimension of the coding space, 5 the distance of the code, and D the dimension of the Hilbert 
space of one qudit. We also consider the distinction between degenerate and nondegenerate codes. 
Our definition of graph codes follows in Sec. 15.3.21 and the techniques we use to find nondegenerate 
codes, which are the main focus of this chapter, are indicated in Sec. I5.3.3[ while various results in 
terms of specific codes are the subject of Sec. 15.41 

In Sec. 15.4.21 we show how to construct graph codes with 5 = 2 that saturate the quantum 
Singleton (QS) bound for arbitrarily large n and -D, except when n is odd and D is even, and we 
derive a simple sufficient condition for graphs to yield such codes. For n odd and D = 2 we have an 
alternative and somewhat simpler method of producing nonadditive codes of the same size found in 
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|SSW07) . For both D = 2 and D ~ 3 we have studied nondcgenerate codes on sequences of cycle 
and wheel graphs, in Sees. 15.4.31 and 15.4.41 These include a number of cases which saturate the QS 
bound for (5 = 2 and 3, and others with 5 = 3 and 4 which are the largest possible additive codes for 
the given n, D, and 5. Section [5.4.41 contains results for a series of hypereube graphs with n 4, 8, 
and 16, and in particular a ((16, 128,4))2 additive code. 

In Sec. 15.51 we show that what we call G-additive codes are stabilizer codes (hence "additive" 
in the sense usually employed in the literature), using a suitable generalization of the stabilizer 
formalism to general D. In this perspective the stabilizer is a dual representation of a code which 
is equally well represented by its codewords. The final Sec. 15.61 has a summary of our results and 
indicates directions in which they might be extended. 

5.2 Pauli operators and graph states 
5.2.1 Pauli operators 

Let j = 0, 1, ... D — 1 be an orthonormal basis for the D-dimensional Hilbert space of a qudit, 

and define the unitary operators 

D-l D-l 

Z=Y,u^^\j){j\, X=Y,\j){j®l\, (5.1) 

with © denoting addition mod D. They satisfy 

Z^=I = X^, XZ = ujZX, uj ■=(?''''I° . (5.2) 

We shall refer to the collection of operators {X^'Z'^}, ^, = 0, 1, . . . , D — 1, as (generalized) 
Pauli operators, as they generalize the well known /, X, Z, XZ (= —lY) for a qubit. Together they 
form the Pauli basis of the space of operators on a qudit. 

For a collection of n qudits with a Hilbert space H = 'Hi®'H2® • ■ ■ Tin we use subscripts to 
identify the corresponding Pauli operators: thus Zi and Xi operate on the space T-Li of qudit I. An 
operator of the form 

P = Lo^X^' Z\- Zl^ ■ ■ ■ Zl- , (5.3) 

where A, and /z; and vi for 1 < ^ < tt., are integers in the range to D — 1, will be referred to as 
a Pauli product. If [ii and vi are both 0, the operator on qudit I is the identity, and can safely be 
omitted from the right side of (|5.3p . The collection Q of all operators P of the form (|5.3p with 
A = 0, i.e., a prefaetor of 1, forms an orthonormal basis of the space of operators on "H with inner 
product (A, B) = Z3^"Tr(^'''i3); we call it the (generalized) Pauli basis Q. 

If P and Q are Pauli products, so is PQ, and hence the collection V of all operators of the form 
(|5.3|) for n fixed form a multiplicative group, the Pauli group. While V is not Abclian, it has the 
property that 

PQ = uj^QP, (5.4) 

where ^ is an integer that depends on P and Q. (When D = 2 and w = — 1 it is customary to also 
include in the Pauli group operators of the form (|5.3I) multiplied by i. For our purposes this makes 
no difference.) 

The base of an operator P of the form (|5.3p is the collection of qudits, i.e., the subset of 
{1,2,... n}, on which the operator acts in a nontrivial manner, so it is not just the identity, which 
is to say those j for which either /ij or Vj or both are greater than 0. A general operator R can 
be expanded in the Pauli basis Q, and its base is the union of the bases of the operators which 
are present (with nonzero coefficients) in the expansion. The size of an operator R is defined as 



^ See IHDM05] for a list of references to work that employs operators of this type. 
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the number of qudits in its base, i.e., the number on which it acts in a nontrivial fashion. For 
example, the base of P = uj^X^Xj^Z^ (assuming £> > 3) is {1,4} and its size is 2; whereas the size 
of i? = + 0.5X3^1^3 + iX^ is 4. 

For two distinct qudits I and m the controlled- phase operation C;„i on 'Hi(E}'Hm, generahzing the 
usual controUed-phase for qubits, is defined by 

D-lD-l D-l 

cim = Y.Y. ® \^)^^\ = E ® (5-5) 

i=0 fc=0 3=0 

Of course, Cim = Cmi, and it is easily checked that (Cim)^ = I- It follows from its definition that 
Cim commutes with Zi and and thus with Zp for any qudit p. 

5.2.2 Graph states 

Let G = (V, E) be a graph with n vertices V, each corresponding to a qudit, and a collection E of 
undirected edges connecting pairs of distinct vertices (no self loops). Multiple edges are allowed, as 
in Fig lS.ll for the case of D = 4, as long as the multiplicity (weight) does not exceed D ~ 1, thus at 
most a single edge in the case of qubits. The Im element F;™ = F,„; of the adjacency matrix T is 
the number of edges connecting vertex I with vertex m. The graph state 

|G) =W|G") =W (!+)«"), (5.6) 

is obtained by applying the unitary operator 

{l,m}eE 

to the product state 

|G°) (5.8) 

where 

D-l 

\+):=D-'/'J2\j) (5.9) 
3=0 

is a normalized eigenstate of X, with eigenvalue 1. In (|5.7p the product is over all distinct pairs of 
qudits, with {Cim)^ = I when / and m are not joined by an edge. Since the G;,„ for different I and 
m commute with each other, and also with Zp for any p, the order of the operators on the right side 
of (|5.7p is unimportant. 

Given the graph G we define the graph basis to be the set of D" states 

|a) := |ai,a2,...,a„) = Z^|G) 

^Z'^^Zr---Z:-\G) (5.10) 

where a = (ai,...a„) is an 71-tuplc of integers, each taking a value between and D — I. The 
original graph state |G) is |0, 0, . . . , 0) in this notation. That this collection forms an orthonormal 
basis follows from the fact that the Zp operators commute with the Cim operators, so can be 
moved through the unitary U on the right side of (|5.6p . As the states Z'^\+), < ly < D ~ 1, 
are an orthonormal basis for a single qudit, their products form an orthonormal basis for n qudits. 
Applying the unitary U to this basis yields the orthonormal graph basis. The n-tuple representation 
in (j5.10p is convenient in that one can define 



|a®b) |ai © &i, 02 © 62, •■ 
|ja) := \jai,ja2, ■ ■ -Jan), 



(5.11) 
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Figure 5.1: Action of X and X^ 



on graph state {D — A). 



where j is an integer between and D — 1, and arithmetic operations are mod D. 

One advantage of using the graph basis is that its elements are mapped to each other by a Pauh 
product (up to powers of w), as can be seen by considering the action of Zi or Xi on a single qudit. 
The result for Zi follows at once from (|5.10p . And as shown in Add. 15. Al and illustrated in Fig. 15. 1[ 
the effect of applying Xi to \G) is the same as applying (Zm)^'"^ to each of the qudits corresponding 
to neighbors of / in the graph. Applying these two rules and keeping track of powers of uj resulting 
from interchanging Xi and Zi , sec (|5.2p , allows one to easily evaluate the action of any Pauli product 
on any |a) in the graph basis. 

5.3 Code construction 
5.3.1 Preliminaries 

Consider a quantum code corresponding to a X-dimensional subspace, with orthonormal basis {|cg)}, 
of the Hilbert space H oi n qudits. When the Knill-Laflamme |KL97) condition 



is satisfied for all q and r between and K—1, and every operator Q onH such that 1 < size{Q) < S, 
but fails for some operators of size i5, the code is said to have distance S, and is an ((rt, K, S)) d code; 
the subscript is often omitted when D = 2. (See the definition of size in Sec. 15. 2.1] The only operator 
of size is a multiple of the identity, so (|5.12p is trivially satisfied.) A code of distance 5 allows the 
correction of any error involving at most [{5 — 1)/2J qudits, or an error on 5 — 1 (or fewer) qudits 
if the location of the corrupted qudits is already known (e.g., they have been stolen). 

It is helpful to regard (j5.12p as embodying two conditions: the obvious off-diagonal condition 
saying that the matrix elements of Q must vanish when r ^ q; and the diagonal condition which, 
since f{Q) is an arbitrary complex-valued function of the operator Q, is nothing but the requirement 
that all diagonal elements of Q (inside the coding space) be identical. The off-diagonal condition 
has a clear analog in classical codes, whereas the diagonal one does not. Both must hold for all 
operators of size up to and including (5—1, but need not be satisfied for larger operators. 

In the coding literature it is customary to distinguish nondegenerate codes for which ,f{Q) ~ 
for all operators of size between 1 and (5—1, i.e., for all q and r 



(cJQIO =/(Q)<5, 



qr 



(5.12) 



{cq\Q\cr) for 1 < size(Q) < 6, 



(5.13) 



5. QUANTUM ERROR CORRECTING CODES USING QUDIT GRAPH STATES 51 



and degenerate codes for which f{Q) 7^ for at least one Q in the same range of sizes. See p. 444 of 
jNCOOj for the motivation behind this somewhat pecuhar terminofogy when 6 is odd. In this chapter 
our focus is on nondegenerate codes. For the most part they seem to perform as well as degenerate 
codes, though there are examples of degenerate codes that provide a larger K for given values of n, 
S, and D than all known nondegenerate codes. Examples are the ((6, 2, 3))2 El and ((25,2,9))2 codes 
mentioned in (CRSS98j . 

5.3.2 Graph codes 

When each basis vector \cg) is a member of the graph basis, of the form (|5.10p for some graph G, 
we shall say that the corresponding code is a graph codes associated with this graph. As noted in 
Sec. 15.11 this differs from the definition employed in jSWOli ISchbl ISchaj . but agrees with that in 
more recent D = 2 studies |YCO[ ICSSZOQj , because we do not concern ourselves with the processes 
of encoding and decoding. In what follows we shall always assume 5 > 2, since 5 = 1 is trivial. As 
the left side of (|5.12p is linear in Q, it sufBces to check it for appropriate operators drawn from the 
Pauli basis Q as defined in Sec. 15.2.11 It is helpful to note that for any Q g Q, any pair \cq) and 
\cr) of graph basis states and any n-tuple a, 

(cg © a| Q\cr e a) = (cq| Z-^QZ^\cr) 

= LO^{Cg\Q\Cr) (5.14) 

for some integer ji depending on Q and a; see (|5.10p . (|5.1ip and (|5.4p . Therefore, if (j5.12p is 
satisfied for sonic Q and a collection {|cq)} of codewords, the same will be true for the same Q and 
the collection {\cq ® a)} (with an appropriate change in f{Q))- Thus we can, and hereafter always 
will, choose the first codeword to be 

|co> = |0,0,...,0) = |G). (5.15) 

Analogous to Hamming distance in classical information theory we define the Pauli distance A 
between two graph basis states as 

A(c„c,) = A(|c,),|c,)) := min{size(g) : (c,|Q|c,) ^0}, (5.16) 

where it suffices to take the minimum for Q Cz Q, the Pauli basis. (Ket symbols can be omitted from 
the arguments of A when the meaning is clear.) Also note the identities 

A{Cg, Cr) = A(C,., Cq) = A{Cq a, © a) 

= A(co,c^ec5), (5.17) 

where a is any n-tuple, and © means difference mod D, see (|5.1ip . The second equality is a 
consequence of (|5.14p . Note that if in (j5.16p we minimize only over Q operations which are tensor 
products of Z's (no X's), A is exactly the Hamming distance between the n-tuples Cq and c^, see 

For the case q = r, where (|5.16|) gives (for Q ~ I), we introduce a special diagonal distance 
A' which is the minimum size of the right side of (|5.16p when one restricts Q to be an element of 
Q of size 1 or more. The diagonal distance does not depend on the particular value of g = r, but 
is determined solely by the graph state |G) — see (|5.14p with r = q — and thus by the graph G. This 
has the important consequence that if we consider a particular G and want to find the optimum 

^ While there seems to be no proof that the ((6, 2, 3))2 degenerate code has a larger K than any nondegenerate code 
with n = 6 and 5 = 3, some support comes from the fact that we performed an exhaustive search of all graphs with 
6 vertices and did not find a nondegenerate graph code with 5 = 3 and K > I. But the notion that this degenerate 
code is superior to nondegenerate codes is undercut by the observation that the well known nondegenerate ((5, 2, 3))2 
code uses only 5 instead of 6 qubits to achieve equal values of K and 5. 
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codes for a given S that is no larger than A', the collection of operators Q G Q for which (|5.12p 
needs to be checked will all have zero diagonal elements, f{Q) = 0, and we can use (|5.13p instead of 
(j5.12p . In other words, for the graph in question and for S < A', all graph codes are nondegenerate, 
and in looking for an optimal code one need not consider the degenerate case. Our computer results 
in Sec. 15.41 are all limited to the range 6 < A' where no degenerate codes exist for the graph in 
question. Any code with S > A' will necessarily be degenerate, since there is at least one nontrivial 
Q for which (|5.12|) must be checked for the diagonal elements. 

A code is G-additive (graph-additive) if given any two codewords \cq) and jc^) belonging to the 
code, \cq © c^) is also a codeword. As shown in Sec. 15.51 this notion of additivity implies the code is 
additive in the sense of being a stabilizer code. For this reason, we shall omit the G in G-additive 
except in cases where it is essential to make the distinction. Codes that do not satisfy the additivity 
condition are called nonadditive. The additive property allows one to express all codewords as 
"linear combinations" of k suitably chosen codeword generators. This implies an additive code must 
have K = Z?'", r an integer, whenever D is prime. We will see an example of this in Sec. 15.41 for 
D = 2. 

The quantum Singleton (QS) bound |KL97| 

n > log^ K + 2{S - 1) or K < D"-2('5-i) (5.18) 

is a simple but useful inequality. We shall refer to codes which saturate this bound (the inequality 
is an equality) as quantum Singleton (QS) codes. Some authors prefer the term MDS, but as it is 
not clear to us how the concept of "maximum distance separable," as explained in (MS77j . carries 
over to quantum codes, we prefer to use QS. 

5.3.3 Method 

We are interested in finding "good" graph codes in the sense of a large K for a given n, S, and 
D. The first task is to choose a graph G on n vertices, not a trivial matter since the number of 
possibilities increases rapidly with n. We know of no general principles for making this choice, 
though it is helpful to note, see App. 15. Al that the diagonal distance A' cannot exceed 1 plus the 
minimum over all vertices of the number of neighbors of a vertex. Graphs with a high degree of 
symmetry are, for obvious reasons, more amenable to analytic studies and computer searches than 
those with lower symmetry. 

Given a graph G and a distance 5, one can in principle search for the best nondegenerate code 
by setting |co) = \G), finding a |ci) with A(co,Ci) > 6, after that |c2) with both A(co,C2) > S and 
A(ci,C2) > i5, and so forth, until the process stops. However, this may happen before one finds 
the largest K, because a better choice could have been made for \cq) at some point in the process. 
Exhaustively checking all possibilities is rather time consuming, somewhat like solving an optimal 
packing problem. 

In practice what we do is to first construct a lookup table containing the D" — 1 Pauli distances 
from \G) to all of the other graph basis states, using an iterative process starting with all Q G Q 
of size 1, then of size 2, etc. This process also yields the diagonal distance A'. As we are only 
considering nondegenerate codes, we choose some (5 < A', so that (|5.13|) can be used in place of 
()5.12|) . and use the table to identify the collection S of all graph basis states with a distance greater 
than or equal to S from |co) = \G). If S is empty there are no other codewords, so K = 1. However, 
if S is not empty then K is at least 2, and a search for the optimum code (largest K) is carried out 
as follows. 

We produce a graph S (not to be confused with G) in which the nodes are the elements of S, and 
an edge connects two nodes if the Pauli distance separating them — easily computed from the lookup 
table with the help of (|5.17|) — is greater than or equal to 5. An edge in this graph signifies that the 
nodes it joins are sufficiently (Pauli) separated to be candidates for the code, and an optimal code 
corresponds to a largest complete subgraph or maximum clique of S. Once a maximum clique has 
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been found, the corresponding graph basis states, including |co), satisfy (|5.13p and span a coding 
space with the largest possible K for this graph G and this S. 

The maximum clique problem on a general graph is known to be NP-complete jGJ79j and hence 
computationally difficult, and we do not know if S has special properties which can be exploited to 
speed things up. We used the relatively simple algorithm described in |CP90| for finding a maximum 
clique, and this is the most time-consuming part of the search procedure. 

The method just described finds additive as well as nonadditive codes. In fact one does not know 
beforehand whether the resultant code will be additive or not. If one is only interested in additive 
codes, certain steps can be modified to produce a substantial increase in speed as one only has to 
find a set of generators for the code. 

5.4 Results 

5.4.1 Introduction 

Results obtained using methods described above are reported here for various sequences of graphs, 
each sequence containing graphs of increasing n while preserving certain basic properties. We used 
a computer search to find the maximum number K of codewords for each graph in the sequence, 
for distances S < A' and for D = 2 or 3, qubits and qutrits, up to the largest number n of qudits 
allowed by our resources (running time). Sometimes this revealed a pattern which could be further 
analyzed using analytic arguments or known bounds on the number of codewords. 

In the ease of distance S = 2 we can demonstrate the existence of QS codes for arbitrarily large 
values of n and D, except when n is odd and D is even, see Part A. In the later subsections we 
report a significant collection of Z? = 2 and 3 codes for S = 2, 3, and 4, including QS codes; codes 
which are the largest possible additive codes for that set of n, D and S; and a new ((16, 128,4))2 
additive code. 

Tables show the K found as a function of other parameters. The meaning of superscripts used 
in the tables is given below. 

• a — Indicates the maximum clique search was terminated before completion. This means the 
code we found might not be optimal, i.e. there might be another code with larger K for this 
graph. We can only say the code is maximal in the sense that no codeword can be added 
without violating (|5.13p . Absence of this superscript implies no code with a larger K exists 
for this d and this graph, either because the program did an exhaustive search, or because K 
saturates a rigorous bound. 

• & - Indicates a nonadditive code. Codes without this superscript are additive. 

• c - Indicates a QS code, one where K saturates the Singleton bound ()5.18|) . 

• d - Indicates this is not a QS code, but the largest possible additive (graph or other) code for 
the given n, 6 and D, This follows from linear programming bounds in |Gra07j for D = 2 and 
[KKKS06] for D = 3, along with the fact. Sec. I5.3.2[ that for an additive code, K must be an 
integer power of D when D is prime. A larger nonadditive code for this graph might still be 
possible in cases flagged with a as well as d. 

5.4.2 Distance S = 2; bar and star graphs 

It was shown in |CRSS98] that for = 2 one can construct 6 = 2 QS codes for any even n, and 
similar codes for larger D are mentioned, without giving details, in |Rai99| . One way to construct 
graph codes with 5 = 2 is to use the method indicated in the proof, App. 15. B| of the following result. 
Partition theorem. Suppose that for a given D the vertices of a graph G on n qudits can be 
partitioned into two nonempty sets Vi and V2 with the property that for each vertex in Vi the sum of 
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(a) (b) (c) 




(d) (e) 



Figure 5.2: Examples from different graph sequences: (a) bar (odd n), (b) star, (c) cycle, (d) wheel, 
(e) n = 16 hypercube. 



the number of edges (the sum of the multiplicities if multiple edges are present) joining it to vertices 
in V2 is nonzero and coprime to D, and the same for the number of edges joining a vertex in V2 to 
vertices in Vi. Then there is an additive QS code on G with distance S = 2. 

A bar graph is constructed by taking n vertices and dividing them into two collections Vi and V2 , 
of equal size when n is even, and one more vertex in V2 when n is odd, as in Fig. I5.2r a). Next pair 
the vertices by connecting each vertex in Vi by a single edge to a vertex in V2, with one additional 
edge when n is odd, as shown in the figure. (Multiple edges are possible for D > 2, but provide 
no advantage in constructing codes.) When n is even the conditions of the partition theorem are 
satisfied: 1 is always coprime to D. For odd n, the last vertex in Vi has 2 edges joining it to V2, 
which is coprime to D when D is odd. Hence bar graphs yield 5 = 2 QS codes for all n when D is 
odd, and for even n when D is even. 

A star graph, Fig. I5.2f b). has a central vertex joined by single edges to every peripheral vertex, 
and no edges connecting pairs of peripheral vertices. Since the diagonal distance A' is 2, nondegen- 
erate star codes cannot have 5 larger than 2. As in the case of bar codes, one can construct additive 
QS codes for any n when D is odd, and for even n when D is evenlfl . For odd n and D = 2 there 
are nonadditive codes with 

see App. 15. CI for details. Codes with these parameters were discovered earlier by Smolin ct al. 
|SSW07| using a different approach. Computer searches show that for all odd n < 7 star graphs 
cannot yield a K larger than (j5.19p . 



5.4.3 Cycle graphs 

We used computer searches to look for graph codes based on cycle (loop) graphs, Fig. I5.2r c). Ta- 
ble 15.11 shows the maximum number K of codewords for codes of distance 5 = 2 and (5 = 3 for both 
D = 2 qubits and Z? = 3 qutrits. In the qutrit case the best codes were obtained by including one 
double edge (weight 2), as in Fig. I5.2r c). though when n is odd equally good codes emerge with only 
single edges. In the qubit case all edges have weight 1. 

^ We omit the details. In some but not all cases one can use the Partition theorem with Vi and V2 the center and 
the peripheral vertices. Allowing some double edges when D > 2 extends the range of n values where the Partition 
theorem can be employed. 



5. QUANTUM ERROR CORRECTING CODES USING QUDIT GRAPH STATES 55 



Table 5.1: Maximum K for qubit and qutrit cycle graphs. Sec Sec. 15.4.11 for detailed meaning of 
superscripts. 
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a Non-exhaustive search 
b Nonadditive code 

c Code saturating Singleton bound (|5.1Sp 
d Largest possible additive code 



The D = 2 entries in Table include for n = 5 the well known ((5,2,3))2, the nonadditive 
((5, 6, 2))2 presented in [RHSS97! , and, for larger n, a ((9, 12, 3))2 code similar to that in |YCLO08j 
and the ((10, IS, 3))2 of ^CSSZ09) based upon the same graph. 

The D = 3, S = 3 entries are interesting because the QS bound is saturated for 4 < 71 < 10 but not 
for n = 11. The ((11,3® = 729, 3))3 code we found, the best possible additive code according to the 
linear programming bound in [KKKS06] . falls short by a factor of 3 of saturating the X = 3^ = 2187 
QS bound, and even a nonadditive code based on this graph must have K < 19900 . 

One can ask to what extent the results for (5 = 2 in Table 15.11 could have been obtained, or 
might be extended to larger n, by applying the Partition theorem of Part A to a suitable partition 
of the cycle graph. It turns out — we omit the details — that when D is odd one can use the Partition 
theorem to produce codes that saturate the QS bound for any n, but when D is even the same 
approach only works when n is a multiple of 4. In particular, the ((6, 16, 2))2 additive QS code 
in Table 15.11 cannot be obtained in this fashion since the cycle graph cannot be partitioned in the 
required way. 

5.4.4 Wheel graphs 

If additional edges are added to a star graph so as to connect the peripheral vertices in a cycle, as in 
Fig. l5.2f d). the result is what we call a wheel graph. Because each vertex has at least three neighbors, 
our search procedure, limited to 6 < A', can yield (5 = 4 codes on wheel graphs, unlike cycle or star 
graphs. The construction of (5 = 2 codes for any D is exactly the same as for star graphs, so in 
Table [52] we only show results for (5 = 3 and 4, for both £) = 2 and 3. The ((16, 128,4))2 additive 
code appears to be new, and its counterpart in the hypercube sequence is discussed below. 

* Since the distance 5 = 3 does not exceed the diagonal distance A' = 3 for this graph, a graph code is necessarily 
nondegenerate, see Sec. 15.3.21 and hence the quantum Hamming bound — see p. 444 of |NCOO| — extended to D = 3 
applies, and this yields an upper bound of K < 1990. 
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Table 5.2: Maximum K for qubit and qutrit wheel graphs. Sec Sec. 15.4.11 for detailed meaning of 
superscripts. 
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5.4.5 Hypercube graphs 

Hypercube graphs, Fig. I5.2r e). have a high symmetry, and as n increases the coordination bound, 
App. 15. Al allows A' to increase with n, unlike the other sequences of graphs discussed above. We 
have only studied the D = 2 case, with the results shown in Table 15.31 Those for 5 = 2 are 
an immediate consequence of the Partition theorem: each hypercube is obtained by adding edges 
between two hypercubes of the next lower dimension, and these are the Vi and V2 of the theorem. 
The generators for the ((16, 128,4))2 additive code are given in Table [531 The 2'^ = 128 codewords 
arc of the form, sec (|5.1ip . © a2g2 © ■ ■ ■ argy), where each aj can be either or 1. 

Table 5.3: Maximum K for qubit hypercube graphs. See Sec. 15.4.11 for detailed meaning of super- 
scripts. 
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Table 5.4: Generators of ((16, 128,4))2 additive code for hypercube graph 

Generator Bit notation 

TsD 10,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1) 

|g2) 10,0,0,0,0,0,0,0,0,0,1,1,0,0,1,1) 

Igg) 10,0,0,0,0,0,0,0,1,1,0,0,0,0,1,1) 

|g4) 10,0,0,0,0,0,1,1,0,1,0,0,0,1,0,0) 

Igs) 10,0,0,0,1,1,0,0,0,0,0,1,0,0,0,1) 

Igg) 10,0,1,1,0,0,0,0,0,1,0,0,0,1,0,0) 

|g7) 11,1,0,0,0,0,0,0,0,0,0,1,0,0,0,1) 



5.5 G- Additive codes as stabilizer codes 

The stabilizer formalism introduced by Gottesman in jGotbj for D ^ 2 (qubits) provides a compact 
and powerful way of generating quantum error correcting codes. It has been extended to cases where 
D is prime or a prime power in |AK01[ IBB[ IKKKS06] . In jSchbj stabilizer codes were extended in 
a very general fashion to arbitrary D from a point of view that includes encoding. However, our 
approach to graph codes is somewhat different, see Sec. 15.11 and hence its connection with stabilizers 
deserves a separate discussion. We will show that for any D > 2 a G-additive (as defined near the 
end of Sec. I5.3.2P code is a stabilizer code, and the stabilizer is effectively a dual representation of 
the code. 

The Pauli group V for general n and D was defined in Sec. 15.2.11 Relative to this group we 
define a stabilizer code (not necessarily a graph code) C to he a, K > 1-dimensional subspace of the 
Hilbert space satisfying three conditions: 

CI. There is a subgroup S oiV such that for every T in S and every \^) in C 

T|V) = IV-) (5.20) 

C2. The subgroup S is maximal in the sense that every T in V for which (|5.20p is satisfied for all 
I?/") G C belongs to S. 

C3. The coding space C is maximal in the sense that any ket \'4)) that satisfies (j5.20|) for every T E S 
lies in C. 

If these conditions are fulfilled we call S the stabilizer of the code C. That it is Abelian follows 
from (|5.4p . since for K > there is some nonzero j^/^) satisfying (|5.20p . One can also replace (j5.20p 
with 

T|c,) = \c,) (5.21) 

where the {|cq)} form an orthonormal basis of C. Note that one can always find a subgroup S oiV 
satisfying CI and C2 for any subspace C of the Hilbert space, but it might consist of nothing but 
the identity. Thus it is condition C3 that distinguishes stabilizer codes from nonadditive codes. A 
stabilizer code is uniquely determined by S as well as by C, since S determines C through C3. 

As we shall see, the stabilizers of G-additive graph codes can be described in a fairly simple way. 
Let us begin with one qudit, n = 1, where the trivial graph G has no edges, and the graph basis 
states are of the form {Z^|+)} for c in some collection C of integers in the range < c < D—1. The 
subgroup S oiV satisfying CI and C2 must be of the form {X'^} for certain values ofs, 0<s<Z?— 1, 
belonging to a collection S. This is because Z and its powers map any state Z'^\+) to an orthogonal 
state, and hence T in (|5.2ip cannot possibly contain a (nontrivial) power of Z. Furthermore, since 



X'Z'\+) =uj^'Z''\+), 



(5.22) 
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sec (|5.2p . will leave unchanged only if o;"^^ — 1, or 

cs = (mod D). (5.23) 

Thus for S to satisfy CI, it is necessary and sufficient that (|5.23|) hold for every c G C, as well as 
every s G S*. Further, S = j^'*} is maximal in the sense of C2 only if S contains every s satisfying 
(|5.23p for each c G C. As shown in App. I5.D[ such a collection S must either (depending on C) 
consist of s = alone, or consist of the integer multiples vsi, with z/ = 0, 1, . . . {D / si — 1), of some 
si > that divides D. In either case, S* is a subgroup of the group of integers under addition 
mod D, and indeed any such subgroup must have the form just described. 

We now take up C3. Given the maximal collection S of solutions to (j5.23p . we can in turn ask 
for the collection of C of integers c in the range to D — 1 that satisfy (|5.23p for every s in S. 
Obviously, C" contains C , but as shown in App. 15. Dl C" = C if and only if C is a subgroup of Z/^, 
i.e., C is G-additive. Next note that every T in 5, as it is a power of X and because of (j5.22[l . maps 
every graph basis state to itself, up to a phase. Thus when (and only when) C is G-additive, the 
codewords are just those graph basis states for which this phase is 1 for every T £ S. To check C3, 
expand an arbitrary ji/j) in the graph basis. Then T\ip) = ji/j) for all T G 5 means that all coefficients 
must vanish for graph basis states that do not belong to C. Hence C3 is satisfied if and only if C is 
G-additive. 

The preceding analysis generalizes immediately to rt > 1 in the case of the trivial graph G'^ with 
no edges. A graph code C has a basis of the form {Z'^|G°)} for a collection C of integer n-tuples 
c G Z^, and is G-additive when the collection C = {c} is closed under component-wise addition 
mod D, i.e., is a subgroup of Z^^. Whether or not C is G-additive, the subgroup S V satisfying 
CI and C2 consists of all operators of the form X^ = Xl'- X'l^ ■ ■ ■ with the n-tuple s satisfying 

n 

c • s := ^ cisi = (mod D) (5.24) 
1=1 

for every c G C. Just as for n ~ 1, S cannot contain Pauli products with (nontrivial) powers of Z 
operators. Let S denote the collection of all such s. The linearity of (j5.24p means S is an additive 
subgroup of Z^ . 

One can also regard (|5.24p as a set of conditions, one for every s E S, that are satisfied by certain 
c G ZJ). The set C' of all these solutions is itself an additive subgroup of Z^, and contains C. 
In App. 15. Dl we show that C" = C if and only if C (the collection we began with) is an additive 
subgroup of Z^ , and when this is the case the sizes of C and S are related by 

|C| • IS*! = (5.25) 

Just as for n = 1, any X^ maps a graph basis state for the trivial graph G^ — they are all product 
states — onto itself up to a multiplicative phase, and the same argument used above for n = 1 shows 
that C3 is satisfied for all T G 5 if and only if C is G-additive. 

To apply these results to a general graph G on ti qubits, note that the unitary l/l defined in (|5.7[) 
provides, through (|5.6p and (|5.10p . a one-to-one map of the graph basis states of the trivial G° onto 
the graph basis states of G. At the same time the one-to-one map UPW carries the iS satisfying CI 
and C2 (and possibly C3) for the G° code to the corresponding S, satisfying the same conditions for 
the G code. (The reverse maps are obtained by interchanging Z//''' and U.) Consequently, the results 
obtained for G" apply at once to G, and the transformation allows the elements of the stabilizer 
for the G graph code to be characterized by integer n-tuples s satisfying (|5.24p . Thus we have 
shown that G-additive codes are stabilizer codes, and for these the coding space and stabilizer group 
descriptions are dual, related by (|5.24p : each can be derived from the other. 
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5.6 Conclusion and discussion 

In this chapter we have developed an approach to graph codes which works for qudits with general 
dimension D, and employs graphical methods to search for specific e xamples o f such codes. It is 
similar to the approaches developed independently in jCSSZ091 lYCOi IHTZ+OS] . We have used it 
for computer searches on graphs with a relatively small number n of qudits, and also to construct 
certain families of graphs yielding optimum distance 6 = 2 codes for various values of D and n which 
can be arbitrarily large. It remains a challenging problem to do the same for codes with distance 
S >2. 

In a number of cases we have been able to construct what we call quantum Singleton (QS) codes 
that saturate the quantum Singleton bound |KL97| : these include the S = 2 codes for arbitrarily 
large n and D mentioned above, and also a number of (5 = 3 codes in the case of I? = 3 (qutrits) , see 
Tables [23] and [221 The results for cycle graphs for D = 3 and (5 = 3 in Table [5?T] are interesting in 
that the QS bound is saturated for n < 10, but fails for n = 11, as it must for nondegenerate codes; 
see the discussion in Sec. 15.4.31 Our results are consistent with the difficulty of finding QS codes for 
larger 6 |Gra07j . but suggest that increasing D may help, as observed in jSWOlj . It is worth noting 
that we have managed to construct many of the previously known nonadditive codes, or at least 
codes with the same {{n, K,S))d, using simple graphs. Some other nonadditive codes not discussed 
here, such as the ((10,24,3))2 code in jYCOj . can also be obtained from suitably chosen graphs. 
While all these results are encouraging, they represent only a beginning in terms of understanding 
what properties of graphs lead to good graph codes, and how one might cfSciently construct such 
codes with arbitrarily large n and S, for various D. 

As noted in Sec. I5.3.2[ all graph codes with distance S < A' , where A' is the diagonal distance 
of the graph, are necessarily nondegenerate, and our methods developed for such codes will (in 
principle) find them all. All codes with S > A' are necessarily degenerate codes, and their systematic 
study awaits further work. It should be noted that our extension of graph codes to D > 2 is based on 
extending Pauli operators in the manner indicated in |HDM05j . Though the extension seems fairly 
natural, and it is hard to think of alternatives when D is prime, there are other ways to approach the 
matter when D is composite (including prime powers), which could yield larger or at least different 
codes, so this is a matter worth exploring. 

The relationship between stabilizer (or additive) codes and G-additive (as defined in Sec. I5.3.2|) 
graph codes has been clarified by showing that they are dual representations, connected through 
a simple equation, (|5.24p . of the same thing. One might suspect that such duality extends to 
nongraphical stabilizer codes, but we have not studied the problem outside the context of graph 
codes. Nonadditive codes, which — if one uses our definition. Sec. 15.51 — do not have stabilizers, are 
sometimes of larger size than additive codes, so they certainly need to be taken into account in the 
search for optimal codes. The graph formalism employed here works in either case, but computer 
searches are much faster for additive codes. 



5. A The X-Z rule and related 

X-Z Rule. Acting with an X operator on the i'th qudit of a graph state \G) produces the same 
graph basis state as the action of Z operators on the neighbors of qudit i, raised to the power given 
by the edge multiplicities r,;,„. 

The operator Xi commutes with Cim when i ^ I and i ^ m, but ii i ^ I (or similarly i = m) one 
can show using <\5.5\i and ()5.ip that 

XlClm = ClmZmXl = Z„iClmXi. (5.26) 

That is, an Xi operator can be pushed from left to right through a Cim with at most the cost of 
producing a Z operator associated with the other qudit: \i i = I one gets Z^, if « = m one gets Zi. 
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Since all Z commute with all C, one can place the resulting either to the left or to the right of 

Clm- 

Now consider pushing Xi from the left to the right through U, the product of Cim operators 
defined in (|5.7p . Using (|5.26p successively for those Cim that do not commute with Xi , one sees that 
this can be done at the cost of generating a Z„i for every edge of the graph connecting i to another 
vertex m. Let the product of these be denoted as Z := Y[(i=i m)eE ^m"^ ■ Then, with definition 
()5.6p . we can show 

X,\G) ^ X,U\G") = ZUX,\G") 

= ZU\G°) ^ Z\G) , (5.27) 

which completes the proof of the X-Z Rule. 

For graph codes satisfying (|5.13p . the X-Z Rule leads to the: 
Coordination bound. The diagonal distance A' for a graph G cannot exceed v + l, where v is the 
minimum over all vertices of the number of neighbors of a vertex, this being the number of vertices 
joined to the one in question by edges, possibly of multiplicity greater than 1. 

To make the counting absolutely clear consider Fig. 15.11 where the vertex on the left has 3 
neighbors, and each of the others has 1 neighbor, so that in this case v = 1. To derive the bound, 
apply X to a vertex which has v neighbors. By the X-Z rule the result is the same as applying 
appropriate powers of Z to each neighbor. Let P be this X tensored with appropriate compensating 
powers of Z at the neighboring vertices in such a way that P\G) = \G). The size of P is -|- 1, and 
A' can be no larger. 

Another useful result follows from the method of proof of the X-Z Rule: 
Paulis to Paulis. Let P be a Pauli product (|5.3p . and for U defined in (|5.7p let 

P'=U'^PU, P"=UPU'^. (5.28) 

Then both P' and P" are Pauli products. 

To see why this works, rewrite the first equality as UP' ~ PIA, and imagine pushing each of 
the single qudit operators, of the form X'^' , making up the product P through U from left to 
right. This can always be done, see the discussion following (|5.26p . at the cost of producing some 
additional Z operators, which can be placed on the right side of lA, to make a contribution to P' . 
At the end of the pushing the final result can be rearranged in the order specified in (j5.3p at the 
cost of some powers of w, see (|5.2p . The argument for P" uses pushing in the opposite direction. 



5.B Partition theorem proof 

Given the partition of the n qudits into sets Vi and V2 containing ni and 712 elements, the code of 
interest consists of the graph basis states |c) = |ci,C2, . . . ,c„) satisfying the two conditions 

^ c, EE (mod D) (5.29) 
^ Cj = (mod D) (5.30) 

This code is additive and contains K — x £)"2-i — £)"-2 codewords. (The counting can be 

done by noting that (|5.29p defines a subgroup of the additive group 17^ , and its cosets arc obtained 
by replacing with some other integer on the right side of (|5.29p .) 

We first demonstrate that this code has 6 > 2 hy showing that any Pauli operator, except the 
identity, applied to a single qudit maps a codeword into a graph basis state not in the code. If Z'^ 
for < < I? is applied to a qudit in Vi, the effect will be to replace on the right side of (|5.29p 
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with ly, so this graph state is not in the code. If Xf^, < /i < £> is apphed to a qudit in Vi the 
resuh according to the X-Z Rule, App. 15. Al will be the same as placing Z operators on neighboring 
qudits in V2 (as well as Vi) in such a way that on the right side of (|5.30p is replaced by gfi, where 
g is the total number of edges (including multiplicities) joining the Vi qudit with qudits in V2. But 
as long as g is coprime to Z?, as specified in the condition for the theorem, gfi cannot be a multiple 
of D, and ()5.30|) will no longer be satisfied. The same is true if Z^A"^ is a applied to a qudit in Vi. 
Obviously the same arguments work for Pauli operators applied to a single qudit in V2. Thus we 
have shown that 6 > 2. 

But 5 > 2 is excluded by the QS bound, so we conclude that we have an additive code of 
K = Z?"^^ elements and distance (5 = 2 that saturates the QS bound. 

5.C Construction of qubit star graph codes 

As noted in See. 15. 4. 21 a star graph for n-qubits consists of a central vertex joined by edges to n — 1 
peripheral vertices. Let Vi be the central vertex and V2 the set of peripheral vertices. When n is 
even and D = 2 the conditions of the Partition theorem. Sec. I5.4.2[ are satisfied, and the 6 = 2 code 
constructed in App. 15. Bl consists of the 2"~^ graph basis states with no Z on the central qubit and 
an even number r of Z's on the peripheral qubits, thus satisfying (|5.29l) and (j5.30p . and yielding an 
additive QS code. 

When n is odd the central vertex is connected to an even number n — 1 of vertices in V2, so the 
conditions of the Partition theorem no longer hold. A reasonably large 5 = 2 nonadditive code can, 
however, be constructed by again assuming no codeword has Z on the central qubit, and that the 
code contains all graph basis states with r Z's on the peripheral qubits for a certain selected set R 
of r values. 

The set R must satisfy two conditions. First, it cannot contain both r and r + 1, because applying 
an additional Z to a codeword with r Z's yields one with r + 1, and one cannot have both of them in 
a code of distance 6 = 2. Second, applying X to the central vertex and using the X-Z rule, App. 15. A} 
maps a codeword with r Z's to one with r' = n — 1 — r; hence R cannot contain both r and n—l — r. 
For example, when n = 7(n— 1 = 6 peripheral qubits) the set R = {0, 2, 5} satisfies both conditions, 
as does R = {1, 4, 6}, whereas R = {1, 2, 6} violates the first condition and R = {1, 3, 5} the second. 

By considering examples of this sort, and noting that the number of such graph basis states with 
r Z's is ("!^^) which is equal to („"7-r)' '-'^^ ^^'^^ that for n odd one can construct in this way a 
nonadditive code with 

(n-3)/2 

E 

codewords. 

5.D Solutions to c • s 

Let A be the collection of all n-component integer vectors (i.e., n-tuples) of the form a = (ai, 02, . . . a„), 
< Oj < D ~ 1, with component- wise sums and scalar multiplication defined using arithmetic oper- 
ations mod D. In particular, ^ is a group of order D" under component-wise addition mod D. We 
shall be interested in subsets C and S of A that satisfy 

n 

c-s:=^cisi = (mod D) (5.32) 
1=1 

for all c G C and s G S". Given some collection C, we shall say that S is maximal relative to C if 
it includes all solutions s that satisfy (|5.32p for every c G C. It is easily checked that a maximal S 



\ ^ ) A{n-l)l2) 

= (mod D) 
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is an additive subgroup of A: it includes the zero vector and — s mod D whenever s S S*. A similar 
definition holds for C being maximal relative to a given S. We use |C| to denote the number of 
elements in a set or collection C. 

Theorem. Let C be an additive subgroup of A, and let S be maximal relative to C, i.e., the set 
of all s that satisfy (j5.32p for every c & C. Then C is also maximal relative to S , and 

|C| • IS*! = L>". (5.33) 



The proof is straightforward when Z) is a prime, since Z^j is a field, and one has the usual rules 
for a linear space. The composite case is more difficult, and it is useful to start with n ~ 1: 

Lemma. Let C be a subgroup under addition mod D of the integers lying between and D — 1, 
and S all integers in the same range satisfying 

cs = (mod D) (5.34) 

for every c G C. Then C consists of all integers c in the range of interest which satisfy (|5.34p . and 
\C\ -151=0. 

When C = {0} the proof is obvious, since |C| = 1 and |5| = D. Otherwise, because it is an 
additive subgroup of Zd, C consists of the multiples {/-tci} of the smallest positive integer ci in C, 
necessarily a divisor of -D, when /i takes the values 0, 1, ... si — 1, where si = D/a. One quickly 
checks that all integer multiples s = vs\ of this s\ satisfy (|5.34p and are thus contained in S. But S is 
also an additive subgroup, and s\ is its minimal positive element (except in the trivial case c\ = 1), 
for were there some smaller positive integer s' in S we would have < c\s' < D, contradicting 
(|5.34|) . Similarly there is no way to add any additional integers to C while preserving the subgroup 
structure under addition mod D without including a positive c less than ci, which will not satisfy 
([EM)) for s = si. 

For n > 1 it is helpful to use a generator matrix F, with components Fri, each between and 
D — 1, with the property that c € C if and only if it can be expressed as linear combinations of rows 
of F, i.e., 

ci = Y^ brFri (mod D) (5.35) 

r 

for a suitable collection of integers {br}. This collection will of course depend on the c in question, 
and for a given c need not be unique, even assuming (as we shall) that < br < D — 1. In particular 
the matrix F for which each row is a distinct c in C, with r running from 1 to |C|, is a generator 
matrix. It is straightforward to show that if F is any generator matrix for C, S consists of all 
solutions s to the equations 

n 

J2Frisi = (mod d) for r = 1,2,.... (5.36) 

1=1 

The collections C and S*, vectors of the form (|5.35p and those satisfying ()5.36p . remain the same 
if F is replaced by another generator matrix F' obtained by one of the following row operations: 
(i) permuting two rows; (ii) multiplying (mod D) any row by an invertible integer, i.e., an integer 
which has a multiplicative inverse mod D\ (iii) adding (mod D) to one row an arbitrary multiple 
(mod D) of a different row; (iv) discarding (or adding) any row that is all zeros, to get a matrix of a 
different size. Of these, (i) and (iv) are obvious, and (ii) is straightforward. For (iii), consider what 
happens if the second row of F is added to the first, so that F'^i = F^i except for 

F[i=F,i+F^i {iTiodD). (5.37) 



Then setting 



6'i = 6i, b'^ = b2- bi (mod d), b'l = 6; for / > 3 



(5.38) 
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leads to the same c in (|5.35|) if b and F are replaced by b' and F' on the right side. Likewise, any c 
that can be written as a linear combination of F' rows can be written as a combination of those of 
F, so the two matrices generate the same collection C, and hence have the same solution set S to 
(|5.36p . Since adding to one row a different row can be repeated an arbitrary number of times, (iii) 
holds for an arbitrary (not simply an invertible) multiple of a row. 

The corresponding column operations on a generator matrix are (i) permuting two columns; (ii) 
multiplying a column by an invertible integer; (iii) adding (mod D) to one column an arbitrary 
multiple (mod D) of a different column. Throwing away (or adding) columns of zeros is not an 
allowed operation. When column operations are carried out to produce a new F' from F, the new 
collections C and S' obtained using (|5.35p and (|5.36|) will in general be different, but C" is an 
additive subgroup of the same size (order), \C'\ = |C|, and likewise \S'\ = \S\. The argument is 
straightforward for (i) and (ii), and for (iii) it is an easy exercise to show that if the second column 
of F is added to the first to produce F' , the collection C is mapped into C by the map 

c'l = ci + C2 (mod D); c[ = ci for I > 2 (5.39) 

whose inverse will map C into C when one generates F from F' by subtracting the second column 
from the first. Thus |C| = \C'\. The same strategy shows that = \S\; instead of (|5.39p use 
S2 = S2 — si (mod D), and sj = s; for I ^ 2. 

The row and column operations can be used to transform the generator matrix to a (non unique) 
diagonal form, in the following fashion. If each Fri is zero the problem is trivial. Otherwise use 
row and column permutations so that the smallest positive integer / in the matrix is in the upper 
left corner r = I = I. Suppose / does not divide some element, say F13, in the first row. Then by 
subtracting a suitable multiple of the first column from the third column we obtain a new generator 
F' with < ^"{3 < /, and interchanging the first and third columns we have a generator with a 
smaller, but still positive, element in the upper left corner. Continue in this fashion, considering 
both the first row and the first column, until the upper left element of the transformed generator 
divides every element in both. When this is the case, subtracting multiples of the first column from 
the other columns, and multiples of the first row from the other rows, will yield a matrix with all 
zeros in the first row and first column, apart from the nonzero upper left element at r = 1 = ^, 
completing the first step of diagonalization. 

Next apply the same overall strategy to the sub matrix obtained by ignoring the first row and 
column. Continuing the process of diagonalization and discarding rows that arc all zero (or perhaps 
adding them back in again), one arrives at a diagonal n x n generator matrix 

Frl = flSrU (5.40) 

where some of the // may be zero. The counting problem is now much simplified, because for each 
I ci can be any multiple mod D of //, and s; any solution to //S/ = (mod Z?), independent of 
what happens for a different /. Denoting these two collections by Ci and S'/, the lemma implies that 
\Ci \ ■ \Si\ = D for every ^, and taking the product over / from 1 to n yields (|5.33p . This in turn 
implies that C consists of all possible c that satisfy (|5.32p for all the s € 5. To see this, note that 
the size |C| of C is D"^ /\S\. If we interchange the roles of C and S in the above argument (using a 
generator matrix for S, etc.), we again come to the result (|5.33p . this time interpreting |C| as the 
number of solutions to (|5.32p with S given. Thus since it cannot be made any larger, the original 
additive subgroup C we started with is maximal relative to S. This completes the proof. 
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6 

Location of quantum information 
in additive graph codes 

6.1 Introduction 

Quantum codes in which quantum information is redundantly encoded in a cohection of code carriers 
play an important role in quantum information, in particular in systems for error correction and in 
schemes for quantum communication |Sch95l ISho95 , BBP~'"96[ ICS96| . They are a generalization of 
the classical codes well known and widely used in everyday communication systems |MS77| . While for 
the latter it is fairly obvious where the information is located, the quantum case is more complicated 
for two reasons. First, a quantum Hilbert space with its non-commuting operators is a more complex 
mathematical structure than the strings of bits or other integers used in classical codes. Second, the 
very concept of "information" is not easy to define in the quantum case. However, in certain cases 
one is able to make quite precise statements. Thus in the five qubit code LMPZ96J that encodes 
one qubit of information, none of the encoded information is present in any two qubits taken by 
themselves, whereas all the information can be recovered from any set of three qubits |Gri05j . 

Similar precise statements can be made, as we shall see, in the case of an additive graph code 
on a collection of n qudits which constitute the carriers of the code, provided each qudit has the 
same dimension D, with D some integer greater than one (not necessarily prime). It was shown in 
[LYGGOSj that all additive graph codes are stabilizer codes, and in |Schb[ IGKR] that all stabilizer 
codes are equivalent to graph codes for prime D. A de tailed disc ussion of non-binary quantum error 
correcting codes can be found in jGotal [RaiMl ISWOl) IhTZ+08[ ILYGGOSj . The five qubit code just 
mentioned is an example of a quantum code that is locally equivalent to an additive graph code 
[SW01| . and the information location has an "all or nothing" character. In general the situation 
is more interesting in that some subset of carriers may contain some but not all of the encoded 
information, and what is present can be either "classical" or "quantum." or a mixture of the two. 
Since many of the best codes currently known are additive graph codes, identifying the location 
of information could prove useful when utilizing codes for error correction, or designing new or 
better codes, or codes that correct some types of errors more efficiently than others [IM07j . Our 
formalism can also be applied to study quantum secret sharing schemes employing graph states and 
can even handle a more general setting where there might be subsets that contain partial information 
and hence are neither authorized (contain the whole quantum secret) nor unauthorized (contain no 
information whatsoever about the secret). 

Our approach to the problem of information location is algebraic, based upon the fact that 
generalized Pauli operators on the Hilbert space of the carriers form a group. Subgroups of this 
group can be associated with different types of information, and the information available in some 
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subset of the carriers can also be identified with, or is isomorphic to, an appropriate subgroup, as 
indicated in the isomorphism theorem of Sec. 16.51 In the process of deriving this theorem we go 
through a scries of steps which amount to an encoding procedure that takes the initial quantum 
information and places it in the coding subspace of the carrier Hilbert space. These steps can in 
turn be transformed into a set of quantum gates to produce an explicit circuit that carries out the 
encoding. This result, although somewhat subsidiary to our main aims, is itself not without interest, 
and is an alternative to a previous scheme |SW01j limited to prime D. 

There have been some previous studies of quantum channels using an algebraic approach similar 
to that employed here. Those most closely related to our work are by Beny et al. [BKKOTa, BKK07b] 
(and see Beny |Ben| ) and Blume-Kohout et al. |BKNPV08) . These authors have provided a set of 
very general conditions under which an algebraic structure is preserved by a channel. In App. 16. Dl 
we show that our results fit within the framework of a "correctable algebra" as defined in [BKKOTal 
IBKK07bl IBenj . See also the remarks in Sec. 16.71 

The remainder of this chapter is organized as follows. Some general comments about types 
of quantum information and their connection with certain ideal quantum channels are found in 
Sec. 16.21 Section [6.31 contains definitions of the Pauli group and of some quantum gates used later in 
the chapter. The formalism associated with additive graph codes as well as our encoding strategy is 
in Sec. l6.4| this along with some results on partial traces leads to the fundamental isomorphism result 
in Sec. 16. 5[ which also indicates some of its consequences for the types of information discussed in 
Sec. 16.21 Section r6.6l contains various applications to specific codes, for both qubit and qudit carriers. 
Finally, Sec. 16.71 contains a summary, conclusions, and some open questions. Appendices 16. Al and 
16. Bl contain longer proofs of theorems, App. 16. Cl presents an efficient linear algebra based algorithm 
for working out the results for any additive graph code, and App. [6^ illustrates the connection with 
related work in |BKK07a! and |BKK07b| . 

6.2 Types of information 

Both classical and quantum information theory have to do with statistical correlations between 
properties of two or more systems, or properties of a single system at two or more times. In the 
classical case information is always related to a possible set of physical properties that are distinct 
and mutually exclusive — e.g., the voltage has one of a certain number of values — with one and only 
one of these properties realized in a particular system at a particular time. For quantum systems 
it is useful to distinguish different types or species of information |Gri07j . each corresponding to a 
collection of mutually distinct properties represented by a (projective) decomposition = {Jj} oi 
the identity / on the relevant Hilbert space H: 

I^T.-^^' J,^j] = Jl JjJk=6jkJj. (6.1) 
j 

Any normal operator M has a spectral representation of the form 

M = ^/.,J„ (6.2) 

j 

where the /ij are its eigenvalues, and the decomposition {Jj} is uniquely specified by requiring 
/ij 7^ /ife when j 7^ k. This means one can sensibly speak about the type of information JT{M) 
associated with a normal operator M. When M is Hermitian this is the kind of information obtained 
by measuring Al. 

This terminology allows one to discuss the transmission of information through a quantum chan- 
nel in the following way. Let £ be the completely positive, trace preserving superoperator that maps 
the space of operators CCH) of the channel input onto the corresponding operator space C{H') of 



6. LOCATION OF QUANTUM INFORMATION IN ADDITIVE GRAPH CODES 



66 



the channel output W (which may have a different dimension from T-L). Provided 

£{Jj)£{Jk) = Oior J ^k, (6.3) 

for all the operators {Jj} associated with a decomposition J of the H identity, we shall say the chan- 
nel is ideal or noiseless for the species of information, or, equivalently, the J7 type of information 
is perfectly present in the channel output Ji' . Formally, each physical property Jj at the input 
corresponds in a one-to-one fashion to a unique property, the support oi£{Jj) (or the corresponding 
projector) at the output. Thus we have a quantum version of a noiseless classical channel, a device 
for transmitting symbols, in this case the label j on Jj, from the input to the output by associating 
distinct symbols with distinct physical properties — possibly a different collection of properties at the 
output than at the input. 

The opposite extreme from a noiseless channel is one in which £{Jj) is independent of j up to a 
multiplicative constant. In this case no information of type JJ is available at the channel output: the 
channel is blocked, or completely noisy; equivalently, the J7 species of information is absent from the 
channel output. Hereafter we shall always use "absent" in the strong sense of "completely absent", 
and the term present, or partially present for situations in which some type of information is not 
(completely) absent but is also not perfectly present: i.e., the channel is noisy but not completely 
blocked for this type of information. 

In some cases all the projectors in {Jj} will be of rank 1, onto pure states, but in other cases 
some or all of them may be of higher rank, in which case one may have a refinement C = {Li\ 
of {Jj} such that each projector Jj is a sum of one or more projectors from the C decomposition. 
It is then clear that if the C information is absent /perfectly present from/in the channel output 
the same is true of the J information, but the converse need not hold. Thus it may be that the 
coarse grained information is perfectly present, but no additional information is available about 
the refinement. A particularly simple situation, which we will encounter later, is one in which the 
output v.' is itself a tensor product, say 'H'i^'H'2, J a, decomposition of H'l, J = { Jj- ® /} and K. a 
decomposition of l-i^, JC = {I ® Kk}- It can then be the case that the information associated with 
the decomposition is perfectly present and that associated with the /C decomposition is (perfectly) 
absent from the channel output. 

Suppose = {Jj} and /C = {Kk} are two types of quantum information defined on the same 
Hilbert space. The species and K. are compatible if all the projectors in commute with all 
the projectors in /C, in which case the distinct nonzero projectors in the collection {JjKk} provide 
a common refinement of the type discussed above. Otherwise, if some projectors in one collection 
do not commute with certain projectors in the other, and K. are incompatible and cannot be 
combined with each other. This is an example of the single framework rule of consistent quantum 
reasoning, [Gri96| or Ch. 16 of jGri02| . The same channel may be ideal for some and blocked 
for some /C, or noisy for both but with different amounts of noise. From a quantum perspective, 
classical information theory is only concerned with a single type of (quantum) information, or several 
compatible types which possess a common refinement, whereas the task of quantum, in contrast to 
classical, information theory is to analyze situations where multiple incompatible types need to be 
considered. 

The term "classical information" when used in a quantum context can be ambiguous or mis- 
leading. Generally it is used when only a single type of information, corresponding to a single 
decomposition of the identity, suffices to describe what emerges from a channel, and other incom- 
patible types can therefore be ignored. Even in such cases it is helpful to indicate explicitly which 
decomposition of the identity is involved if that is not obvious from the context. The contrasting 
term "quantum information" can then refer to situations where two or more types of information 
corresponding to incompatible decompositions are involved, and again it is helpful to be explicit 
about what one has in mind if there is any danger of ambiguity. 

An ideal quantum channel is one in which there is an isometry V from JL to "H' such that 

£{A) = VAV^ (6.4) 
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for every operator A on H. In this case the superoperator £ preserves not only sums but also 
operator products: 

£iAB)^£{A)£{B). (6.5) 

Conversely, if (|6.5p holds for any pair of operators, one can show that the quantum channel is ideal 
[BKK07al IBKK07b] . i.e. £ has the form (|6.4p . As the isometry maps orthogonal projectors to 
orthogonal projectors, (|6.3p will be satisfied for every species of information, and we shall say that 
all information is perfectly present at the channel output. The converse, that a channel which is 
ideal for all species, or even for an appropriately chosen pair of incompatible species is an ideal 
quantum channel, is also correct; see jGri05| [Gri07| . 

The preservation of operator products, (j6.5p . can be a very useful tool in checking for the presence 
or absence of various types of information in the channel output, as we shall see in Sec. 16.51 When 
(|6.5[) holds for arbitrary A and B belonging to a particular decomposition of the identity, this suffices 
to show that the channel is ideal for this species. However, note that this sufficient condition is not 
necessary, since (|6.3p could hold without the £{Aj) being projectors, in which case £{A'j) is not 
mapped to £{Aj)^. 

We use the term ideal classical channel for a type of information — {Jj} to refer to a situation 
where ()6.3p is satisfied and, in addition, 

£{JjAJk) = for j ^ k, (6.6) 

where A is any operator on the input Hilbert space H. That is, not only is type perfectly 
transmitted, but all other types are "truncated" relative to this type, in the notation of |Gri96| . 



6.3 Preliminary remarks and definitions 
6.3.1 Generalized Pauli operators on n qudits 

We generalize Pauli operators to higher dimensional systems of arbitrary dimension D in the follow- 
ing way. The X and Z operators acting on a single qudit are defined as 

D-l D-l 

Z=J2c.^j){j\, + (6.7) 

j=0 j=0 

and satisfy 

X" = Z^ = I, XZ = ujZX, uj = c^'^'/-^, (6.8) 

where the addition of integers is modulo D, as will be assumed from now on. For a collection of n 
qudits we use subscripts to identify the corresponding Pauli operators: thus Zi and Xi operate on 
the space of qudit i. The Hilbert space of a single qudit is denoted by H, and the Hilbert space of 
n qudits by T-Ln ■, respectively. Operators of the form 

^A^x^z ._ ^^x^'Z^' (g) X^'Z^^ ® • • • ® X^'^Z^" (6.9) 

will be referred to as Pauli products, where A is an integer in Z^i and x and z are ?7-tuplcs in ZJ,, the 
additive group of n-tuple integers mod D. For a fixed n the collection of all possible Pauli products 
(|6.9p form a group under operator multiplication, the Pauli group Vn- If p is a Pauli product, then 

= / is the identity operator on Tin, and hence the order of any element of Vn is either D or else 
an integer that divides D. While Vn is not abelian, it has the property that two elements commute 
up to a phase: piP2 = w'*"P2Pi, with A12 an integer in Zd that depends on pi and p2- 

The collection of Pauli products with A = 0, i.e. a pre-factor of 1, is denoted by Q„ and forms an 
orthonormal basis of C{'Hn), the Hilbert space of linear operators on with respect to the inner 
product 

7^Tr[9jg2] = '^qi,,^, Vqi,(j2eQ„. (6.10) 
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Note that Q„ is a projective group or group up to phases. There is a bijcctive map between Q„ 
and the quotient group Vn/i^^^I} for A € Zd where {uj'^I}, the center of Vn, consists of phases 
multiplying the identity operator on n qudits. 

6.3.2 Generalization of qubit quantum gates to higher dimensions 

In this subsection we define some one and two qudit gates generalizing various qubit gates. The 
qudit generalization of the Hadamard gate is the Fourier gate 



D-l 

For an invertible integer q € (i.e. integer for which there exists q S Z/j such that qq = 1 mod D), 
we define a multiplicative gate 

D-l 

:= 1^') ' (6.12) 

where qj means multiplication mod D. The requirement that q be invertible ensures that Sg is 
unitary; for a qubit Sg is just the identity. 

For two distinct qudits a and b we define the CNOT gate as 

D-l D-l 

CNOT,, := \J) 01 ® ^6 = E 1^') OL ® + , (6.13) 

the obvious generalization of the qubit Controlled-NOT, where a labels the control qudit and b labels 
the target qudit. Next the SWAP gate is defined as 

D-l 

SWAP„b:= E \k){j\a^\j){k\b- (6-14) 

],k=0 

It is easy to check that SWAP gate is hermitian and does indeed swap qudits a and b. Unlike the 
qubit case, the qudit SWAP gate is not a product of three CNOT gates, but can be expressed in 
terms of CNOT gates and Fourier gates as 

SWAP,, = CN0T,,(CN0Tb,)tCN0T,,(F2 ® /,), (6.15) 

with 

(CNOTba)t = (CNOTbJ^-i = {la ® F2)CNOT6a(/a ® F^). (6.16) 
Finally we define the generalized Controlled-phase or CP gate as 

CPa, = E 1^') 01a ® = Y ^'"IJ) OL ® \k) {k\, . (6.17) 

j=0 j,k=0 

The CP and CNOT gates are related by a local Fourier gate, similar to the qubit case 

CNOTab = (/a ® Fb)CPab(/a ® Fb)\ (6.18) 

since F maps Z into X under conjugation (see Table I5TT|) . 

The gates F, Sg, SWAP, CNOT and CP are unitary operators that map Pauli operators to Pauli 
operators under conjugation, as can be seen from Tables 16.11 and 16.21 They are elements of the so 
called Clifford group on n qudits |Gota[ IHDMOS] , the group of n-qudit unitary operators that leaves 
Vn invariant under conjugation, i.e. if O is a Clifford operator, then Vp G Vn, OpO^ G Vn- From 
Tables [01 and one can easily deduce the result of conjugation by F, Sg, SWAP, CNOT and CP 
on any Pauli product. 
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Pauli operator 




F 


Z 


Zi 


X 


X 


Xt 





Tabic 6.1: The conjugation of Pauli operators by one-qudit gates F and Sg {q is the multiplicative 
inverse of q mod D) . 



Pauli product 


CNOTaf, 


SWAP.b 


GVab 




Za®Zb 


Za(Slb 


Ia®Zb 


Za<E)Ib 


Za®h 




Za®Ib 


Ia®Xb 


Ia®Xb 


Xa®Ib 


z^-' ® Xb 


Xa®h 


Xa®X^-' 


Ia®Xb 


Xa ® Z^-' 



Table 6.2: The conjugation of Pauli products on qudits a and b by two-qudit gates CNOT, SWAP 
and CP. For the CNOT gate, the first qudit a is the control and the second qudit b the target. 



6.4 Graph states, graph codes and related operator groups 
6.4.1 Graph states and graph codes 

Let G = {V,E) be a graph with n vertices V, each corresponding to a qudit, and a collection E 
of undirected edges connecting pairs of distinct vertices (no self loops). Two qudits can be joined 
by multiple edges, as long as the multiplicity does not exceed D — 1. The graph G is completely 
specified by the adjacency matrix F, where the matrix element Tab represents the number of edges 
that connect vertex a with vertex b. The graph state 

|G) =[/|Go) =(7(1+)^") (6.19) 

is obtained by applying the unitary (Clifford) operator 

U= H {GPabf-', (6.20) 

(a,b)eE 

where each pair (a, b) of vertices occurs only once in the product, to the trivial graph state 

\Go) ■■= 1+)^", (6.21) 

with 

1+) ^= ^ E l^')- (6-22) 

^ i=o 

Define to be the stabilizer of |G) , the subgroup of operators from Vn that leave |G) unchanged. 
The stabilizer Sq' of the trivial graph state |Go) is simply the set of all X-type Pauli products with 
no additional phases, 

= {X^ :x = (a;i,a;2,...,a:„)}, (6.23) 

where Xj are arbitrary integers between and D — 1. Since |G) is related to |Go) through a Clifford 
operator (see (|6.19[) and (|6.20[) 1. it follows at once that the stabilizer S'^ of |G) is related to the 
stabilizer Sq of the trivial graph through the Clifford conjugation 



5*^ = USl^U\ 



(6.24) 
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with U defined in ([Qg]) . 

A graph code C ean be defined as the if-dimensional subspacc He of Hn spanned by a collection 
of K mutually orthogonal codewords 

|c,) = ZC^|G), j^l,2,...,K (6.25) 

where 

Cj = (cji,Cj2, ■ • • ,Cj«) (6.26) 

is for each j an n-tuple in Z^). The Cjk notation suggests a matrix c with K rows and n columns, of 
integers between and D — 1, and this is a very helpful perspective. In this chapter we are concerned 
with additive graph codes, meaning that the rows of this matrix form a group under component- wise 
addition mod D, isomorphic to the abelian coding group C, of order |C| = if, of the operators Z^^ 
under multiphcation. We use (C,|G)) to denote the corresponding graph code. For more details 
about graph states and graph codes for arbitrary £>, see |LYGG08] . 

Note that the codeword (0,0, ... ,0) is just the graph state |G), and in the case of the trivial 
graph |Go) this is the tensor product of |+) states, (|6.2ip . not the tensor product of |0) states which 
the n-tuple notation (0,0, ... ,0) might suggest. Overlooking this difference can lead to confusion 
through interchanging the role of X and Z operators, which is the reason for pointing it out here. 



6.4.2 The encoding problem 

A coding group C can be used to create an additive code starting with any n qudit graph state, 
including the trivial graph |Go), because the entangling unitary U commutes with Z^ for any z; 
thus 

|c,) = Z'^^C/IGo) = C/Z'^^IGo) = C/|c^") (6.27) 

where the |c^) span the code (C, \Gq)). But in addition the coding group C is isomorphic, as explained 
below to a trivial code Co, 

Co = (Z^"^Z2"^...,Z^"^) (6.28) 

which is generated by, i.e., includes all products of, the operators inside the angular brackets ( ). 
Here k is an integer less than or equal to n, and each ruj is 1 or a larger integer that divides D. The 
simplest situation is the one in which each of the mj is equal to 1, in which case Co is nothing but 
the group, of order D^ , of products of Z operators to any power less than D on the first k qudits. 
One can think of these qudits as comprising the input system through which information enters the 
code, while the remaining n — k qudits, each initially in a 1+) state, form the ancillary system for 
the encoding operation. 

If. however, one of the nij is greater than 1, the corresponding generator Z™^ is of order 

d-i = D/in-j, (6.29) 

and represents a qudit of dimensionality dj rather than D. Thus for example, if Z) = 6 and mi = 2, 
applying Z^ and its powers to |+) will produce three orthogonal states corresponding to a qutrit, 
di = 3. (Identifying operators Z and X on these three states which satisfy (|6.8p with D = 3 is not 
altogether trivial, and is worked out in Sec . 16 . 4 below . ] In general one can think of the group Co in 
(|6.28p as associated with a collection of k qudits, the j'th qudit having dimension dj, and therefore 
the collection as a whole a dimension oi K = did2 ■ ■ ■ dk, equal to that of the graph code. If one 
thinks of the information to be encoded as initially present in these k qudits, the encoding problem 
is how to map them in an appropriate way into the coding subspacc Ji of the n (Z3-dimensional) 
carriers. 

We address this by first considering the connection between C and Co in a simple example with 
71 = 3, I? = 6, and 

C^{ztzlzlzlzl), (6.30) 
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a coding group of order 6. The two generators in (|6.30p correspond, in the notation introduced in 
(|6.26p . to the rows of the 2x3 matrix 



4 3 3 
3 3 



.31) 



By adding rows or multiplying them by constants mod D one can create 4 additional rows which 
together with those in (|6.3ip constitute the 6 x 3 c matrix. 

Through a sequence of elementary operations mod D — a) interchanging of rows/columns, b) 
multiplication of a row/column by an invertible integer, c) addition of any multiple of a row/column 
to a different row/column — a matrix such as f can be converted to the Smith normal form [New721 
[5to96] 

s = vf-w, (6.32) 

where v and w are invertible (in the mod D sense) square matrices, and s is a diagonal rectangular 
matrix, as in (|6.33[) . It is proved in |Sto96| that a. K x n matrix can be reduced to the Smith form 
in only 0{K^^^n) operations from Z^, where 9 is the exponent for matrix multiplication over the 
ring Zd; i.e. two m x m matrices over Zd can be multiplied in 0{m^) operations from Zd. Using 
standard matrix multiplication 9 = 3, but better algorithms |CW87] allow for 9 = 2.38. 
For the example above, the sequence 



4 3 3\ /4 0\ /4 0\ /2 
3 3/ ^ lo 3 3/^ 1 3 0/ ^ lo 3 



(6.33) 



proceeds by adding the second row of f to the first (mod 6), then the second column to the third 
column, and finally multiplying the first row by 5 (which is invertible mod 6). The final step is 
needed so that the diagonal elements divide D: mi = 2, m2 = 3, so that c?i = 3 and d2 = 2. Thus 
we arrive at the trivial coding group 

Co = (z2,Z|), (6.34) 

isomorphic to C in (|6.30|) . 

Since the procedure for reducing a matrix to Smith normal form is quite general, the procedure 
illustrated in this example can be applied to any coding group C, as defined following (|6.26p . to find 
a corresponding trivial coding group Cq. The row operations change the collection of generators but 
not the coding group that they generate; i.e., the final collection of K rows is the same. The column 
operations, on the other hand, produce a different, but isomorphic, coding group, and one can think 
of these as realized by a unitary operator W which is a product of various SWAP, CNOT and Sq 
gates, so that 

C = WCoW\ (6.35) 

that is, conjugation by W maps each operator in Co to its counterpart in C. In our example, 
W = CNOT32 is the only column operation, the second arrow in (|6.33p . and represents the first 
step in the encoding circuit for this example, Fig. l6.1f bV It is left as an exercise to check that this 
relates the generators in (|6.30p and (|6.34p through (j6.35p . Table indicates how different matrix 
column operations are related to the corresponding gates in the encoding circuit. 
The overall encoding operation 

|c,) = UW\c^) (6.36) 

starting with the trivial code on the trivial graph (Cq, |Go)) and ending with the desired code (C, \G)) 
is shown for our example in Fig. I6.ir b) for the case of a graph indicated in (a) in this figure. It is 
important to notice that both W and J7, and therefore their product, are Clifford operators, unitaries 
that under conjugacy map Pauli products to Pauli products. This follows from the fact that the 
gates in Table [6?3l are Clifford gates, and will allow us in what follows to extend arguments that are 
relatively straightforward for trivial codes on trivial graphs to more general additive graph codes. 
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(a) (b) 



Figure 6.1: (a) The graph state used in the example; (b) The encoding circuit: the input states 
^Ci"i.i^C2m2| _^ ^^^^ correspond to the trivial code Cq arc mapped by W to C, then U entangles 
the qudits. Here mi = 2, 7712 = 3 and Q are integers such that ^ Q ^ dj — 1, with di = 3, ^2 = 2. 



Matrix operation in Zd 


Clifford conjugation 


Interchange of columns a and b 


SWAPab 


Multiplication of column a 
by invertible integer q 


Sq on qudit a 


Addition of m times column b to 
column a 


{CNOTabY" 



Table 6.3: The correspondence between matrix column operations in Z^j and conjugation by Clifford 
gates. For the CNOT gate, the first qudit a is the control and the second qudit b the target. 



6.4.3 The information group 

In this section we define the information group that plays a central role in the isomorphism theorem 
in Sec. 16.51 below. The basic strategy is most easily understood in terms of Co = (Co, |Go)), the 
trivial code on the trivial graph. However, because the overall encoding map UW in (j6.36p is a 
Clifford operation mapping Pauli products to Pauli products, various results that apply to Co can 
be immediately translated to the general graph code C = (C, |C)) we are interested in, and for this 
reason most of the formulas valid for both are written in the form valid for C even if the derivations 
are based on Co. 

The pointwise stabilizei0 of Cq , the subgroup of operators from Vn that leave every codeword 
|c^) unchanged, is given by 

So = {X^ : X = {r]idi,r]2d2, . . . ,r]kdk,Xk+i, ■ . . ,x„)} , (6.37) 

where the dj are defined in ()6.29|) . rjj is any integer between and rrij — 1, and the Xj for j > k 
are arbitrary integers between and D — I. That this is correct can be seen as follows. First, Pauli 
products belonging to Sq cannot contain Zj operators, for such operators map each codeword onto 
an orthogonal state. On the other hand, every X^^ leaves |Go), (|6.2ip . unchanged, so it belongs 
to Sq if and only if it commutes with Z™'^ , which means XjUij must be a multiple oi D, or Xj 
a multiple of dj, see (|6.29p . Thus elements of Sq commute with elements of Cq, (|6.28p . Since its 
operators cannot alter the phases of the codewords, no additional factors of lu^ are allowed, and 

^Also called the "fixer" or "fixator". It is important to distinguish this subgroup from the group theoretical notion 
of the stabilizer of the coding space in the sense of the subgroup of P„ that maps the coding space onto itself without 
necessarily leaving the individual vectors fixed. As we shall not employ the latter, it should cause no confusion if we 
hereafter follow the usual convention in quantum codes and omit "pointwise," even though retaining it would add 
some precision. 
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thus Sq is given by (|6.37p . The stabihzcr of the (nontrivial) code C is then the isomorphic group S 
obtained using the unitary UW of (|6.36p : 

S = {UW)So{UW)'' = {{UW)s{UW)^ : s S So}, (6.38) 

a coUection of Pauh products because the unitary UW, as remarked earher, is a Chfford unitary. 
The order of Sq, and thus of S, is given by 

h 

\S\ = i?"-^- n m, = — r = Twr = — • (6-39) 

Next define the subgroup W of Vn 

W=(5^,C) (6.40) 

generated by operators belonging to the stabihzer S'^ of the graph state or to the coding group C, 
and denote it by Wo = {S^,Co) in the case of the trivial code. The elements of Sq commute with 
those of (both are abelian and the former is a subgroup of the latter), and also with those of Co, 
as noted above. As group properties are preserved under the UW map, as in (j6.38p . we conclude 
that all elements in S commute with those in W, even though W is not (in general) abelian, and 
hence 5 is a normal subgroup of W. Now define the abstract information group as the quotient 
group 

g = W/S ^ {S^,C)/S (6.41) 

consisting of cosets of S, written as gS or Sg for g in W. Note that because any element g of W is 
a Pauli product, = I is the identity, and the order of g is either D or an integer that divides D. 
Consequently the order of any element of G is also D or an integer that divides D. 
To understand the significance of Q consider a trivial code on a single qudit, with 

Co = (zr), S^ = {X,), So = {xt') (6.42) 

The elements of Qq can be worked out using its identity / and the generators X and Z: 

I^So = {Ii,xt\X^,''\...} 
X = XiSo = {Xi,X^^+\X^'^^+\ . . .} 

Z = 1^0 = {Zf S ^XfS . . .}. (6.43) 

It is evident that the cosets X, X^ ~ XfSo and so forth up to X'^'^~^ are distinct, whereas X"^^ = 
I = So- The same is true for powers of Z. Furthermore, 

XZ = XiZ"''So = w^^Z^^Xi^o = LoZX, (6.44) 

with w = a;"'i = (j2;ri/di^ Thus Qo is generated by operators X and Z that satisfy (|6.8p with D 
replaced by di, which is to say the corresponding group is what one would expect for a qudit of 
dimension di . The same argument extends easily to the trivial code on k carriers produced by Co , see 
(|6.28|) : Go is isomorphic to the group of Pauli products on a set of qudits of dimension di,d2, ■ ■ ■ ,dk. 
The same structure is inherited by the abstract information group G for the code C = {C,\G)) 
obtained by applying the UW map as in (|6.38p . 

The abstract information group G is isomorphic to the information group G of information 
operators acting on the coding space He a-nd defined in the following way. Its identity is the 
operator 

P=|5ril](5) = 151-1^3, (6.45) 

sG5 
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where S(^) denotes the sum of the operators that make up a collection A. In fact, P is just the 
projector onto "He, as can be seen as follows. Since 5 is a group, P^ = P; and since a group contains 
the inverse of every element, and s € 5 is unitary (a Pauli product), P^ = P. These two conditions 
mean that P is a projector onto some subspace of Hn- Since S is the (pointwise) stabilizer of the 
coding space each s in 5 maps a codeword onto itself, and thus P maps each codeword to itself. 
Consequently, all the codewords lie in the space onto which P projects. Finally, the rank of P is 

Tr[P] = D"/\S\ = |C| = K (6.46) 

(see (|6.39p ). since the trace of every s in S is zero except for the identity with trace (Note that 
while Vn contains the identity multiplied by various phases, only the identity operator occurs in S.) 
Therefore P projects onto "He, and is given by the formula 

K 

^ = El'^j)(^^l- (6.47) 

i=i 

The other information operators making up the information group G = {g} are formed in a 
similar way from the different cosets making up W/S: 

g = \S\-'^igS) = gP = PgP = PgP. (6.48) 

That is, for each coset form the corresponding sum of operators and divide by the order of the 
stabilizer S. The second and third equalities in (|6.48p reflect the fact that the product of the cosets 
S and gS in either order is gS, which is to say P forms the group identity of Q. They also tell us 
that the operators that make up G act only on the coding space, mapping He onto itself, and give 
zero when applied to any element of Jin in the orthogonal complement of He- Because 5 is a normal 
subgroup of W, products of operators of the form ()6.48|) mirror the products of the corresponding 
cosets, so the map from the abstract G to the group C/ is a liomomorphism. That it is actually an 
isomorphism is a consequence of the following, proved in App. 16. Al 

Lemma 6.1. Let TZ be a linearly independent collection of Pauli product operators that form a 
subgroup of Vn, o,nd for a Pauli product p let pTZ ~ {pr : r G TZ}. Then 
i) The operators in pTZ are linearly independent. 

a) If p and q are two Pauli products, one or the other of the following two mutually exclusive 
possibilities obtains: 
a) 

pTZ = e">'qTZ (6.49) 

in the sense that each operator in pTZ is equal to e"^ times an operator in qTZ 
P) The union pTZU qTZ is a collection of 2\TZ\ linearly independent operators. 

Since the collection of Pauli products Q„ with fixed phase forms a basis of C{Tin), a collection 
of Pauli products can be linearly dependent if and only if it contains both an operator and that 
operator multiplied by some phase. As the (pointwise) stabilizer S leaves each codeword unchanged, 
the corresponding operators are linearly independent, and the lemma tells us that distinct cosets 
^iS ^ hS give rise to distinct operators g ^ h. Either gS = e^'^hS, in which case g = e"^/i ^ h 
(since if e"^ = 1 the cosets are identical) . Or else the gS operators are linearly independent of the 
hS operators, and therefore g and h are linearly independent. Consequently the homomorphic map 
from to is a bijection, and the two groups are isomorphic. 

The single qudit example considered in (j6.42p provides an example of how G and G are related. 
In this case the projector 

Po = (l/mi)(/i+Xfi + ---) (6.50) 
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projects onto the subspace spanned by |+), Z™^ |+), Z^™^ |+), While each of the operators that 

make up a coset such as X in (|6.43|) is unitary, their sum, an operator times Pq, is no longer unitary, 
though when properly normalized acts as a unitary on the subspace onto which Pq projects. That 
the different sums of operators making up the different cosets are distinct is in this case evident from 
inspection without the need to invoke Lemma l6.f I 

Let us summarize the main results of this subsection. For an additive graph code C we have 
defined the information group Q of operators acting on the coding subspace "He, whose group 
identity is the projector P onto He- It is isomorphic to the group of Pauli products acting on a 
tensor product of qudits of dimensions di, d2, . . . , dk, which can be thought of as the input to the 
code, see Sec. 16.4.21 Each element g of C/ is of the form PgP, so as an operator on Hn it commutes 
with P and yields zero when applied to any vector in the orthogonal complement of He- The 
dimension of He is K = did2 ■■• dk, the size of the code, and hence the elements of G span the space 
of linear operators /^("Hc) onHc- 

6.5 Subsets of carriers and the isomorphism theorem 
6.5.1 Subsets of carriers 

Before stating the isomorphism theorem, which is the principal technical result of this chapter, let 
us review some facts established in Sec. 16.41 The additive graph code (C, \G)) we are interested 
in can be thought of as arising from an encoding isometry that carries the channel input onto a 
subspace Jic of the rt-qudit carrier space "Hn, as in Fig 16.11 This isometry, as explained in Sec. 16.21 
in connection with (|6.4p . constitutes a perfect quantum channel, and thus all the information of 
interest can be said to be located in the Jic subspace, where it is represented by the information 
group G, a multiplicative group of operators for which the projector P on He is the group identity, 
and which as a group is isomorphic to the abstract information group G defined in (j6.4ip . 

We are interested in what kinds of information are available in some subset B of the carriers, 
where B denotes the complementary set. For this purpose it is natural to consider the partial traces 
over B, i.e., the traces down to the Hilbert space Hb, of the form 

gB ^ N-^Trsig], (6.51) 

where g is an element of the information group G, and the positive constant N is defined in (j6.58p 
below. In those cases in which gs = the J^{g) information has disappeared and is not available in 
the subset B, so we shall be interested in those g for which the partial trace does not vanish, that 
is to say in the elements of the subset information group 

G"" = {geG:Ti-s[g]^0}. (6.52) 

We show below that G^ is a subgroup of G, thus justifying its name, and that it is isomorphic to 
the group Gb of nonzero operators of the form gs defined in (|6.51|) . To actually determine which 
elements of G belong to G^ one needs to take partial traces of the g G G to see which of them do not 
trace down to zero. In App. 16. CI we present an efficient linear algebra algorithm based on solving 
systems of linear equations mod D that can find G^ in 0{K^n^) operations from "Lo where 9 is 
defined in Sec. 

If an operator A on the full Hilbert space Hn of the n carriers can be written as a tensor product 
of an operator on Hb times the identity operator Iq on Hq we shall say that A is based in B. Let 
B be the collection of all operators on Hn that are based in B. Obviously, B is closed under sums, 
products, and scalar multiplication. In addition the partial trace Tr^[A] of an operator A in B is 
"essentially the same" operator, apart from normalization in the sense that 

A^ D-\^\ ■Ttb[A](E)Ib. (6.53) 
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li A ^ B is a. Pauli product, then its partial trace over B vanishes, since Tr[X] and Tr[Z] and their 
powers (when not equal to /) arc zero. Consequently the partial trace over B of T,{gS) in (|6.48p is 
the same as the partial trace oiY,[{gS)r\B], which suggests that it is useful to consider the properties 
of collections of Pauli operators of the form (gS) n B with g an clement of W. The following result, 
proved in App. 16. Al turns out to be useful. 

Lemma 6.2. Let g, h be two arbitrary elements of W, and B the collection of operators with base 
in B. 

i) The set {gS) (1 B is empty if and only if (g^^S) D B is empty. 

ii) Every nonempty set of the form {gS) H B contains precisely 

M ^\SC]B\>1 (6.54) 

elements. 

Hi) Two nonempty sets [gS) n B and {hS) n B are either identical, which means gS = hS and 
ll[{gS) n S] = Yi[[hS) f^B], or else they have no elements in common and the operators Y,[[gS) H B\ 
and Ti[(hS) n S] are distinct. 

iv) If both (gS) D B and (hS) D B are nonempty, their product as sets, including multiplicity, is 
given by 

[{gS) n B] ■ [{hS) nB]^ M[{ghS) n B]. (6.55) 

By (|6.55p we mean the following. The product (on the left) of any operator from the collection 
(gS) n B with another operator from the collection {hS) D B belongs to the collection (ghS) D B (on 
the right), and every operator in (ghS) D B can be written as such a product in precisely M different 
ways. 

We are now in a position to state and prove our central result: 
6.5.2 Isomorphism theorem 

Theorem 6.3 (Isomorphism). Let C be an additive graph code with information group Q, P the 
projector onto the coding space He B be some subset of the carrier qudits. Then the collection 
of members of Q with nonzero partial trace down to B, (|6.52p . is a subgroup of the information 
group Q, and the mapping g ^ gs in (|6.51|) carries to an isomorphic group Qb of nonzero 
operators onHs- Furthermore, 

i) If g and h are any two elements of , then 

Tve[gh] - TrsL?] Tr^i/il/iV or {gh)B = gshB (6.56) 

ii) If g ^ h are distinct elements of , gs are distinct elements of Qb- 
Hi) The identity element 

Pb :=TrB[P]/7V, (6.57) 

of Qb is a projector onto a subspace ofHB (possibly the whole space) with rank equal to Tr[P]/iV = 
K/N. 

The normalization constant N is given as 

N ■.= \SC^B\■D\^\/\S\ (6.58) 

where B are the operators based in B. 

Proof. The proof is a consequence of Lemma 16.21 and the following observations. The trace Tr^[g] 
in ()6.51|) is, apart from a constant, the trace of Yi[{gS) n B], and is zero if [gS] n S is empty. 
If the collection [gS] n K is not empty, then by Lemma 16.11 it consists of a collection of linearly 
independent operators, and the trace of its sum cannot vanish. Thus there is a one-to-one, see part 
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(iii) of Lemma [n?21 correspondence between nonempty sets of the form {gS) n B and the elements g 
in . Then (i) and (iv) of Lemma 16.21 imply both that is a group, and also that the map from 
to Gb is a homomorphism, whereas (ii) shows that this is actually an isomorphism: gs ~ hs 
is only possible when gS ~ hS. That in (|6.58p is the correct normalization follows from (|6.54p . 
(16351) . and ([OSl) . □ 

A significant consequence of Theorem 16.31 is the following result on the presence and absence of 
information in the subset _B, using the terminology of Sec. 16.21 

Theorem 6.4. Let C be an additive graph code on n carrier qudits, with information group Q . Let 
B he a subset of the carrier qudits, the corresponding subset information group, and J{g) the 
type of information corresponding to g (as defined in Sec. \6.^ ). Then 

i) The J{g) type of information is perfectly present in B if and only if g E . 

ii) The J{g) type of information is absent from B if and only if g^ ^ for all integers k between 
1 and D ^ 1. 

iii) All information is perfectly present in B if and only if = Q. 

iv) All information is absent from B if and only if consists entirely of scalar multiples of the 
identity element P of Q . 

The proof of the theorem can be found in App. 16. Bl Statement (iii) is useful because the check 
of whether there is a perfect quantum channel from the input to B involves a finite group Q\ one 
does not have to consider all normal operators of the form (|6.2p . Statement (ii) deserves further 
comment. If D is prime then the order of any element of the Pauli group (apart from the identity) is 
D, see the remark following (|6.9p . The same is true of elements of the quotient group Q, (|6.4ip . and 
thus of members g of the isomorphic group Q. Consequently, for any k in the interval \ < k < D, 
there is some m such that 1 = km mod £>, which means g = {g^Y'^ ■ And since is a group, 
~k g gB jj^pijgg g g gB 'pjjyg -^vhen D is prime, g ^ is equivalent to g^ ^ for all integers 
k between 1 and D — 1, and the latter can be replaced by the former in statement (ii). However, 
when D is composite it is quite possible to have Tr^ [g] = but Tr^ [g'^ ] ^ for some k' larger than 
1 and less than D; see the example below. In this situation we can still say that J^{g'' ) is perfectly 
present, but it is not true that J{g) is absent. One can regard the type JT{g) as a refinement of 
J'{g'' ), and as explained in Sec. 16. 2[ although the coarse-grained J^{g'' ) information is perfectly 
present in B, the additional information associated with the refinement is not. 

As an example, suppose g has a spectral decomposition 

g = Jo + iJi-J2-U3, (6.59) 
with the Jj orthogonal projectors such that 

Tr^ [Jo] = Trs [J2] ^ Tr^ [ Ji] = Tr^ [ J3] • (6.60) 

Then Tr^[g] = 0, whereas 

5' = (Jo + J2)-(Ji + J3), (6.61) 
and thus Trg[g^] 7^ 0. Thus cp' is an element of , whereas g is not, and so the coarse grained 
Cfiff') information corresponding to the decomposition on the right side of (|6.61|) is present in B, 
while the further refinement corresponding to the right side of (|6.59p is not. Precisely this structure 
is produced by a graph code on two carriers of dimension D = A, with graph state \G) = | + +), 
coding group C = {ZxZ^j, information group Q = {XiP, Z1Z2P) , coding space projector 

P={I + X,Xl + XfXi + XiX2)/A, (6.62) 

and 

g = XiP = |00) (00| + i|12) (12| - |20) (20| - i|32) (32| , (6.63) 
where |j) = Z^\+) are the eigenvectors of the X operator. 
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6.5.3 Information flow 

At this point let us summarize how we think about information "flowing" from the input via the 
encoding operation into a subset B of the code carriers. At the input the information is represented 
by the quotient group Gq ^ Wq/Sq, see (|6.4ip . or more concretely by the isomorphic group Go of 
operators generated by the cosets. as in (j6.48p . The encoding operation UW, see (j6.36p and (|6.38p . 
maps Go to the analogous G = W/S associated with the code C, and likewise Go to the group of 
operators G acting on the coding space He- Tracing away the complement B of B maps some of 
the g operators of G to zero, and the remainder form the subset information group G^ ■ Applying 
the inverse UW map to G^ gives Go a subgroup of Go that tells us what types of information at 
the input (i.e. before the encoding) are available in the subset of carriers B. This is illustrated by 
various examples in the next section. 

6.6 Examples 

6.6.1 General principles 

In this section we apply the principles developed earlier in the chapter to some simple [[n,fc, (J]]^) 
additive graph codes, where n is the number of qudit carriers, each of dimension D, the dimension 
of the coding space Tic is K = D^ , and 8 is the distance of the code; see Chapter 10 of |NCOO| for 
a definition of 5. We shall be interested in the subset information group G^ ^ (|6.52p . that represents 
the information about the input that is present in the subset B of carriers. Rather than discussing 
G^ or its traced down counterpart Gb, it will often be simpler to use Go i the subset information 
group referred back to the channel input, see Sec. 16.5.31 above, and in this case we add an initial 
subscript to operators: Xqi means the X operator on the first qudit of the input. Since all three 
groups are isomorphic to one another, the choice of which to use in any discussion is a matter of 
convenience. (In the examples below for the sake of brevity we sometimes omit a term e"^/ from 
the list of generators of Go ■) 

Before going further it is helpful to list some general principles of quantum information that 
apply to all codes, and which can simplify the analysis of particular examples, or give an intuitive 
explanation of why they work. In the following statements "information" always means information 
about the input which has been encoded in the coding space through some isometry. 

1. If all information is perfectly present in B, then all information is absent from B. 

2. If all information is absent from B then all information is perfectly present in B. 

3. If the information about some orthonormal basis (i.e., the type corresponding to this decom- 
position of the identity) is perfectly present in B, then the information about a mutually-unbiased 
basis is absent from B. 

4. If two types of information that are "sufficiently incompatible" are both perfectly present in 
_B, then all information is perfectly present in B. In particular this is so when the two types are 
associated with mutually unbiased bases. 

5. For a code of distance 5 all information is absent from any B if \B\ < S, and all information 
is perfectly present in B if \B\ > n — 5. 

Items 1, 2, 3 and 4 correspond to the No Splitting, Somewhere, Exclusion and Presence theorems 
of |Gri07j ■ which also gives weaker conditions for "sufficiently incompatible." The essential idea 
behind 5 is found in Sec. Ill A of |(;BP97j 0. 
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(a) 



Figure 6.2: (a) The graph state for the [[5,1,3]]d code; (b) The graph state for Steane [[7, 1,3]]2 
code 




1 3 5 2p-l 



2 4 6 2p 
(b) 



Figure 6.3: (a) Complete graph (on 6 qudits); (b) Bar graph with n = 2p carriers and p bars 



6.6.2 One encoded qudit 

It was shown in |Rai99| that a [[5, 1, 3]]d code exists for aU D. Here we consider the graph version 
[SWOlj where the coding group is 

C = (Z1Z2Z3Z4Z5) (6.64) 

and the graph state is shown in Fig. I6.2r a). Our formahsm shows that, whatever the value of D, 
there are only two possibilities. When is 1 or 2 C/-^ is the just the group identity, the projector 
P on the coding space, so all information is absent whereas if \B\ is 3, 4 or (obviously) 5, = G, 
so the subsystem B is the output of a perfect quantum channel. To be sure, these results also follow 
from principle 5 in the above list, given that 5 = 3 for this code. 

The Steane [[7, 1, 3]]2 code, a graphical version of which jYCOj has a coding group 

C = {Z3Z5Z7) (6.65) 

for the graph state shown in Fig. 16.21 (b), is more interesting in that while principles 5 ensures that 
all \B\ < 2 = 6 — 1 subsets of carriers contain zero information and all \B\ > 5 = n — S + 1 subsets 
contain all the information, one qubit, it leaves open the question of what happens when \B\ = 3 or 
4. We find that all information is perfectly present when B is {1,2,5}, {1,3,6}, {1,4,7}, {2,3,4}, 
{2, 6, 7}, {4, 5, 6}, or {3, 5, 7} — representing three different symmetries in terms of the graph in the 
figure — and absent for all other cases of \B\ = 3. Therefore all information is absent from the \B\ = 4 
subsets which arc complements of the seven just listed, and perfectly present in all others of size 
\B\ = 4. So far as we know, generalizations of this code to _D > 2 have not been studied. 

2 It is shown in |GBP97j that if noise only affects a certain subset B of the carriers with |B| < 6, then the errors 
can be corrected using the complementary set B. In our notation this is equivalent to saying that all the information 
is in B. 
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Figure 6.4: The graph state of the [[4, 2, 2]]d code 



A simple code in which a specific type of information is singled out is [[n, 1, 1]]_d generated by 

C = (Z1Z2 • • • Zn) (6.66) 

on the complete graph, illustrated in Fig. I6.3r a) for n ~ 6. Whereas all information is (of course) 
present when \B\ ~ n, it turns out that for any subset B with 1 < \B\ < n one has Qq = {XqiZqi), 
i.e., the abelian group consisting of all powers of the operator XiZi on the input qudit. Thus the 
information is "classical," corresponding to that decomposition of the input identity that diagonalizes 
XiZi. The intuitive explanation for this situation is that this XiZi type of information is separately 
copied as an ideal classical channel, see (|6.6p . to each of the carrier qudits, and as a consequence 
other mutually unbiased types of information are ruled out by principle 3. This, of course, is typical 
of "classical" information, which can always be copied. 

A more interesting example in which distinct types of information come into play is the bar 
graph. Fig. 16.31 (b), in which n qudits are divided up into p = n/2 pairs or "bars," and the code is 
generated by 

C = (ZiZ2---Z„). (6.67) 

Let us say that a subset of carriers B has property I if the corresponding subgraph contains at least 
one of bars, and property II if it contains at least one qudit from each of the bars. Then: 

(i) If B has property I but not II, Qq = (A'oi), an abehan group. 

(ii) If B has property II but not I, Qq = {Xq-^Zqx), another abelian group 

(iii) If B has both property I and property II, all information (1 qudit) is perfectly present. 

(iv) When B has neither property I nor II, all information is absent. 

While both (i) and (ii) are "classical" in an appropriate sense and indeed represent an ideal 
classical channel, the two abelian groups do not commute with each other, so the two types of 
information are incompatible, and it is helpful to distinguish them. Case (iii) illustrates principle 4, 
since Xqi and Xq^Zqi (whatever the value of p) correspond to mutually unbiased bases. In case (iv) 
the complement B of B possesses both properties I and II, and therefore contains all the information, 
so its absence from B is an illustration of principle 1. 

6.6.3 Two encoded qudits 

Consider a [[4, 2, 2]]d code based on the graph state shown in Fig. 16.41 whose coding group 

C = (ZiZ2,Z3Z4), (6.68) 

employs two generators of order D, and thus encodes two qudits. Note that while the graph state 
has the symmetry of a square the coding group has a lower symmetry corresponding to the different 
types of nodes employed in the figure. 

Let us begin with the qubit case D = 2. Our analysis shows that when \B\ = 1 all information 
is absent, and thus for > 3 all information is present, consistent with the fact that this code has 
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(5 = 2 [LYGG08] ■ see principle 5. Thus the interesting eases are those in which \B\ = \B\ = 2, for 
which one finds: 



B 


= {l,3},5 = 


{2,4}: 




rB 

= ^0 


= {^aiZaiZQ2, XqiXq2}', 


(6.69) 


B 


= {1,4},B = 


{2,3}: 


rB 

^0 


rB 




(6.70) 


B 


= {1,2}, B = 


{3,4}: 


rB 


rB 

= ^0 


= {NqiZqi, XQ2Z02) ■ 


(6.71) 



In each case the generators commute and thus the subgroup Gq is abelian. Hence the information is 
"classical", and the same type is present both in B and B, not unlike the situation for the complete 
graph considered earlier. However, the three subgroups do not commute with each other, so the 
corresponding types of information are mutually incompatible, a situation similar to what we found 
for the bar graph. 

For Z? > 2 it is again the case that all information is absent when = 1 completely present 
for \B\ > 3. And (|6.69p and (|6.70|) remain correct (with each generator of order D), and these 
subgroups are again abelian. However, when B = {1,2} and B = {3,4}, (|6.7ip must be replaced 
with 

So = {ZoiXq2, Z02), ^(f = (^01' ^01 ■^02)- (6.72) 

In each case the two generators do not commute with each other, so neither subgroup is abelian. 
However, all elements of Go commute with all elements of Gq ■ Also, the two subgroups are isomor- 
phic (interchange subscripts 1 and 2). 

For odd D > 3 one can use for Gq an alternative pair of generators 

Gl^ = {Z^,X,2,Zq2), m:={D + l)/2, (6.73) 

whose order is D and whose commutator is 

(Z™Xo2) Z02 = ^Zo2 {Zi,'lX,2) ■ (6.74) 

This means — see (|6.8p — that Gq , and thus also the (isomorphic) Gq , is isomorphic to the Pauli 
group of a single qudit. Since Gq and Gq commute with each other, it is natural to think of the 
pair as associated with the tensor product of two qudits with the same D. That this is correct can 
be confirmed by explicitly constructing a "pre-encoding" circuit embodying the unitary 

(Fi®F2)^CPrr(Fi®F2), (6.75) 

expressed in terms of the Fourier and CP gates defined in Sec. I6.3.2[ that carries the Pauli groups 
on "pre- input" qudits 1 and 2 onto Gq and Gq , respectively. 

Things become more complicated for even D > 4, where Gq (and also Gq) are no longer isomor- 
phic to the Pauli group of a single qudit. 

6.7 Conclusion 

We have shown that for additive graph codes with a set of n carrier qudits, each of the same 
dimension D, where D is any integer greater than 1, it is possible to give a precise characterization 
of the information from the coding space that is present in an arbitrary subset B of the carriers. 
This information corresponds to a subgroup G^ of a group G, the information group of operators on 
the coding space, that spans the coding space and provides a useful representation of the information 
that it contains. We discuss how what we call a trivial code, essentially a tensor product of qudits 
of (possibly) different dimensions, can be encoded into the coding space in a manner which gives 
one a clear intuitive interpretation of G- The subgroup G^ is then simply the subset of operators in 
G whose trace down to B is nonzero, and the traced-down operators when suitably normalized form 
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a group Q b that is isomorphic to . The information present in those operators in Q that arc not 
in disappears so far as the subsystem B is concerned, as their partial traces are zero. This is the 
central result of our chapter and is illustrated by a number of simple examples in Sec. 16.61 We also 
provide in App. l6.Cl a relatively simple algorithm for finding the elements of . 

These results can be extended to arbitrary qudit stabilizer codes even if they are not graph 
codes, by employing appropriate stabilizer and information groups, as in Sec. 16.41 Here, however, 
the concept of a trivial code, and thus our perspective on the encoding step, may not apply. The 
extension of these ideas, assuming it is even possible, to more general codes, such as nonadditive 
graph codes, remains an open question. 

As shown in App. 16. Dl our formalism can be fitted within the general framework of invariant 
algebras as discussed in jBKKOTal IBKK07bl IBenl IBKNPV08] . The overall conceptual framework 
we use is somewhat different from that found in these references in that we directly address the 
question of what information is present in the subsystem of interest, rather than asking whether 
there exists some recovery operation (the TZ in App. 16. Pp that will map an algebra of operators 
back onto its original space. Thus in our work the operator groups on the coding space and 
Qb on the subsystem are isomorphic but not identical. Hence, even though there is, obviously, a 
close connection between our "group approach" and the "algebraic approach," the algebra of interest 
being generated from the group of operators, further relationships remain to be explored. The fact 
that the arguments in App. 1 6. Dl are not altogether straightforward suggests that the use of groups in 
cases where this is possible may provide a useful supplement, both mathematically and intuitively, to 
other algebraic ideas. In particular the additional structure present in an additive graph code allows 
one to determine in 0{n^ + K^n'^), App. 16. C[ as against 0{K^) for the algorithm presented 
in |BKNPV08] for a preserved matrix algebra, where K is the dimension of the input and output 
Hilbert space. 

6. A Proof of Lemmas 6.1 and 6.2 

Proof of Lemma 16.11 

The operators in pTZ are linearly independent when those in TZ are linearly independent, since 
p is unitary and thus invertible. This establishes (i). For (ii), consider the case where q is the 
identity /. As the collection TZ is linearly independent, there is at most one r £ TZ such that pr 
is a multiple of the identity. If such an r exists, p is of the form e"^?'^^, and since 7?. is a group, 
pTZ = e^'^r~^TZ = e^'^TZ, we have situation (a), with the collection pTZ U TZ linearly dependent. Next 
assume the collection pTZ U 7?. is linearly dependent, which means there are complex numbers {ar} 
and {br}, not all zero, such that 



This is not possible if all the are zero, since this would mean p^^b^r = 0, thus '^^b^r = 
implying b^ ~ Q for every r, since the TZ collection is by assumption linearly independent. Thus at 
least one a^, say as is not zero. Multiply both sides of (|6.76p by on the right and take the trace: 



implying there is at least one r for which Tr[prs~^] ^ 0. But then p is of the form e"^sr~^ = e"^f~^ 
for f = rs~^ G TZ, so we are back to situation (a). Hence the alternative to (a) is (/3): the operators 
in pTZ U TZ are linearly independent. Finally, if q is not the identity /, simply apply the preceding 
argument with p = q^^p in place of p. 




(6.76) 




(6.77) 



Proof of Lemma 16.21 
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Statement (i) is a consequence of the fact that if an invertible operator is in B, so is its inverse, 
and since 5 is a group, gS consists of the inverses of the elements in g^^S. 

Statements (ii) and (iv) follow from a close examination of (j6.55p . Assume both sets on the left 
side arc nonempty. If gsi and hs2 are both in B, so is their product gsihs2 = ghsiS2, where we 
use the fact that g and h commute with every element of S. If, on the other hand, (ghS) D B and 
(gS) nB are nonempty, any element, say ghsi, in the former can be written using a specific element, 
say gs, in the latter, as 

ghsi = igs){hs2) (6.78) 

where S2 = Sis~^ is uniquely determined by this equation, and the fact that both ghsi and gs are 
(by assumption) in B means the same is true of hs2. Thus not only can every element of (ghS) D B 
be written as a product of elements of (gS) D B, but there is a one-to-one correspondence between 
(ghS) n B and (gS) D B, which must therefore be of equal size. A similar argument shows that 
[ghS] n B and [hS) n B are of the same size. This establishes both (ii) and (iv). 

As for (iii), use the fact that the cosets gS and hS are either identical or have no elements 
in common, so the same is true of their intersections with B. If gS and hS have no elements in 
common. Lemma 16.11 with TZ = S tells us that either gS = e"^(/i5) for some nonzero </>, in which 
case mgS) n 6] = e''^S[(/i5) n B] is distinct from S[(/i5) n K], or else the collection (gS) U (hS) is 
linearly independent, which means that its intersection with B shares this property and the operators 
T,[{gS) n B] and T,[{hS) n B] are linearly independent. 

6.B Proof of Theorem 6.4 

The proof of Theorem 16.41 makes use of the following: 

Lemma 6.5. Let g = PgP be an information operator in Q with spectral decomposition 

m—l 

9=Y.^iJo^ (6.79) 

where the mutually orthogonal projectors Jj sum to P. Then each projector Jj can be written as a 
polynomial in g with g^ = P: 

m— 1 
fc=0 

Proof. The proof consists in noting that 

m—l m—l 

j=o j=0 

is a linear equation in the Jj with /3kj ~ Xj a.n m x m Vandermonde matrix whose determinant is 
UjykiP] - fJ-k) (see p. 29 of |HJ99]). As the fij are distinct the matrix Pkj has an inverse ajk- □ 

To prove (i) of Theorem 16.41 first assume that g is in . Since is a group with identity 
P, this means that all powers of g, including g^ = P, are also in . Consequently, the projectors 
entering the spectral decomposition (|6.79p of g satisfy 

N-'Tts[Jj] TvsiJk] = TrsiJjJk] = '5,fcTr5[J,], (6.82) 

with the first equality obtained by expanding Jj and Jk in powers of g, (|6.80p . and using (|6.56p along 
with the linearity of the partial trace. This orthogonality of the partial traces of different projectors. 
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see (|6.3p . implies that the J^{g) type of information is perfectly present in B. Conversely, if the J"{g) 
type of information is perfectly present in B then the partial traces down to B of the different Jj , 
which cannot be zero, are mutually orthogonal and thus linearly independent. Therefore by (j6.79p . 
Tr^[^] cannot be zero, and g is in . 

The prove (ii) note that g'^ absent from for 1 < k < D means that Tr5[g'=] for these 
values of k, and thus by taking the partial trace of both sides of (|6.80[) and using (|6.57|1 . 

Ttb[JA= Na.oPB- (6.83) 

Since these partial traces are identical up to a multiplicative constant there is no information of the 
J{g) type in B. For the converse, if there is no J^{g) information in B then there is also no J{g^): 
Jiff^), etc. information in S, since the projectors which arise in the spectral decomposition of g*^ 
are already in the spectral decomposition of see (|6.81|) . Consequently, by (i), these g'^ must be 
absent from . 

To prove (iii) , note that if all information is perfectly present in B this means that for every g G G 
the J^{g) information is present in B, and therefore, by (i), g € , so G = G^ ■ For the converse, 
let Qi and Q2 be two orthogonal but otherwise arbitrary projection operators on subspaces of the 
coding space He- Because the elements of the information group G form a basis for the set of linear 
operators on He, see comments at the end of Sec. I6.4.3[ Qi and Q2 can both be written as sums of 
elements g in G, and the same argument that was employed in (|6.82p shows that the orthogonality 
of Qi and Q2 implies the orthogonality of Trg[(5i] and Tr^[(52]- 

To prove (iv), note that if G^ consists entirely of scalar multiples of P, the partial trace down 
to B of any projector Q on a subspace of He, since it can be written as a linear combination 
of the partial traces of the g in G, most of which vanish, will be some multiple of Pb, and thus 
all information is absent from B. Conversely, if G^ contains a g which is not proportional to P 
the corresponding J7{g) type of information will be present in B by (i), so it is not true that all 
information is absent from B, a contradiction. 



6.C Algorithm for finding Q 

Here we present an algorithm for determining the subset information group G^ by finding the 
elements g of G whose partial trace down to B is nonzero. If two or more elements differ only by 
a phase it is obviously only necessary to check one of them. For what follows it is helpful to adopt 
the abbreviation 

^(x|z) _ ^x^z (g 84) 

with (x|z) an ri-tuple row vector pair, and thus a 2n-tuple of integers between and D—1. Arithmetic 
operations in the following analysis are assumed to be mod D. 

First consider the trivial code on the trivial graph, Sec. l6.4.2| with information group Gq consist- 
ing of elements of the form go = goPo, see (|6.48|) . with go = £;(^o|Zo) gome element of Wq = tCq), 
and 

Po = \So\-' J2 (6.85) 
xeA'o 

where Xq denotes the collection of n-tuples that enter the stabilizer Sq, (|6.37p . By choosing xq and 
Zq to be of the form 

xo = (a,6,-.-efc,o,o,...o), 

zo = (CiTii, C2TO2, . . . Cfc"ifc,0,0, . . .0), (6.86) 

using integers in the range 

0<ij<{dj-l), 0<Q<idj-l), (6.87) 
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we obtain a single representative go — ijC^ol^o) fQj. each coset goSo in W/Sq- The corresponding 
information operator, which depends only on the coset, is 

go = £;(^"l^«)Po = \So\-' (6.88) 

where the addition of x and Xg is component-wise mod D, and Zgx denotes the scalar product of Zq 
and X mod D (multiply corresponding components and take the sum mod D). 

Elements of the information group of the nontrivial code of interest to us are then of the 
form 

g = {UW)goiUW)^ 
= \So\-^ (6.89) 

where we use the fact that because the conjugating operator UW, ()6.36p . is a Clifford operator there 
is a 2n X 2n matrix Q over Zj^, representing a symplectic automorphism [HDMOSj . such that 

(t/VF)i;(^l^)(C/VP^)t (6.90) 

with (x.\z)Q the 2ri-tuple, interpreted as an n-tuple pair, obtained by multiplying (x|z) on the right 
by Q, and i/(x, z) an integer whose value does not concern us. The explicit form of Q can be worked 
out by means of the encoding procedure presented in Sec. 16.4.21 using Tables EH] and [621 

The operators appearing in the sum on the right side of (j6.89p are linearly independent Pauli 
products, since Q is nonsingular. The trace down to B of such a product is nonzero if and only if 
its base is in B, and when nonzero the result after the trace is essentially the same operator: see 
(j6.53p and the associated discussion. Consequently gs = N~^TTg[g] is nonzero if and only if the 
trace down to B of at least one operator on the right side of (|6.89p is nonzero. A useful test takes 
the form 

Trg[£;(^l^)] (x|z)J = 0, (6.91) 

where is the zero row vector, and J is a diagonal 2n x 2n matrix with 1 at the diagonal positions 
j and 2j whenever qudit j belongs to B, and elsewhere. Therefore the g associated with xq and 
Zq through (|6.88p and (|6.89|) is a member of if and only if there is at least one x G Xq such that 

(x + xo|zo)QJ = or (x|0)gj = -(xo|zo)QJ. (6.92) 

The X that belong to Xq are characterized by the equation 

xA/ = 0, (6.93) 

where M is an n x fc matrix that is everywhere except for Mjj ~ mj for 1 < j < k, using the rrij 
that appear in (|6.28p . Consequently, instead of asking whether (|6.92p has a solution x belonging to 
Xq one can just as well ask if there is any solution to the pair (|6.92p and (|6.93p . or equivalently to 
the equation 

xT = uo (6.94) 

where T is an n x (2n + k) matrix whose first 2n columns consist of the top half of the matrix QJ, 
(upper n elements of each column), and whose last k columns are the matrix M in (j6.93p : while 
Uq is a row vector whose first 2ti elements are —{xo\zo)QJ and last k elements are 0. Deciding if 
(j6.94p has a solution x becomes straightforward once one has transformed T to Smith normal form, 
including determining the associated invertible matrices, see (j6.32p . As this needs to be done just 
once for a given additive code and a given subset B, the complexity of the algorithm for finding 
is 0{vP) for finding the Smith form plus 0{n^K^) for testing the elements of G once the Smith 
form is available. By using the group property of one can construct a faster algorithm, but that 
is beyond the scope of this Dissertation. 
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6.D Correctable *-algebra 

The counterpart in |BKK07b] of our notion of information perfectly present at the output of a 
quantum channel, see Sec. 16. 2[ is that of a correctable * -algebra A of operators acting on a Hilbert 
space. The * (sometimes denoted C*) means that A, as well as being an algebra of operators in the 
usual sense, contains whenever it contains a. Let the channel superoperator £ be represented by 
Kraus operators, 

£ip)^Y.^iP^h (6.95) 

j 

satisfying the usual closure condition ^j^j ^ ^^^^ ^'^^ P he a projector onto some subspace 
PTi of the Hilbert space H. Then a *-algebra A is defined in |BKK07b] to be correctable for £ on 
states in P'H provided a = PaP for every a in A, and there exists a superoperator TZ (the recovery 
operation in an error correction scheme) whose domain is the range of £, whose range is C{'H)^ and 
such that 

P[{no£)\a)]P ^ PaP (6.96) 
for all a ^ A. Here the dagger denotes the adjoint of the superoperator in the sense that 

Tr[6((7^of)(c))] =Tr[((7eo£)t(6))c] (6.97) 

for any b and c in C{7i). In [BKK07b] . see Theorem 9 and Corollary 10, it is shown that any 
correctable algebra in this sense is a subalgebra of (what we call) a maximal correctable algebra 

Am = {ae C{PH) ■.[a,PEjEjP]=0 V t,j}. (6.98) 

We can apply this to our setting described in Sees. 16.4] and 16.51 where P is the projector on the 
coding space He ^-nd £b is the superoperator for the partial trace down to the subset B of carriers, 

£b{p) = Trsip] = E,pE] for p e £(H) (6.99) 

with Kraus operators 

E,:=lB®{j\B, (6.100) 
where \j)B is any orthonormal basis oiHg, so 

EjE,=lB®\i){j\s- (6.101) 

We shall now show that collection of operators in (defined in Thcorcm l6.3|) spans a ^-algebra 
which is correctable for £b on states in P'H ~ He, and is the maximal algebra of this kind, i.e. 
span(^^) = Am- First note that span{Q^) is indeed a *-algebra: every g Q Q is a unitary operator 
and Q contains the adjoint of each of its elements; replacing g with g'^ in (|6.48p yields g^ . Of course 
Ti'siff] = if and only if Trj3[5^] ~ and in addition, a = PaP for a G span(fj^) because g = PgP, 
(I05D . 

By definition Trg[5] ^ for g E , and this means that the partial trace down to B of at least 
one element in the corresponding coset gS, see (j6.48p . must be nonzero. Let h be such an element; 
since it is a Pauli product it must be of the form h = hs ® !§■ As a consequence, 

[g,PElE,P] = [h,PElE,P] - P[h,EjE,]P 

^P[hB® IbJb ® \i) (jls ]P = 0, (6.102) 

where the successive steps are justified as follows. Since g depends only on the coset gS and h 
belongs to this coset, hS = gS and h = Ph = hP ~ g. This means we can move the projector P 
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outside the commutator bracket, and once outside it is obvious that the latter vanishes for every i 
and j. Thus any g in belongs to the maximal Am defined in (|6.98|) . as do all linear combinations 
of the elements in . 

To show that Am is actually spanned by we note that any a belonging to Am can be written 

as 

a = b + c, (6.f03) 

where 6 is a linear combination of elements of and c of elements of G that do not belong to , 
so Trg[c] — Tr^[c^^] = 0. Thus it is the case that 

P{no£B)Hb)P = b, P{TZo£b)\c)P ^c, (6.104) 

where the first follows, see (|6.96p . from the previous argument showing that the span of G^ is a 
subalgebra oi Am, and the second from linearity and the assumption that a belongs to Am- Multiply 
the second equation by and take the trace: 

Tr[c1"c] = Tr [c^P ((7^ o £b)\c)) P] 

= Tr [(7^ o £b{c^)) c] = 0, (6.105) 

where we used the fact that Pc^P = c\ and £b{c^) = Tr^[c^] = 0. Thus c = and any element of 
Am is a linear combination of the operators in G^ ■ 

In conclusion, we have shown for any additive graph code C and any subset of carrier qudits B, 
the *-algebra spanned by operators in G^ is exactly the maximal correctable algebra Am defined in 
(j6.98p . In App. 16. CI we outline an algorithm that enumerates the elements in G^ for any He and 
£b, which in light of the result above is an operator basis of Am- 
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Bipartite equientagled bases 



7.1 Introduction 

We present two different solutions to the problem posed by Karimipour and Memarzadeh in [KM06| 
of constructing an ortlionormal basis of two qudits with the following properties: (i) The basis 
continuously changes from a product basis (every basis state is a product state) to a maximally 
entangled basis (every basis state is maximally entangled), by varying some parameter t, and (ii) 
for a fixed t, all basis states are equally entangled. As mentioned in [KM06| . such a family of bases 
may find applications in various quantum information protocols including quantum cryptography, 
optimal Bell tests, investigation of the enhancement of channel capacity due to entanglement and 
the study of multipartite entanglement. For a more detailed motivation the interested reader may 

consult }KM06| . 

The chapter is organized as follows : In Sec. 17.21 we summarize the main r esuhs of |KM06j and 
then introduce the concept of Gauss sums and some useful related properties. Next we provide 
an explicit parameterization of a family of equientangled bases and we prove that it interpolates 
continuously between a product basis and a maximally entangled basis, for all dimensions. We 
illustrate the behaviour of our solution with explicit examples. In Sec. 17.31 we construct another 
such family using a completely different method based on graph states, describe a simple extension 
of it to multipartite systems, and then illustrate its behaviour with specific examples. Finally in 
Sec. I7.4l we compare the two solutions and make some concluding remarks. 

7.2 Construction based on Gauss sums 
7.2.1 Summary of previous work 

Let us start by summarizing the main results of |KM06j . Consider a bipartite Hilbert space 7i(E)'H, 
where both Hilbert spaces have the same dimension D. The authors first defined an arbitrary 
normalized bipartite state 




(7.1) 



fc=0 



Next for TO, 71 = 0, 1, — 1, they considered the collection of "shifted" states 



lV'o,o) 



D-1 




(7.2) 
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where 

D-l 



X:=^|fc®l)(fc| (7.3) 



fe=0 



is the generahzed Pauh (or shift) operator and denotes addition modulo D. They noted that all 
states have the same value of entropy of entanglement |NC00j given by the von-Neumann entropy 



D-l 



i?(l^rn,„)) = £^(|V'0,0)) = -5I I«fcl'l0gl3l«fcl'' (^-4) 



fc=0 



where the logarithm is taken in base D for normalization reasons so that all maximally entangled 
states have entanglement equal to one regardless of D. 
Demanding the states in (|7.2p be orthonormal yields 



D-l 

(afc)*afc®„ = S„i,o, Vm = 0, . . . , D - 1, (7.5) 



fc=0 



and the authors proved (see their Eqn. (36)) that (|7.5|) is satisfied if and only if the coefEcients Uk 
have the form 



D-l 

i 



^E^"^'^''' (7.6) 



D 

where 9j are arbitrary real parameters and lu = c^'^'f^ is the Z?-th root of unity. 

Therefore the authors found a family of orthonormal states, all having the same Schmidt 
coefficients and hence the same value of entanglement. To ensure it interpolates from a product 
basis to a maximally entangled basis, it is sufficient to find a set of parameters {d^}^SQ^ for which 
the magnitude of aj. is \ak\ = 1/Vd for all k. Then the problem is solved by defining 

ak{t):=^j:e^''°u:''^ (7.7) 

where t e [0, 1] is a real parameter. When t = we have ak ~ Sk^ so the basis states are product 
states and when < = 1, the basis is maximally entangled by assumption. We also observe there is a 
continuous variation in between these two extremes as a function of t. 

Karimipour and Mcmarzadeh considered the existence of such a set {^^ }^^^ in arbitrary dimen- 
sions (see the last paragraph of Sec. V in jKM06j ) . They found particular solutions for D ^ 5, but 
did not find a general solution for arbitrary D. 

7.2.2 Quadratic Gauss Sums 

We now define the basic mathematical tools we will make use of in the rest of this section. The most 
important concept is that of a quadratic Gauss sum, defined below. 

Quadratic Gauss Sums. Let p, m be positive integers. The quadratic Gauss sum is defined as 

p-i 

^^,2^i/m/p^ (7.8) 

3=0 

The quadratic Gauss sums satisfy a reciprocity relation known as 
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Landsberg-Schaar Identity. Letp,m be positive integers. Then 

p-l ^ 2m-l 

The quadratic Gauss sums can be generalized as follows. 

Generalized Quadratic Gauss Sums. Let p, to, n be positive integers. The generalized quadratic 
Gauss sum is defined as 

^^g27ri(i^m+jn)/p_ (7-10) 

i=o 

Finally the following reciprocity formula for generalized Gauss sums holds. 

Reciprocity Formula for Generalized Quadratic Gauss Sums. Letp,m,n be positive integers 
.such that mp ^ and mp + n is even. Then 

p-l m-l 

g7riO"^m+j"n)/p _ ^Tri(mp-n^) / imp Q-TrKfp+jn)/m -|^|\ 

,/p ^ Jm ^ ' \ ■ ) 

i=0 V 

The definitions of the Gauss sums (|7.8p and (|7.10p as well as the Landsberg-Schaar's identity 
(|7.9[) can be found in standard number theory books [HW08[ IEW051 INatOO] . The reciprocity formula 
for the generalized quadratic Gauss sum is not as well-known, and can be found in |BE81) . 

7.2.3 Explicit Solution 

We now show that a family of equientangled bases that interpolates continuously between the product 
basis and the maximally entangled basis exists for all dimensions Z?, as summarized by the following 
Theorem. 

Theorem 7.1. The collection of normalized states 

D-l 

= "-kWlk © to)|A; © TO © n), (7.12) 

TO, 71 = 0, . . . , D — 1, 

indexed by a real parameter t G [0, 1] with 

3=0 

with the particular choice of 



TTp/D if D is even 
2np/D if D is odd, 



defines a family of equientangled bases that continuously interpolates between a product basis att = Q 
and a maximally entangled basis at t ~ I. 

That (|7.12p defines a family of equientangled bases that consists of a product basis at i = 
follows directly from the remarks of Sec. 17.11 afc(O) = Sk.o. Next note that a continuous variation of 
t in the interval [0, 1] corresponds to a continuous variation of the Schmidt coefficients of the states 
in the basis. The latter implies that no matter which measure one uses to quantify the entanglement. 
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the measure will vary continuously with t, since any pure state entanglement measure depends only 
on the Schmidt coefficients of the state |Vid99j . 

The only thing left to show is that the basis states in Theorem [73] are maximally entangled when 
t — 1, or, equivalently, that |afe(l)| = 1/Vd for all k. We prove this by explicitly evaluating the 
value of afe(l) in the following Lemma. 

Lemma 7.2. Let ak(t) and be as defined by Theorem \7.1\ Then for all k 

g7ri/4 ( o;"'"'^/^, if D is even 

Ofcfl) = ^X< 1.2/,/, :2k + D\ „^ ,, . (7.15) 

^' Vd I c<j-^ (i^^ j , if D IS odd ^ ' 

Lemma [7.21 implies at once that |aA;(l)| = I/a/D, and therefore proves Theorem 17. II 

Proof, (of Lemma [7.2p Note first that the expression for 0^,(1) in (|7.13p with defined in (|7.14p 
resembles the generalized quadratic Gauss sum (|7.10p . We will use the reciprocity formula (|7.1ip to 
prove Lemma [721 There are two cases to be considered: Even D and odd D. 

Even D. Note that one can rewrite afc(l) in (|7.13p with 6*° defined in (|7.14p as 



3=0 i=o 

Applying the reciprocity formula ()7.1ip to last term in ()7.16p with m — l,n ~ "^.k^p ~ D (noting 
that mp + n = D + 2k is even) yields 

n,(^] - _J_„^i{D-4k~)/4D _ 2 -7Tik~/D 

ak{l)- r^C - [ L) C 



g7ri/4 2 

= ^^w-'^ since = 1 for even D. (7.17) 

Odd D. The proof is essentially the same as in the even D case, but we explicitly write it below 
for the sake of completeness. Using a similar argument we rewrite afe(l) in ()7.13p with 9j defined in 
dm as 

Applying again the reciprocity formula (|7.1ip to last term in (|7.18p with m = 2,n = 2k, p = D 
(noting that mp + n = 2D + 2k is even) yields 

where we used (— 1)2*^+^ = _(^[2k+D^ since 2k + D is odd. This concludes the proof of Lemma 17^ 
and implicitly of Theorem 17.11 □ 
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D=5 

|a^(t)l 




0.0 0.2 0.4 0.6 0.8 1.0 ' 



Figure 7.1: The variation of |afc(t)| with t for _D = 5. Note how at t = all coefficients but one are 
zero, and how at t = 1 all coefficients are equal in magnitude to l/\/5, with a continuous variation 
in between. The dashed line represents the l/Vb constant function. 



7.2.4 Examples 

In this section we present some examples that illustrate the behaviour of the solution we provided 
in Theorem 17. 11 for various dimensions. First we consider D = 5 and we plot the absolute values of 
the ak (t) coefficients as a function of t in Fig. 17.11 It is easy to see that indeed the basis interpolates 
between a product basis and a maximally entangled one in a continuous manner. We observe that 
all coefficients are non-zero for i > and we believe that this is probably also the case for all odd 
D's. 

In Fig. 17.21 we perform the same analysis as above, but now for D = 8. We observed that some 
coefficients vanish for some values of t, which seems to be true in general for even D. 

In Fig. 17.31 we plot the entropy of entanglement of the states in the basis as a function of t 
for dimensions D = 2, 3, 5, 8 and 100. We see how the entanglement varies continuously but not 
monotonically between and 1. 

Finally in Fig. I7.4l we display a parametric plot of the variation of the second Schmidt coefficient 
ai{t) in the complex plane as t is varied from to 1 for D = 51, so that the reader can get an idea 
of how the coefficients defined in (|7.13p look in general. The other coefficients at look similar. 



7.3 Construction based on Graph States 
7.3.1 Explicit solution 

We provide below another solution to the problem that uses qudit graph states. Again having in 
mind a bipartite Hilbert space 'H®'H, both local spaces having dimension D, we define a one-qudit 
state ^ 

It is easy to see that the collection of D states 

|m):=Z'"i+), m = 0,...,i:>-l (7.21) 
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D=8 




Figure 7.2: The variation of |afe(i)| with t for D = 8. Again note how at i = all coefficients but 
one are zero, and how at t = 1 all coefficients are equal in magnitude to 1/V8, with a continuous 
variation in between. The dashed line represents the l/VS constant function. 



Normalized 
Entropy of 
Entanglement 




Figure 7.3: The entropy of entanglement as a function of t for various dimensions. Note that the 
variation is not monotonic (except for D — 2), although for large D the oscillations tend to be 
smoothed out. 
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D=51 
Im[ai(t)] 




Figure 7.4: Parametric plot of ai{t) in the complex plane as t is varied from to 1. Note that 
ai(0) = and ai(l) = J"^~j^ e'^'^^~2^), the value provided by Lemma [^21 The starting point t ~ 
and the ending point t ~ 1 are marked by solid disks. 



defines an orthonormal basis of Ji (also known as the Fourier basis), {m\n) = 8mm where 



(7.22) 



fc=0 

with bj ~ c^^^/^ being the D-th root of unity. It then follows at once that the collection of states 

|m)|n) = (Z"'® Z")|+)|+), m,n^{),...,D-l (7.23) 

defines an orthonormal product basis of the bipartite Hilbert space T-L®T-L. 
Next we define the generalized controUed-Phase gate as 

D-l D-l 

C:=Y, \k) {k\ ®Z^=Y, {j] ® \k) {k\ (7.24) 

fc=0 

and note that C is a unitary operator that commutes with ® Z" , for all m, n = 0, . . . , -D — 1. 
The state ^ 

\G):=C\+)\+) = ^ oj^'^mk) (7.25) 
j.k={) 

is an example of a two-qudit graph state and it is not hard to see that \G) is maximally entangled. 
Then the collection of states 

(Z" Z'')\G) = (Z™ ® Z")C|+)|+) 
= C(Z'"®Z")|+)|+), m,n = 0,...,i?-l (7.26) 

defines an orthonormal basis of the bipartite Hilbert space H^H, which we call a graph basis. Since 
(g) Z" are local unitaries and \G) is a maximally entangled state, then all the other graph basis 
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states must also be maximally entangled. For more details about graph states of arbitrary dimension 
see |HDE+[lLYGG08llGLG10| . 

We now have all the tools to construct a continuous interpolating family of equientanglcd bases, 
as summarized by the Theorem below. 

Theorem 7.3. The collection of normalized states 

\G^At)) = ® (7.27) 
m, n = 0, . . . , D — 1, 

indexed by a real parameter t G [0, 1] where 

C(t)=^..^'^-*b-)(j|®|fc)(fc| (7.28) 

defines a family of equientanglcd bases that continuously interpolates between a product basis att — 
and a maximally entangled basis at t ~ I. 

Proof. Wc make the crucial observation that C{t) commutes with Z"^®Z'^ for all m, n = 0, — 1 

and all t G [0,1] which implies that {|G„i_„(t))}^^^^Q defines an orthonormal basis since it differs 
from the orthonormal basis in (|7.23|) only by the unitary operator G{t). All states in the basis 
are equally entangled, and moreover, share the same set of Schmidt coefficients since any two basis 
states are equivalent up to local unitaries of the form Z™ ® Z". 

Finally note that C(t = 0) = / (g) / and C{t = 1) = C (defined in (fTM]) ). and therefore at t = 
the basis is product, see ()7.23|) . and at i = 1 the basis is maximally entangled, see (|7.26|) . The 
operator C(t) can be viewed as a controUed-Phase gate whose "entangling strength" can be tuned 
continuously. The Schmidt coefficients of the states in the basis vary continuously with t and hence 
the entanglement also varies continuously with t, regardless of which entanglement measure one uses 
(see the remarks following Theorem 17. ip . □ 

Our construction above can be expressed in the framework described in the last two paragraphs 
of Sec. II of |KM06j . by setting = Z" and K = Z^". 

Next we prove that the Schmidt coefficients of the basis states in Theorem 17.31 are all non-zero 
for any t > 0, so all bases consist of full Schmidt rank states whenever t > 0. 

Lemma 7.4. The equientanglcd family of bases {|Gm,n(t))}m7i=o defined in Theorem \7.3\ consists 
of full Schmidt rank states, for any < t ^ 1. 

Proof. We will show that the product of the Schmidt coefficients is always non-zero, which implies 
that no Schmidt coefficient can be zero, whenever < t ^ 1. 
Let 

D-l 

1^) = E ^^i^im (7-29) 

be an arbitrary normalized pure state in a bipartite Hilbert space 'H®'H and let {A^} denote the set 
of Schmidt coefficients of \tp) satisfying Afc = 1; note that they are equal to the squares of the 
singular values of the coefficient matrix in (|7.29p . The product of the squares of the singular values 
is just the product of the eigenvalues of the latter product being equal to det(51il''') — \ dct(J7)p, 
so we conclude that 

D-l 

[] Afc = |det(f))p. (7.30) 

fc=0 



7. BIPARTITE EQUIENTAGLED BASES 



96 



The states \Gm.n{t)) of Theorem 17.31 share the same set of Schmidt coefficients (they are all 
related by local unitaries) so it suffices to show that the product of the Schmidt coefficients is 
non-zero only for the state \Go,o{t))- Recall that 

jkt 

\GoAt)) = Cit)\+)\+) = ^b->|fc). (7.31) 

Expressing the coefficients u!-'^*/D as a matrix il{t), one can easily see that D • f2(i) is a D x D 
Vandermonde matrix whose determinant is 

det [D ■ Q{t)] = Yl i^^' ~ ^^') = n W^'""^' - l] (7.32) 
]>k j>k 

(see p. 29 of [HJ99| for more details on Vandermonde matrices). For a given t, the product above is 
zero if and only if at least one term is zero, i.e. there must exist integers j, k, with ^ k < j ^ D—1, 
such that 

(i - k)t = nD ^t = n -, (7.33) 

j — k 

for some positive integer n ^ 0. Note that < j — fc ^ D — 1, so D/{j ~ k) > 1 and the above 
equation can never be satisfied for < i ^ 1. We have therefore proved that det[r2(t)] ^ for 
< t ^ 1, which, in the light of (|7.30p . is equivalent to saying that the product of the Schmidt 
coefficients is non-zero for < t ^ 1, and this concludes the proof of the Lemma. □ 

For this family of equientangled bases, we do not have an analytic expression for the Schmidt 
coefficients nor the entropy of entanglement for general D though they can be easily found by 
numerically diagonalizing the coefficient matrix n{t)n{ty . Having said that, we derived a simple 
analytic expression for the product of all Schmidt coefficients, see ()7.30|) and (|7.32p . which is simply 
related to an entanglement monotone called G-concurrence (first introduced in jGouOSj ) which is 
defined for a pure bipartite state (|7.29p in terms of its Schmidt coefficients {A^} as 

/D-i 

Ccim) ^dIUxA ^D\ det(0)|2/^ (7.34) 

\k=0 ) 

where X]fc -^fe = 1- ^^'^ G-concurrcnce is zero whenever at least one Schmidt coefficient is zero and 
is equal to one if and only if the state is maximally entangled. Unlike the entropy of entanglement, 
we are able to show two analytical results that are true for all D expressed in Lemmas 17.51 and 17.61 

Lemma 7.5. The G-concurrence of the equientangled basis states, {\G„i,n{t))}^'^n=o defined in 
Theorem \7.3\ is 

2D-1 D-i 

Gait) = — ^ Yl [sin^iT^rt/D)]^"^''^^" for all m,n,D. (7.35) 



Proof. Every basis state has the same G-concurrence since they all share the same set of Schmidt 
coefficients (recall that they differ only by local unitaries). Invoking the definition of G-concurrence, 
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wc have 



Cait) = D\det{n{t))\^/° ^ D 

2/D 



— det{D-n{t)) 



2/D 



1 



j>k 
D-1 r 



D-1 



in 



r=l 
D-1 



2-2- 



1 -Q 2iD~r)/D _ cosi2nrt/D)]^ 



D 



r=l 
D-1 



± Y[ 22(^--)/^ [sin^nrt/D)] — 



r=l 
D-1 D-1 



D 



(7.36) 
(7.37) 

(7.38) 
(7.39) 
(7.40) 
(7.41) 



where in (|7.36p we used the fact that dct(cM) = det(M) ior a. D x D arbitrary matrix AI and 
an arbitrary constant c. The first equahty in (j7.37p follows at once from (|7.32p . whereas the second 
equality in (j7.37p follows from a simple counting argument in which one replaces j — fc by r, making 
sure that the different pairs (j, k), j > k, that give rise to the same r are counted; for a given r there 
are D — r such pairs. □ 



It turns out that the G-concurrence has the following nice property: 

in (j7.35p is strictly increasing 



Lemma 7.6. The G-concurrence of the basis states {\Gm.n{t))} 
in the open interval t G (0, 1), for all dimensions D. 



D-1 
m.n— 



Proof. We prove this by showing that the first derivative of the Gait) with respect to t is strictly 
positive. First note that since Gait) > is a positive function in the interval t £ (0, 1). This means 
showing j^Goit) > is equivalent to showing that ^ log CG(i) > 0. The derivative of the logarithm 
of (TMTI) is 

log Gait) = ;^ E '^(^ - cotinrt/D), (7.42) 



dt 



where cot(-) denotes the cotangent function. We differentiate again to get 



D-1 



dt^ 



log Gait) 



^ 1 
2^ r^D-r) 



D3 



r=l 



sin^ {irrt/D) 



(7.43) 



Note that the right hand side of (|7.43|) is strictly negative whenever t > 0, which implies that the 
first derivative of the logarithm (|7.42p is a strictly decreasing function of t and hence achieving its 
minimum value at t = 1, which is given by 



D-1 



^ log Gait) Ui= E '^(^ - ^ot{nr/D) = 0, 



(7.44) 



where the last equality follows from symmetry considerations (terms cancel one by one). We have 
shown ^ log Gait) > 4^ 'Si^cit) > and therefore we conclude that the G-concurrence is strictly 
increasing for i e (0, 1). □ 
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D=5 



[Schmidt(t)]"2 
1.0 




Figure 7.5: The square roots of the Schmidt coefficients as functions of t for D — 5. Note how 
at t = all coefficients but one arc zero, and how at t = I all coefficients are equal in magnitude 
to l/\/5, with a continuous variation in between. The dashed line represents the l/Vb constant 
function. 



7.3.2 Extension to multipartite systems 



The construction presented in Theorem 17.31 can be easily generalized to multipartite systems of 
arbitrary dimension. The concept of maximally entangled states is not defined for three parties 
or more, so in this case, the family continuously interpolates between a product basis and a qudit 
graph basis. It is still true that for a fixed t, all basis states constructed this way have the same 
entanglement (as quantified by any entanglement measure) since they only differ by local unitaries. 

As a specific example, consider the tripartite GHZ state (|000) + |lll))/-\/2- This state is a 
stabilizer state jNCOOj and therefore is local- unitary equivalent |Schbj to a graph state \G) = 
(Ci2C23Ci3)|-|-)i|-|-)2|+)3, where the subscripts on C indicates which pair of qubits the C gate 
is applied to. By varying the "strength" of the controUed-Phase gate, one can now construct a fam- 
ily of equally entangled basis for the Hilbert space of 3 qubits that continuously interpolates between 
a product basis and the GHZ-likc graph basis. This GHZ construction can be easily generalized to 
higher dimensions and also to n parties by using the complete graph given by 



(7.45) 



= 1 ]>i 



where |+) is defined in (|7.20p and C{t) is defined in (j7.28p . Finally note that this construction 
works for any graph state of any dimension, and not just for GHZ-like graph states. Such bases with 
tunable entanglement may be of use in the study of multipartite entanglement. 



7.3.3 Examples 

In this subsection we perform a similar analysis as the one in Sec. 17.2.41 so that one can easily 
compare the behaviour of both solutions. 

We consider again a D = 5 example, for which we plot in Fig. 17.51 the square root of the Schmidt 
coefficients as functions of t. It is easy to see that indeed the basis interpolates between a product 
basis and a maximally entangled one in a continuous manner. As proven in Lemma l7.41 all Schmidt 
coefficients are non-zero for t > 0. 
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In Fig. 17.61 we plot the same quantities for Z? = 8. 

D=8 

[Schmidtft)!'/^ 




Figure 7.6: The square roots of the Sehmidt coeffieients as funetions of i for D = 8. Again note how 
at t — all coefficients but one are zero, and how at t = I all coefficients are equal in magnitude 
to l/VS, with a continuous variation in between. The dashed line represents the l/VS constant 
function. 
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Figure 7.7: The entropy of entanglement as function of t for various dimensions. Note that the 
variation seems to be monotonically increasing for all Z?, a statement we did not prove. 



Observe how in both examples above the variation of the Sehmidt coefficients is not oscillatory, 
as in the examples of Sec. 17.2.41 

In Fig. 17.71 we plot the entropy of entanglement of the basis states as a function of t for dimen- 
sions D = 2, 3, 5, 8 and 100. We observe that the entropy of entanglement varies continuously and 
monotonically between and 1. It is not known if the entropy of entanglement is always strictly 
increasing for all D although we verified this by visual inspection for all D ^ 10. In Fig. 17. 8[ we 
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Figure 7.8: The G-concurrence as function of t for various dimensions. The variation is strictly 
increasing in t for all D as shown in Lemma 17.61 



plot the G-concurrence for the same dimensions. We see how the curves are strictly increasing and 
this is true for all D as proven in Lemma 17.61 

7.4 Conclusion 

We have solved the problem posed in |KM06j by providing two families of equicntangled bases for 
two identical qudits for arbitrary dimension D. The construction of the first solution is based on 
quadratic Gauss sums and follows along the lines of jKM06] , whereas the second family is constructed 
using a different method based on qudit graph states. 

The first solution based on quadratic Gauss sums has an explicit analytic expression for the 
Schmidt coefficients that is easy to evaluate since they are just sums with D terms (see (|7.13|) and 
(|7.14p '). However some Schmidt coefficients can be zero and the entropy of entanglement of the 
states in the basis varies non-monotonically with t ior D > 2. 

The second solution based on graph states consists entirely of full Schmidt rank states for any 
< t ^ 1 that seem to have an entropy of entanglement that is strictly increasing as t increases. 
Unfortunately we did not find a simple analytic expression for the Schmidt coefficients, but they 
can be computed numerically without much difficulty. We found a simple analytic expression for 
another pure state entanglement measure, the G-concurrence, which we proved is strictly increasing 
as t increases. Finally we remark that one can extend this construction to equally entangled bases 
of more than two parties that interpolate continuously between a product basis and a graph basis 
even if the concept of maximally entangled states is not defined for more than two parties. This 
construction may be of interest in studying multipartite entanglement. 
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